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Introduction 


This  report,  the  thirty-fifth  in  a  series  of  quarterly  progress 
reports  issued  by  the  Research  Laboratory  of  Electronics,  contains 
a  review  of  the  research  activities  of  the  Laboratory  for  the  three- 
month  period  ending  August  31,  1954.  Inasmuch  as  this  is  a  pre¬ 
liminary  report,  no  results  should  be  considered  final. 


I.  PHYSICAL  ELECTRONICS 


Prof.  W.  B.  Nottingham  D.  H.  Dickey 

S.  Aisenberg  J.  D.  Hobbs 

W.  M.  Bullis 


J.  M.  Houston 
W.  J.  Lange 
L.  E.  Sprague 


A.  ELECTRON  EMISSION  PROBLEMS 

1.  A  Redetermination  of  the  Crystallographic  Variation  of  Electron  Field 
Emission  from  Tungsten 

This  work  has  been  completed  and  will  be  presented  in  Technical  Report  No.  286. 

A  correction  is  necessary  to  the  discussion  of  experimental  results  given  in  the 
Quarterly  Progress  Report  of  July  15,  1954.  The  discontinuity  in  slope  at  approximately 
800°K  noted  in  the  preliminary  emission  vs.  temperature  measurements  was  found  to 
be  nonexistent  when  a  more  accurate  temperature  scale  was  used.  The  remainder  of 
the  discussion  is  still  correct. 

J.  M.  Houston 


2.  Growth  of  Single  Crystals  of  Tantalum 


KOVAR 


Some  rather  large  single  crystals  of  tantalum  have  been  grown  by  resistance  heating 
and  a  moving  temperature  gradient.  The  moving  temperature  gradient  is  achieved  by 
using  mercury  as  the  electrode,  as  shown  in  Fig.  1-1.  The  tantalum  ribbon  is  fastened 
to  the  upper  electrode.  The  lower  end  of  the  ribbon  is  welded  to  a  5-cm  length  of 
120-mil  tungsten  rod  which  is  submerged  in  the  mercury  and  serves  to  weight  the  ribbon 
sufficiently  to  overcome  surface  tension  and  maintain  the  ribbon  vertical. 

With  the  metal  high -vacuum  valve  closed,  enough  air  is  admitted  to  the  mercury 
reservoir  through  the  stopcock  to  push  mercury  into  the  tube  containing  the  ribbon. 

The  metal  valve  is  used  as  an  adjustable 
slow  leak  to  the  forepump,  which  slowly 
exhausts  the  air  from  the  mercury  reser¬ 
voir,  thus  lowering  the  level  in  the  tube  at 
the  desired  rate. 

The  tantalum  ribbon  is  1  mil  thick 
and  2-3  mm  wide.  We  have  obtained  the 
largest  crystals  after  a  3  per  cent  cold 
work,  as  first  reported  by  Mrowca  (1).  We 
have  cold-worked  the  tantalum  by  cold 
rolling  rather  than  cold  stretch.  The 
percentage  of  cold  roll  is  then  measured 
by  the  elongation  of  the  ribbon.  With  a 


Fig.  1-1 

Tantalum  crystal-growing  setup. 
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temperature  of  about  2400 °K  and  a  mercury  lowering  rate  of  approximately  8  mm/hr, 
at  least  one  single  crystal  of  approximately  1  cm  in  length  is  obtained  on  nearly  every 
run  of  6  to  8  hours. 

W.  J.  Lange 

References 

1.  B.  A.  Mrowca,  J.  Appl.  Phys.  14,  684  (1943). 

B.  EXPERIMENTAL  STUDIES 

1.  Ionization  Gauge  Control  Circuit  and  Ionization  Gauge 

The  ionization  gauge  control  circuit  has  been  constructed,  modified,  and  used.  The 
modifications  include  an  additional  negative  feedback  loop  to  reduce  the  input  impedance, 
the  time  constant,  and  the  zero  drift  (caused  by  line  voltage  variation).  A  diode  limiter 
was  also  included  to  protect  the  electron  current  meter  against  accidental  shorts.  The 
circuit  is  still  free  from  batteries.  When  the  drawings  and  operating  instructions  have 
been  completed,  a  more  detailed  description  of  the  circuit  will  be  given. 

In  order  to  test  the  operation  of  the  control  circuit,  it  was  used  to  investigate  the 
variation  with  geometry  of  the  gauge  sensitivity  of  the  screened  Bayard -Alpert  type  of 
ionization  gauge.  The  sensitivity  of  the  ion  gauge  is  measured  by:  K  =  l/p  (i+/i  ), 
where  i+  and  i_  are  the  ion  and  electron  currents,  and  p  is  the  pressure.  The  sensitiv¬ 
ity,  K,  depends  upon  the  nature  of  the  gas,  the  electrode  potentials,  the  electrode  geom¬ 
etry,  and  is  not  supposed  to  depend  upon  electron  current.  The  Bayard-Alpert  ion 
gauge  in  our  laboratory  has  two  supposedly  identical  filaments:  one  used  as  the  electron 
emitter,  and  the  other  as  a  flash  filament  and  a  spare.  In  operation,  the  flash  filament 
is  at  electron  collector  potential,  while  the  emitting  filament  is  100  volts  negative  with 
respect  to  the  electron  collector.  Using  the  modified  ion  control  circuit  plus  an  appro¬ 
priate  switching  circuit  to  quickly  interchange  the  two  hot  filaments  electrically,  the 
relative  values  of  the  gauge  constants  for  the  two  filaments  were  measured  for  various 
identical  electron  currents.  Not  only  was  a  discrepancy  found  between  the  two  gauge 
constants,  where  none  was  expected,  but  the  discrepancy  was  a  function  of  the  elec¬ 
tron  current,  and  varied  from  approximately  20  per  cent  at  10  pa  to  10  per  cent  at 
10  ma  for  a  sealed-off  ion  gauge.  Incomplete  measurements  on  another  gauge  indi¬ 
cate  that  the  discrepancy  can  be  as  great  as  40  per  cent,  and  may  cause  a  corresponding 
error  in  the  pressure  measurement.  The  discrepancy  was  reduced  at  high  electron 
currents,  probably  because  the  higher  space  charge  influenced  the  electron  orbits  more 
than  the  electrode  geometry  and  made  the  difference  in  geometry  less  important. 

The  results  indicate  that  the  gauge  constant  K  depends  upon  the  electron  current, 
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and  that  the  ion  current  is  therefore  not  exactly  a  linear  function  of  electron  current, 
as  is  commonly  supposed.  A  method  has  been  developed  for  measuring  the  dependence 
of  ion  current  upon  electron  current  in  spite  of  the  changes  of  pressure  caused  by 
changes  of  current  to  the  electron  collector. 

S.  Aisenberg 


2.  Infrared  Radiation  Pyrometer 

Preliminary  measurements  taken  with  an  optical  bench  model  of  the  revised  optical 
system  described  in  the  Quarterly  Progress  Report,  July  15,  1954,  indicate  that  the 
sensitivity  and  response  are  substantially  the  same  as  previously  reported. 

A  pilot  model  of  the  complete  revised  optical  system  has  been  constructed;  the 
prism  has  been  procured,  and  when  it  is  prepared  this  system  will  be  aligned  and  further 
tests  made. 

W.  M.  Bullis 
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II.  MICROWAVE  GASEOUS  DISCHARGES 


Prof.  S.  C.  Brown  A.  L.  Gilardini  Helen  R.  Land 

Prof.  W.  P.  Allis  E.  Gordon  J.  J.  McCarthy 

A.  MICROWAVE  BREAKDOWN  IN  HYDROGEN  AT  HIGH  PRESSURES 

In  preparation  for  extending  data  on  the  continuous  wave  breakdown  of  hydrogen  at 
high  pressure  by  the  use  of  pulse  techniques,  a  study  was  made  of  the  effect  of  various 
pulse  widths  and  pulse  repetition  frequencies  on  the  breakdown  of  hydrogen.  This  was 
done  to  determine  what  combinations  of  pulse  widths  and  pulse  repetition  frequencies 
would  give  values  of  breakdown  electric  field  equal  to  that  obtained  for  the  continuous 
wave.  The  heavy  lines  in  Fig.  II- 1  represent  the  data  taken  of  breakdown  for  various 
values  of  pulse  repetition  frequency  as  a  function  of  pulse  width  for  a  given  value  of 
pressure  of  hydrogen.  In  all  instances,  as  the  pulse  width  increases,  the  value  of  break¬ 
down  electric  field  approaches  that  for  continuous  wave  breakdown.  This  effect,  of 
course,  occurs  at  smaller  values  of  pulse  width  for  large  repetition  rates  than  for  small 
repetition  rates. 

Because  the  plate  dissipation  of  the  magnetron  is  limited,  the  increased  peak  output 
of  the  magnetron  must  be  accompanied  by  an  appropriate  change  in  the  duty  cycle  of  the 
pulse.  If  the  peak  power  is  to  be  twice  the  continuous  wave  power,  the  duty  cycle  must 
not  exceed  0.5.  On  the  curves  of  Fig.  II- 1  the  pulse  duty  cycle  is  plotted  as  a  dotted 
line.  From  these  curves  it  can  be  seen  that  there  are  values  of  duty  cycles  which  make 
it  impossible  to  obtain  the  combinations  of  pulse  width  and  pulse  repetition  frequency 
which  would  give  an  equivalent  continuous  wave  breakdown.  However,  for  duty  cycles 
greater  than  0.  1  there  are  combinations  of  pulse  width  and  repetition  frequency  which 
give  an  equivalent  continuous  wave  breakdown.  Data  obtained  at  different  pressures 
indicate  that  the  percentage  deviation  of  pulse  breakdown  from  continuous  wave  break¬ 
down  decreases  as  the  pressure  increases;  thus  it  is  reasonable  to  assume  that  a  duty 
cycle  that  permits  equivalent  continuous  wave  breakdown  at  the  pressure  indicated  in 
the  figure  will  continue  to  do  so  as  the  pressure  increases. 

The  trend  of  these  constant  duty  cycle  curves  can  be  understood  through  the  reali¬ 
zation  that  the  time  required  to  obtain  breakdown  in  the  gas  is  the  time  required  for  the 
electron  density  in  the  gas  to  reach  the  ambipolar  transition  density  of  approximately 
10  electrons/cm  .  There  are,  however,  an  initial  number  of  electrons  left  over  from 
the  previous  breakdown  which  can  influence  the  time  required  to  obtain  the  necessary 
build-up  of  electron  density.  An  equivalent  continuous  wave  breakdown  is  obtained  when 
the  electron  density  builds  up  within  the  limits  of  a  given  pulse  width  to  the  ambipolar 
transition  density  from  an  original  density  determined  by  the  time  elapsed  since  the  pre¬ 
vious  breakdown.  The  build-up  of  electron  density  is  determined  by  the  free  diffusion  of 
electrons  in  hydrogen;  the  decay  of  electron  density  after  a  breakdown  is  determined  by 
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Fig.  II- 1 

The  variation  of  breakdown  field  as  a  function  of  pulse  width  and  pulse 
repetition  frequency  showing  the  effect  of  the  pulse  duty  cycle. 


the  ambipolar  diffusion  of  electrons  in  hydrogen  from  the  breakdown  density  to  the 
ambipolar  transition  density.  If  the  pulse  repetition  frequency  is  smaller  than  a  critical 
frequency,  practically  no  influence  is  felt  from  the  electrons  of  the  previous  discharge 
and  the  pulse  breakdown  for  a  given  pulse  width  should  no  longer  be  a  function  of  pulse 
frequency.  The  decay  of  electron  density  in  the  ambipolar  region  is  given  by 


where  nQ  =  initial  electron  density,  n  =  density  at  time  t,  D&  =  ambipolar  diffusion  coef¬ 
ficient,  A  =  diffusion  length  of  cavity,  and  t  =  time.  For  a  pressure  of  85  mm  of  Hg, 

and  values  of  (D  p)  =  700,  A  =  0.  4  cm,  n  =  10^,  and  n  =  10^,  a  diffusion  time  of  0.  13 

a  o 

sec  is  obtained  which  corresponds  approximately  to  a  repetition  frequency  of  7  to  8 
cycles  per  second. 

J.  J.  McCarthy 


B.  PLASMA  OSCILLATION 

To  further  the  investigations  of  D.  H.  Looney  (l),  a  tube  has  been  constructed  and 
put  into  operation  (2);  oscillations  and  a  standing- wave  pattern  have  been  observed.  The 
plasma-beam  interaction  region  was  three  times  the  maximum  length  used  by  Looney. 
The  purpose  of  the  new  tube  is  to  allow  independent  measurements  of  the  oscillation  fre¬ 
quency,  plasma  electron  density,  and  standing-wave  pattern  in  the  interaction  region. 
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The  variation  of  phase  velocity  with  electron  density  can  now  be  determined,  and  the 
assumption  that  the  plasma  electrons  behave  as  a  gas  which  propagates  an  acoustical 
pressure  wave  can  be  tested  using  various  models  of  the  electron  gas. 

E .  Gordon 
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C.  MICROWAVE  DETERMINATION  OF  THE  PROBABILITY  OF  COLLISION  OF 
SLOW  ELECTRONS  IN  NEON 

The  collision  probability  for  momentum  transfer  of  slow  electrons  in  Ne  will  be 
measured  with  the  microwave  method  of  Gould  and  Brown  (1).  Cataphoresis  will  be  used 
to  purify  neon  (spectroscopically  pure  Airco  neon)  from  the  other  rare  gas  impurities, 
and  activated  uranium  to  eliminate  the  impurities  that  are  not  rare  gases  (2). 

A.  L.  Gilardini 
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III.  SOLID  STATE  PHYSICS 


Prof.  W.  P.  Allis  Prof.  G.  G.  Harvey  J.  M.  Goldey 

Prof.  S.  C.  Brown  Prof.  L.  Tisza  J.  B.  Thomas 

I 

A.  DYNAMIC  REVERSIBILITY  IN  THE  KINETIC  THEORY 

Boltzmann's  integral  equation  is  generally  assumed  to  give  a  correct  account  of  the 
kinetic  processes  in  gases  when  the  intermolecular  interactions  can  be  described  in 
terms  of  binary  collisions.  The  purpose  of  the  present  investigation  is  to  draw  attention 
to  an  additional  limitation. 

The  fundamental  equations  of  physics  satisfy  the  principle  of  dynamic  reversibility: 
the  transformation  reversing  the  direction  of  time  t  —  -t  carries  over  every  solution  into 
another  dynamically  admissible  solution.  In  contrast,  the  phenomenological  description 
of  such  processes  as  electrical  conductivity  or  viscous  flow  is  irreversible. 

It  would  be  desirable  to  have  a  statistical  theory  that  contains  both  the  reversible  and 
irreversible  aspects  and  indicates  clearly  how  the  latter  arises  from  the  former. 

It  is  easily  shown  that  the  Boltzmann  equation  falls  short  of  this  requirement.  The 
theory  predicts  that  any  inhomogeneity  existing  at  an  initial  time  in  an  isolated  system 
decays  exponentially  until  the  quiescent  state  is  reached.  This  prediction  is  in  conflict 
with  the  principle  of  dynamic  reversibility  that  requires  spontaneous  departures  from  the 
quiescent  state. 

From  the  experimental  point  of  view  the  implications  of  the  Boltzmann  equation  are 
correct  for  signals  above  noise  level,  but  the  fluctuations  (thermal  noise)  provide  us  with 
a  clear  indication  that  spontaneous  deviations  from  equilibrium  have  been  too  drastically 
averaged  out  from  the  classical  theory. 

The  origin  of  this  shortcoming  has  been  submitted  to  a  detailed  analysis,  the  results 
of  which  may  be  helpful  for  the  formulation  of  a  more  general  kinetic  theory,  in  which 
resistivity  and  noise  are  accounted  for  within  the  same  formalism. 

L.  Tisza 


B.  MICROWAVE  STUDY  OF  SEMICONDUCTORS 


Measurements  of  the  carrier  contribution  to  the  dielectric  coefficient  of  germanium 
have  been  completed.  The  carrier  contribution  is  given  by  o-./coeo.  In  order  to  obtain 
an  effective  mass  from  these  measurements,  it  is  necessary  to  have  an  expression  for 
cr..  The  general  expression  for  cr  in  a  crystal,  from  which  cm  may  be  obtained,  is 


cr 


(V,  E)  2  9f°  , 

v  k  'a.  e _ o  ,3 

(vc  +  »  9E 


(1) 


where  f°  is  the  equilibrium  distribution  function  and  (V^E)q  is  the  gradient  of  the 
energy  along  the  direction  a  in  k  space.  In  order  to  obtain  an  expression  for  cr, 
therefore,  something  must  be  known  about  the  energy  surfaces.  We  have  obtained  an 
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expression  for  or  for  the  conduction  band  by  assuming  that  the  band  edge  occurs  along 
the  111  axis  and  that  the  surfaces  of  constant  energy  are  ellipsoids  of  revolution  about 
the  minimum.  This  is  the  model  recently  proposed  by  Lax  and  his  associates  (1). 

For  this  model  we  obtain 


2 

o-  =  - 

i  _  */  2  .  2' 


7  2  ,  2\ 

m  lvc  +  “  ) 


(2) 


where  m  is  given  by 


*  3m1m2 

m  = 


2m  ^  +  m^ 


(3) 


where  m^  is  the  longitudinal  mass,  and  m->  the  transverse  mass. 

For  the  same  energy  surfaces  Herring  (2)  has  recently  obtained  an  expression  which 
differs  slightly  from  Eq.  2.  His  result  may  be  approximated  for  the  present  case  to 


n^m** 


cr  .  = 
l 


1  - 


2  /  *x2  2 
gj  (m  )  p 


Be 


(4) 


where 


■d  _  3r(2+r) 

(1  +  2r)2 

where  r  =  m^/m^  and 


E 


<E> 


C  = 


E 


the  angular  brackets  indicating  Maxwellian  averages. 

The  valence  band  of  germanium,  according  to  recent  evidence  (3),  is  two-fold  degen¬ 
erate  with  the  band  edge  occurring  at  the  origin.  In  this  case,  the  holes  are  distributed 
between  the  two  degenerate  bands.  Here,  integration  of  Eq.  1  leads  to  Eq.  2  again,  but 
the  effective  mass  is  obtained  from  the  relation 


n 


m 


n. 


m 


n. 


+ 


m. 


where  n  =  n^  +  n^,  and  m^  and  m^  are  suitable  averages  of  the  masses  in  the  two 
bands. 
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The  measurements  were  made  at  two  temperatures.  For  the  N  material  the  con¬ 
ductivity  followed  the  law 

<r  =  o'  (100)1'66 

T  300  V  T  / 

This  value  of  the  exponent  is  in  agreement  with  that  obtained  by  other  workers  (4).  For 
the  P  type  we  obtained 


<T 


T 


=  o-  m] 

300  \  T  ) 


2.  17 


whereas  the  accepted  value  of  the  exponent  is  2.33  (4). 

The  measured  values  of  the  effective  mass  are  as  follows.  For  electrons  we  obtain 


m  =  0.  10  +  0.  05  m  where  we  have  used  Eq.  2 
-  o  ^ 


If  we  use  Eq.  4  we  have  m  =  0.  08 


+  0.  04  m  ,  where  we  have  used  a  mass  ratio  r  of  17  as  obtained  by  Lax  et  al.  (l)  in  cal¬ 
culating  B.  For  holes,  the  measured  value  of  the  mass  is  m*  =  0.33  +  0.  05  mQ  when 
we  have  used  a  value  of  2.  17  for  the  temperature  exponent.  If  the  temperature  exponent 
is  taken  as  2.33,  we  get  for  the  mass  m  =  0.28+0.05  mQ. 

These  results,  for  both  N  and  P  types  are  in  agreement  with  results  obtained  by 
cyclotron  resonance  experiments  (1,3)  and  magnetic  susceptibility  measurements  (5). 
They  are  not  in  agreement,  however,  with  the  results  of  an  experiment  similar  to  the 
present  one  (6,  7). 

J.  M.  Goldey 


References 

1.  B.  Lax,  H.  Zeiger,  R.  N.  Dexter,  and  E.  S.  Rosenblum,  Phys.  Rev.  93,  1418 
(1954). 

2.  C.  Herring,  personal  communication. 

3.  R.  N.  Dexter,  H.  Zeiger,  and  B.  Lax,  Phys.  Rev.  95,  557  (1954). 

4.  F.  J.  Morin,  Phys.  Rev.  93,  62  (1954). 

5.  J.  H.  Crawford  and  D.  K.  Stevens,  Phys.  Rev.  94,  1415  (1954). 

6.  T.  S.  Benedict  and  W.  Shockley,  Phys.  Rev.  89,  1152  (1953). 

7.  T.  S.  Benedict,  Phys.  Rev.  91,  1565  (1953). 


-9- 
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Prof.  M.  A.  Herlin 
Dr.  C.  W.  Garland 


R.  P.  Cavileer 
L.  D.  Jennings,  Jr. 


W.  M.  Whitney 
J.  W.  Wright 


H.  H.  Kolm 


A.  THE  VISCOSITY  OF  LIQUID  HELIUM 

Measurements  made  in  liquid  helium  indicate  that  thermomechanical  flow  caused  by 
heat  input  to  the  rotor  is  now  entirely  responsible  for  the  excessive  damping  observed; 
all  damping  from  the  magnetic  bearing  itself  has  been  eliminated.  The  results  of  two 
measurements  are  shown  in  Fig.  IV- 1,  which  represents  rotor  speed  as  a  function  of 
angular  displacement.  It  will  be  seen  that  the  timing  points  fall  on  straight  lines  to 
within  0.5  per  cent. 

The  first  curve,  taken  at  2.  17°K,  indicates  a  viscosity  of  41  micropoise,  which  is 
about  twice  the  value  observed  by  other  investigators.  The  measurement  was  continued 
until  the  level  of  liquid  helium  dropped  below  the  rotor;  the  slope  of  the  curve  is  seen 
to  change  discontinuously  at  this  point,  the  new  value  indicating  a  viscosity  of  5  micro¬ 
poise,  the  accepted  value  for  helium  gas  at  this  temperature.  This  observation  elimi¬ 
nates  the  possibility  of  instrumental  damping  in  the  most  direct  possible  way. 

Calculations  indicate  that  the  excessive  damping  in  liquid  helium  would  be  accounted 
for  by  a  heat  input  of  only  6  mw  to  the  rotor.  This  same  heat  input  would  produce 
an  apparent  viscosity  of  45.8  micropoise  at  the  lower  temperature  of  1.97°K  (the 
value  observed  by  other  investigators  is  12  micropoise),  which  is  seen  to  be  in  good 
agreement  with  the  observed  value  of  44.4  micropoise.  It  is  of  interest  to  note 
that  the  thermomechanical  damping  mechanism  in  question,  being  proportional  to  l/TS, 
increases  very  rapidly  below  1.5°K  and  might  account  for  the  observed  discrepancy 
between  all  oscillating  disk  measurements  on  the  one  hand,  and,  on  the  other,  some 
approximate  preliminary  results  obtained  by  a  static  method  (1). 

Figure  IV-2  shows  the  result  of  a  measurement  made  at  higher  speeds  in  an  effort 
to  determine  the  nature  of  the  nonlinear  effects  observed  by  other  investigators.  It  is 
seen  that  the  damping  line  does  not  show  curvature  but  changes  its  slope  nearly  discon¬ 
tinuously,  indicating  a  new  value  of  apparent  viscosity  at  higher  speeds.  This  observa¬ 
tion  is  in  agreement  with  the  "saturation  effect"  of  the  excess  decrement  observed  by 
other  authors  using  oscillating  disk  methods  (2). 

It  is  clear  from  measurements  made  thus  far  that  the  major  source  of  heat  input  is 
eddy  current  heating  owing  to  corrections  made  by  the  servomechanism  to  counteract 
mechanical  noise.  Major  changes  are  necessary  to  remedy  this  situation;  most  expe¬ 
diently,  removal  of  the  apparatus  to  a  more  favorable  location.  Further  work  will  not 
be  undertaken  until  this  difficulty  can  be  overcome. 


H.  H.  Kolm 
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Fig.  IV-1 


The  observed  decay  of 
gaseous  helium 


rotation  in  helium  II  and  in 
at  two  temperatures. 


Fig.  IV -2 

The  decay  of  rotation  in  helium  II  at  higher  speeds. 
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B.  HALL-EFFECT  AND  MAGNETO-RESISTANCE  EFFECT  MEASUREMENTS  IN 

RESISTANCE  MINIMUM  METALS:  DISCUSSION  OF  FEASIBILITY 

The  "resistance  minimum"  effect  in  magnesium  has  been  investigated  through 
measurement  of  electric  and  thermal  conductivities  and  thermoelectric  power.  Electric 
conductivity  has  been  investigated  through  a  mutual-inductance  technique  so  that  bulk 
samples  could  be  used.  In  order  to  separate  effects  of  the  mean  free  path  and  the  elec¬ 
tron  energy  level  density,  it  is  desirable  to  measure  the  Hall  effect.  Also,  the  magneto¬ 
resistance  effect  has  already  proved  enlightening  in  other  resistance  minimum  metals. 
Because  of  the  high  electrical  conductivity,  the  conventional  Hall-effect  measurement 
method  requires  fine  films  of  metal  to  provide  sufficiently  high  Hall  emfs,  and  the 
metallurgical  state  of  the  film  is  difficult  to  control.  It  is  desirable,  therefore,  to 
devise  a  mutual-inductance  method  for  this  measurement,  also.  Mutual  inductance 
techniques  for  measurements  on  magnetic  salts  have  been  improved  in  this  laboratory 
by  L.  D.  Jennings  and  the  writer.  It  is  now  possible  to  apply  these  methods  to  Hall- 
effect  measurements  with  adequate  accuracy. 

In  essence  the  mutual  inductance  as  a  function  of  applied  dc  magnetic  field  for  two 
different  coil  configurations  (described  below)  determines  both  the  magneto-resistance 
and  Hall  coefficients.  For  mathematical  simplicity  a  spherical  sample  is  chosen.  In 
the  first  configuration,  the  axes  of  the  primary  and  secondary  coils  are  aligned  parallel 
to  the  applied  dc  magnetic  field.  The  Hall  effect  is  then  radial  and  contributes  nothing 
to  the  emf  induced  in  the  secondary  coil.  The  effect  of  unavoidable  misalignment  of  the 
secondary  coil  has  been  computed  and  is  unimportant.  This  is  essentially,  then,  the 
experiment  of  Rorschach  and  Herlin;  the  conductivity  is  obtained  as  a  function  of 
magnetic  field. 

In  the  second  configuration  the  primary  and  secondary  coils  and  the  magnetic  field 
are  mutually  perpendicular.  The  circulating  currents  induced  by  the  primary  field  are 
diverted  by  the  dc  field,  producing  an  induced  emf  in  the  secondary.  In  the  conventional 
Hall  experiment,  the  current  lines  are  constrained  to  flow  parallel  to  edges  of  the  metal 
by  the  Hall  field  caused  by  edge  charges.  In  this  method  the  current  is  actually  diverted, 
and  the  change  in  mutual  inductance  gives  the  product  of  the  conductivity  and  the  Hall 
constant.  The  conductivity  measurement  is  therefore  a  necessary  companion  experi¬ 
ment  . 
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In  the  second  configuration,  the  alignment  of  the  primary  with  respect  to  the 
secondary  is  critical,  but  the  alignment  of  either  coil  with  respect  to  the  magnetic  field 
is  not.  The  effect  of  primary-secondary  misalignment  may  be  corrected  by  rotation  of 
the  dc  field  ninety  degrees  into  the  direction  of  the  primary.  The  mutual  inductance  thus 
measured  will  be  that  without  the  Hall  effect  and  will  constitute  the  major  correction  due 
to  misalignment. 

With  realizable  coils,  the  anticipated  change  of  mutual  inductance  over  the  whole 
range  of  magnetic  field  up  to  10,  000  gauss  is  of  the  order  of  0.  1  to  1  microhenry.  The 
bridge  to  be  used  is  sensitive  to  0.001  microhenry,  so  that  a  precision  of  one  per  cent 
should  prevail. 

J.  W.  Wright 
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A.  A  PARAMAGNETIC  RESONANCE  STUDY  OF  AMMONIUM  CHROME  ALUM 


A  detailed  solution  of  the  fine  structure  splitting  of  the  trivalent  chromium  ground 
state  has  been  calculated.  The  effect  of  the  neighbors  of  the  paramagnetic  ion  was 
assumed  representable  by  an  electric  potential  of  cubic  symmetry  with  a  slight  trigonal 
(axial)  distortion  ( 1).  Matrix  elements  of  the  orbital  splitting  were  obtained  by  the 
method  of  Racah  (2)  assuming  only  that  the  expression  for  the  equivalent  crystalline 
potential  obeys  the  same  commutation  rules  with  the  angular  momentum  that  spherical 
harmonics  do.  The  effect  on  energy  levels  and  intensities  of  spin-orbit  coupling  and  of 
a  magnetic  field  at  an  arbitrary  angle  to  the  axial  distortion  was  calculated  to  second 
order.  Otherwise,  orientation  of  the  magnetic  field  in  the  cubic  coordinates  does  not 
affect  the  energy  levels,  but  does  alter  line  intensities  for  specimen  orientations  of  two¬ 
fold  (and  less)  symmetry. 

The  degeneracy  removal  by  orientation  dependent  off-diagonal  terms  and  the  diffi¬ 
culty  in  assigning  uniformly  suitable  designations  to  the  energy  levels  result  in  three 
discontinuities  in  matrix  elements  and  in  field  derivatives  of  energy  levels  as  a  function 
of  orientation  and  field. 

Line  positions  and  intensities  were  calculated  and  checked  with  experiment.  Posi¬ 
tions  corresponded  within  0.  2  per  cent  except  where  measurement  was  unusually  diffi¬ 
cult  because  of  very  weak  absorptions  or  interference  between  two  overlapping  lines. 
Intensities  generally  correlated  within  a  factor  of  two  (within  much  less  when  the 
lines  were  well  separated)  and  had  "regular"  shape.  In  each  case,  correction  was  made 
for  linewidth  and  for  observation  of  absorption  at  constant  frequency  instead  of  constant 
magnetic  field.  Linewidths  and  shapes  can  be  accounted  for  by  assuming  a  distribution 
of  crystalline  electric  fields  and  correcting  to  a  constant  field  observation.  Certain  lines 
have  largely  the  character  of  AS  =  2,  3.  These  have  been  known  as  "forbidden"  (3)  or 
"double-jump"  (4)  lines.  The  check  included  these  absorptions,  whose  intensities  were 

calculated  on  the  basis  of  the  extent  of  their  AS  =  1  character. 

z 

The  magneto-gyric  ratio  was  studied  as  a  function  of  crystal  dilution  and  tempera¬ 
ture  (5).  At  great  dilution  ( 1:50)  it  was  1 . 974  ±  0.  002;  at  full  strength,  1.988  ±0.002 
independent  of  temperature  within  the  experimental  accuracy.  This  implies  a  smaller 
"cubic  field"  in  the  aluminum  lattice  than  in  the  chromium  lattice. 

The  quartet  ground  state  (three  unpaired  electrons)  can  be  split  only  as  far  as 
two  doublets.  A  theorem  due  to  Kramers  (6)  shows  that  these  doublets  cannot  be 
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decomposed  by  an  electric  field.  The  separation  of  these  energy  levels  (zero  field 
splitting,  6)  can  easily  be  measured  experimentally  and  amounts  to  0.  098  cm  1  for  the 
very  dilute  material,  0.  135  cm  ^  for  full  strength,  at  25°C.  These  values  decrease  at 
about  the  rate  of  0.  0005  cm  ^ /°K  to  80°K  (7).  This  implies  a  strong  decrease  in  axial 
distortion  with  lower  temperature. 

At  about  80°  K,  the  full  strength  alum  and  alums  of  low  dilution  undergo  a  second - 
order  phase  transition.  Characteristics  of  this  transition  have  been  studied  elsewhere 
in  optical,  thermodynamic,  and  other  paramagnetic  resonance  work  (8,  9,  10).  Powder 
x-ray  diffraction  patterns  were  compared  above  and  below  the  transition.  Line  spacing 
was  identical  within  the  resolution  of  the  lines,  but  intensities  changed  somewhat.  The 
heavier  atoms  undoubtedly  do  not  change  position  appreciably.  Below  the  transition  two 
microwave  splittings  have  been  observed  which  can  only  be  accounted  for  by  the  exist¬ 
ence  of  the  chromium  ions  in  two  possible  environments. 

The  magnet,  microwave,  control,  and  presentation  apparatus  used  in  this  research 
will  be  described  in  two  forthcoming  Review  of  Scientific  Instruments  articles.  This 
problem  has  been  the  subject  of  a  thesis  submitted  to  the  M.  I.  T.  Department  of  Physics. 

C.  F.  Davis,  Jr. 
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B.  HINDERED  ROTATION  IN  SEMISYMMETRIC  MOLECULES 

A  molecule  which  may  be  described  as  a  semisymmetric  hindered  rotor  is  one 
that  consists  of  a  symmetric  top  attached  at  an  arbitrary  angle  to  a  general  asymmetric 
framework,  and  allowed  to  rotate  about  its  symmetry  axis  relative  to  the  frame. 


-15- 


(V.  MICROWAVE  SPECTROSCOPY) 


A  well-known  example  of  this  type  of  molecule  is  methyl  alcohol,  CH^-OH.  If  the  two 
groups  are  taken  as  rigid,  then  the  microwave  spectrum  will  be  due  to  the  interaction 
of  hindered  rotation  with  over-all  rotation. 

A  theory  is  being  developed  for  this  general  type  of  molecule  which  will  give  the 
Hamiltonian  in  terms  of  principal  moments  of  the  molecule  and  the  angle  of  the  top 
relative  to  principal  axes.  The  pure  hindered-rotation  spectrum  will  then  be  calculated 
and  compared  to  the  experimental  results  for  sample  molecules. 

P.  R.  Swan,  Jr. 
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A.  NUCLEAR  ORIENTATION  IN  MERCURY  VAPOR 

A  new  attempt  is  being  made  to  obtain  nuclear  orientation  in  mercury  vapor  by 
"optical  pumping"  with  circularly  polarized  light  (1,  2,  3).  Although  this  has  been  done 
in  a  beam  (4,  5),  the  successful  use  of  a  vapor  would  facilitate  the  orientation  of  small 
amounts  of  radioactive  isotopes. 

To  prevent  condensation  of  atoms  containing  oriented  nuclei,  a  hot  quartz  cell  was 

sealed  off  with  unsaturated  vapor  at  the  vapor  pressure  corresponding  to  ice  tempera- 

198 

ture.  The  cell  will  be  irradiated  with  light  from  an  Hg  arc  in  a  very  weak  field.  In 

such  a  field,  the  transition  to  the  F  =  3/2  level  of  the  ^P.  state  of  Hg^^  overlaps  with 

198  ^  201 
the  single  line  of  Hg  which  can  be  used,  therefore,  to  orient  Hg  .  Since  a  large 

magnetic  field  is  not  needed,  the  arc  can  be  brought  very  close  to  the  cell,  giving  con¬ 
siderably  more  light  intensity  than  was  available  in  previous  attempts.  The  use  of  a  198 
arc,  rather  than  natural  mercury,  will  considerably  reduce  the  background  from  the  even 
isotopes.  It  is  planned  to  detect  orientation  by  the  method  of  Hawkins  and  Dicke  (5),  that 
is,  by  observing  the  degree  of  polarization  of  the  scattered  light  as  a  function  of  a  mag¬ 
netic  field  parallel  to  the  axis  of  the  light  beam.  The  degree  of  polarization  is  a  function 
of  the  degree  of  orientation.  The  component  of  the  earth's  field  transverse  to  the  light 
beam  axis  will  cause  a  disorienting  precession  of  the  nuclear  magnetic  moment  which 
should  have  a  maximum  effect  when  the  applied  field  just  cancels  the  component  of  the 
earth's  field  parallel  to  the  light  beam  axis.  Calculations  indicate  that  with  the  light 

intensity  available,  it  should  be  possible  to  detect  orientation  if  the  thermal  relaxation 

201 

time  in  the  ground  state  of  Hg  is  greater  than  approximately  0.  2  sec. 


P.  L.  Sagalyn,  J.  R.  Stanley 
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B.  NUCLEAR  MAGNETIC  RESONANCE  IN  GASES 


In  the  nuclear  magnetic  resonance  of  B"  in  BF^  gas  at  various  pressures  with 
various  added  gases,  simultaneous  measurements  of  linewidth  AH  and  signal  amplitude 
vs.  applied  rf  amplitude  were  made.  Theory  gives 


signal^ 
>  ) 


rf— 0 


3 

where  N-d  is  the  number  of  atoms  of  B"  per  cm  and  AH  is  the  linewidth  as  rf  —  0.  In 
r> 

these  experiments  AH  was  principally  due  to  field  inhomogeneity.  At  an  rf  value  such 
that 


( signal \  _  1  /signal 

V  rf  P  rf  . 


rf-0 


theory  gives 

<*h1>2tiW=1 

In  pure  BF^,  T^  is  predicted  to  be  proportional  to  the  collision  frequency.  From 
kinetic  gas  theory,  in  pure  BF^,  collision  frequency  is  proportional  to  Ng.  In  a  mixture 
one  expects  a  similar  dependence  on  all  the  types  of  collisions,  although  there  may  be 
some  gas  components  of  greater  effectiveness  than  others  . 

One  can  thus  define  an  equivalent  Bn  density  Ng  such  that  pure  BF^  of  that  density 
would  produce  the  observed  Tj. 

The  following  quantities  were  obtained  from  measurement: 

W  a  AH  E  ac  H, 
o  cl 


Pg  =  BF^  partial  pressure,  P  =  total  gas  pressure 
We  see  from  the  above  that  in  pure  BF^>  T^  is  proportional  to  Ng  but  also,  in  any  case. 


AU  w 

T  a  ac  — — 

1  2  2 


H 


1 


E' 

c 


Hence  we  expect  that 
W 


E*  2  Nd 
c  B 


AoWo 

-  =  constant.  Also,  we  see  that  — ^ -  =  constant. 


N 


B 


Taking  the  ratio,  we  have  AqE^  =  C^.  Now,  assuming  Boyle's  law,  Ng  is  proportional 
to  Pg,  we  get 
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Wq  AW 

~2 -  =  constant  =  and  — p -  =  constant  =  C^. 


E'  P  -r-, 

c  B 


B 


Thus  C3  =  C2/Cr 

For  mixtures  we  still  have  (AQWo)/Pg  =  C2-  And  by  definition, 
W 

- - —  =  C 

2  1 
E'  P„  1 
c  E 


Now 


C  ,3  =  A  E'2^ 
3  C,  o  c  P 


T 


B 


We  can  evaluate  from  pure  BF^  runs,  that  is. 


C,  =  (a 

3  Vo 


E' 

c 


averaged  over  the  pure  BF^  runs.  Thus 


P  2 

B  A  E' 
o  c 


We  can  now  compare  the  effectiveness  of  an  added  gas  to  that  of  BF^: 


We  find  for 

H2  :  p  ~  0.  4  :  one  run,  Pg  ~  1400  psi,  Pg  ~  300  psi 

2 

N2  :  p  ~  0.  8  :  two  runs,  P^  ~  700  psi,  Pg  ~  1000  psi 

P^  ~  300  psi,  Pg  ~  400  psi 
2 

On  a  purely  collision  frequency  basis,  one  expects,  because  of  the  mass  differences  for 
H2,  p  ~  4;  for  N2,  p  ~  1.3.  Thus  some  other  effect  must  be  present,  perhaps  the  dif¬ 
ficulty  of  transfer  of  angular  momentum  from  a  light  to  a  heavy  rotator. 

As  a  check,  the  constant  C2  was  evaluated  for  each  of  eleven  runs.  In  all  but  one, 
C2  fell  in  the  range  1.4  -  2.  2.  This  seems  to  be  the  experimental  error.  In  spite  of 
an  even  larger  spread  in  obtained  from  data  on  pure  BF^,  a  pressure  range  of  a 
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factor  3  clearly  shows  pressure  narrowing  as  expected.  The  value  of  T^  obtained  is 
2.  3  msec/psi  for  pure  BF^.  The  pressures  used  were  around  1000  psi  (4  runs)  and 
around  300  psi  (2  runs).  The  linewidth  was  approximately  7  0  milligauss  for  all  runs, 
which  corresponds  to  a  T^  of  3  l/2  msec.  One  run  was  made  with  about  7  0  psi  of  O^ 
and  1000  psi  of  BF^.  Oxygen,  because  of  its  strong  paramagnetism,  causes  a  decrease 
of  Tj  with  increase  in  oxygen  pressure.  The  predicted  change  was  negligible  against 
the  experimental  error.  No  change  was  found.  This  seems  to  preclude  small  traces 
of  air  being  able  to  produce  variations  in  the  results. 

N.  I.  Adams  III 


C.  STEPWISE  EXCITATION  OF  MERCURY 

The  structure  of  the  mercury  spectral  line  X4047A  is  being  studied  with  methods 

described  in  previous  reports  from  this  group.  This  line  arises  from  transitions  from 

3  3 

the  state  S^  (valence  electrons  6s  and  7s)  down  to  the  metastable  state  PQ  (6s,  6p). 

(See  Fig.  VI- 1.)  Our  procedure  is  to  measure  the  absorption  of  the  mercury  line  X4047A 

3 

in  mercury  vapor  containing  atoms  in  the  Pq  state  as  a  function  of  the  frequency  of 
this  incident  radiation. 

To  provide  a  density  of  metastable  atoms  sufficient  for  this  experiment,  a  quartz 

resonance  cell,  containing  a  drop  of  mercury,  is  either  sealed  off  with  a  known  pressure 

of  nitrogen,  or  connected  to  a  vacuum  system  so  that  the  atmosphere  may  be  controlled. 

1  3 

Atoms  are  excited  from  the  ground  state  Sn  (6s,  6s)  up  to  the  state  P.  (6s,  6p)  by 

u  3*3 

illuminating  with  the  resonance  line  X2537A.  Since  the  energy  level  P^  to  PQ  is  in 
resonance  with  some  vibrational  levels  of  the  nitrogen  molecule,  collisions  while  the 
mercury  atom  is  excited  knock  some  mercury  atoms  into  the  metastable  state.  It  has 
been  found  that  the  intensities  and  lifetimes  involved  are  such  that  a  population  of  meta¬ 
stable  atoms  which  is  sufficiently  large  and  constant  may  be  maintained. 

A  source  of  X4047A  light  of  variable  frequency  is  provided  by  a  microwave  discharge 

198 

lamp  containing  only  the  isotope  Hg  .  The  light  is  taken  parallel  to  a  variable  mag¬ 
netic  field,  and  one  of  the  Zeeman  components  is  removed.  Absorption  of  the  X4047A 
beam  in  the  mercury  vapor  is  detected  and  measured  by  a  photomultiplier  that  has  been 
filtered  to  be  sensitive  only  to  the  line  X5461A.  This  is  one  of  the  possible  wavelengths 

3 

emitted  as  the  atom  subsequently  falls  from  the  level. 

The  absorption  vs.  frequency  patterns  thus  observed  (Fig.  VI-2)  show  the  same 

qualitative  features  as  have  been  found  by  high- re  solution  spectroscopy.  (For  a  spectral 

plate  of  this  line,  see  reference  1.)  The  even  isotopes  fall  in  a  broad  central  peak. 

Since  they  do  not  have  nuclear  spin,  they  have  no  corresponding  magnetic  hyperfine 

structure,  and  the  small  shift  owing  to  isotopic  volumes  has  not  yet  been  resolved.  A 

199 

tentative  identification  as  hfs  components  corresponding  to  the  odd  isotopes  Hg  and 
201 

Hg  has  been  made  for  all  of  the  subsidiary  absorption  peaks  observed.  The  positions 
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Fig.  VI- 1 

Schematic  term  diagram  of  those  mercury  levels  giving  rise 
to  the  hyperfine  structure  of  the  lines  \2537A  and  \4047A. 


Fig.  VI-2 

Preliminary  determination  of  the  structure  of  the  line  \4047A  obtained 
by  magnetic  scanning.  Vertical  markers  indicate  positions  of  lines  to 
be  expected  from  previous  spectroscopic  data.  This  curve  has  not  been 
corrected  for  variation  of  background  signal  with  the  leakage  field  from 
the  magnet,  which  may  be  seen  in  the  curving  base  line. 
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of  these  peaks  have  been  found  to  satisfy  the  interval  rule  and  their  centers  of  gravity 
fall  approximately  where  expected. 

An  attempt  was  made  to  determine  which  isotopes  should  be  excited  into  the  meta¬ 
stable  state,  with  a  pure  isotope  lamp  as  the  source  of  the  first  resonance  excitation. 

However,  even  if  the  frequency  of  the  \2537A  radiation  brings  only  atoms  of  isotopes 

3 

198  and  201  into  the  P  state,  absorption  peaks  are  found  in  the  metastable  atoms  that 

1  199 

correspond  to  the  hfs  of  Hg  .  When  the  mercury  pressure  is  that  of  the  saturated 
vapor  at  room  temperature,  collisions  and  imprisonment  of  resonance  radiation  act  to 
redistribute  the  excitation  energy  to  all  the  isotopes. 

H.  H.  Plotkin 
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A.  EQUIVALENT  CIRCUIT  FOR  A  THREE -TERMINAL  PAIR  NETWORK 


The  study  of  microwave  amplifiers  of  the  beam  type  requires  at  times  (1)  an  equiva¬ 
lent  circuit  representation  of  nonreciprocal  three-terminal  pair  networks.  The  equiva¬ 
lent  circuit  is  useful  if  it  leads  to  mathematical  simplifications. 

A  simple  equivalent  circuit  for  the  most  general  lossless  nonreciprocal  three- 
terminal  pair  network  can  be  constructed  with  the  aid  of  a  lossless  three-terminal  pair 
gyrator.  The  gyrator  is  defined  by  the  impedance  matrix  Z. 


Z  = 


0 

-x 


x  y 

0  z 


-y 


-z 


0 


The  transfer  impedances  x,  y,  and  z  are  real  quantities. 

The  most  general,  lossless,  nonreciprocal,  three-terminal  pair  network  has  an 
impedance  matrix  of  the  form  (2) 


Z  =  j 
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The  X's  are  all  real  quantities.  The  matrix  has  thus  9  adjustable  parameters.  The 
equivalent  circuit  shown  in  Fig.  VH-1  also  has  9  adjustable  parameters.  It  includes, 
besides  the  gyrator  and  lumped  reactances,  ideal  two-terminal  pair  phase  shifters  (2) 
characterized  by  the  matrix  of  generalized  circuit  parameters  (3) 


Fig.  VH-1 

The  equivalent  circuit. 
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e^*  0 

0  e^* 


(1) 


The  arrows  above  the  ideal  phase  shifters  in  Fig.  VII- 1  point  from  the  input  to  the  output 
terminals  of  the  phase  shifters.  The  voltage  and  current  at  the  input  of  a  phase  shifter 
with  the  matrix  (Eq.  1)  lead  the  voltage  and  current  at  the  output  by  a  phase  angle  4>  . 

Below,  the  elements  of  the  equivalent  circuit  are  given  in  terms  of  the  9  parameters 
of  the  three-terminal  pair  network: 


X1  =  Xn  x  =  X1?  y  =  X 


12 


T  3 


*1  =  *13  +  2  U  *2  *13  '  *12 


X23 
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3  33  X 


12 


cos  (<|>23  +  4>12  -  4>13) 
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B.  THE  MINIMUM  NOISE  FIGURE  OF  LONGITUDINAL  BEAM  AMPLIFIERS 

Some  general  relations  were  derived  concerning  the  minimum  noise  figure  of  micro- 
wave  amplifiers  of  the  beam  type.  A  general  model  of  a  microwave  amplifier  with  arbi¬ 
trary  noise  smoothing  schemes  was  taken  as  a  starting  point.  The  model  assumes  that 
the  amplifier  can  be  split  into  an  input  and  an  output  section.  The  rf  coupling  between 
the  input  and  output  section  is  provided  solely  through  the  electron  beam.  A  klystron, 
a  velocity  jump  amplifier,  and  a  traveling-wave  tube  with  a  severed  helix  are  all  special 
cases  of  this  model.  The  interaction  between  the  beam  and  the  rf  circuit  of  such 
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Fig.  VII- Z 

The  equivalent  circuit  of  a  microwave  amplifier. 


amplifiers  can  be  represented  by  an  equivalent  circuit  shown  in  Fig.  VII-2;  that  is,  by 
two  linear  nonreciprocal  passive  three-terminal  pair  networks  connected  in  cascade. 
(The  reader  is  referred  to  reference  1  for  the  definitions  and  concepts  used  in  this 
section.)  The  terminals  marked  (1),  (2),  (1'),  and  (2')  are  the  "beam  terminals."  The 
voltage  and  current  at  those  terminals  are  the  ac  beam  voltage  and  beam  current.  The 
terminals  marked  (3)  and  (3')  are  the  input  and  output  terminals  of  the  amplifier. 

A  nonreciprocal  passive  three-terminal  pair  network  is  characterized  by  9  complex 
quantities.  However,  for  our  purpose  it  is  sufficient  to  characterize  the  output  network 
by  two  complex  quantities,  a  and  b.  They  are  defined  as  follows:  When  an  ac  beam 
voltage  V L ( 1 ')  and  an  alternating  beam  current  1^1')  are  applied  to  the  terminal  pair  1' 
of  the  output  circuit,  the  voltage  across  the  output  load  is  because  of  linearity 

V  ( 3 ' )  =  aV  x  ( 1 ' )  +  bl  L  ( 1 1 )  (1) 

The  subscript  "1"  of  the  voltage  and  current  indicates  small  signal  beam  voltages  and 
currents. 

One  of  the  relations  derived  concerns  the  general  form  of  the  minimum  noise  figure 
of  the  amplifier  described  above.  The  following  definition  of  the  noise  figure  is  used: 

Noise  power  output  per  unit  frequency  from  beam  noise 

F  =  1  +  - — - - - ; - - —  (2) 

Noise  power  output  per  unit  frequency  from  input  circuit  noise 

The  above  definition  is  applicable  if  the  input  circuit  noise  and  the  beam  noise  are 
uncorrelated.  The  denominator  of  Eq.  2  is  proportional  to  kT,  where  T  is  the  tempera¬ 
ture  of  the  input  circuit. 

The  numerator  of  Eq.  2  can  be  evaluated  by  techniques  of  four-terminal  network 
transformations  of  beam  noise.  The  Thevenin  impedance  of  the  input  circuit  is  connected 
across  the  terminal  pair  3  of  Fig.  VII-2  with  the  result  that  the  three-terminal  pair 
input  network  can  be  treated  as  a  two-terminal  pair  network  so  far  as  its  action  upon  the 
beam  noise  is  concerned.  Under  the  assumption  that  a  lossless  four-terminal  beam 
transducer  between  the  cathode  and  the  amplifier  is  optimized  for  lowest  noise  figure, 
we  find  for  the  noise  figure  the  general  expression 
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SQ  and  II  are  two  parameters  of  the  beam  noise  at  a  plane  beyond  the  potential  minimum 
in  the  electron  gun.  K  is  a  constant  dependent  upon  the  characteristic  parameters  of 
the  amplifier.  The  constant  K  can  be  greater  or  less  in  magnitude  than  unity.  It  is  of 
great  interest  to  find  the  dependence  of  K  upon  the  particular  amplifying  structure  used. 
This  has  been  done  under  the  assumption  that  the  input  three-terminal  pair  network  of 
Fig.  VII-2  is  lossless.  The  equivalent  circuit  for  the  three -terminal  pair  network  of 
the  preceding  section  has  been  used.  It  is  found  that  two  cases  have  to  be  distinguished 

Case  A:  Re(ab*)  <  0  (4) 


Inequality  4  can  be  interpreted  as  follows:  Eq.  1  shows  that  a  modulation  on  the  beam 
at  the  point  (l1)  with  a  small  signal  impedance 


Vjfl') 
Ix(l ')  = 


b 

a 


(5) 


does  not  deliver  any  power  to  the  output  circuit.  According  to  inequality  4,  the  critical 
impedance  5  must  have  a  positive  real  part  for  case  A  to  apply.  A  long  helix  termi¬ 
nated  in  a  matched  load  is  a  known  example  of  such  an  output  structure. 

If  case  A  applies,  it  can  be  shown  that  the  minimum  possible  value  of  K,  K  .  ,  is 


K  .  =  1 

min 


(6) 


Case  B:  Re(ab  )  >  0 


The  critical  impedance  5  for  case  B  has  a  negative  real  part.  A  cavity  with  a  long  gap 
(a  transit  angle  0  in  the  range  0  <  0  <  2tt)  is  a  known  example  of  such  a  structure. 

If  case  B  applies 

K  <  1  (7) 

and  indeed  can  be  zero  or  become  negative. 

It  is  obviously  of  great  interest  to  determine  whether  any  useful  amplifier  structures 
exist  with  a  K  appreciably  less  than  unity.  This  question  has  been  investigated  under 
the  additional  assumption  that  the  output  three-terminal  pair  network  of  Fig.  VU-2  is 
lossless.  Under  this  assumption  it  was  possible  to  show  that  the  available  power  gain 
G  of  the  amplifier  is  limited  by  the  magnitude  of  K  by 
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Thus,  a  structure  with  a  K  appreciably  less  than  unity  and  a  correspondingly  low  noise 
figure  would  be  useless  as  an  amplifier  because  of  its  low  gain.  It  should  be  emphasized 
that  relations  6,  7,  and  8  have  been  obtained  under  the  restrictive  assumptions  that 
either  the  input  three -terminal  pair  network  or  the  input  and  output  networks  are  loss¬ 
less.  A  generalization  of  the  proof  to  include  loss  encounters  great  mathematical  dif¬ 
ficulties. 

H.  A.  Haus 
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C.  BEAM  NOISE  MEASUREMENTS 


Some  measurements  were  made  on  a  gun  of  higher  perveance  than  was  heretofore 
used.  This  gun  had  a  convergence  half  angle  of  13°,  a  perveance  of  0.35  x  10  ,  and 

produced  a  beam  about  0.  08  inch  in  diameter.  The  heater  construction  was  such  as  to 
produce  an  observable  magnetic  field  at  the  cathode  surface  (a  few  gauss).  As  with  the 
previous  gun,  a  strong  growing  noise  wave  was  observed;  however,  when  viewed  on  an 
oscilloscope,  it  fluctuated  violently  at  60  cps.  This  was  soon  traced  to  the  heater.  It 
was  observed  that  the  growing  noise  wave  might  vary  by  about  10  db,  depending  upon 
the  polarity  of  the  heater  current,  and  that  this  bore  a  definite  relation  to  the  direction 
of  the  main  focusing  field;  reversing  one  required  reversing  the  other  to  produce  a 
given  effect.  Previous  measurements  had  indicated  that  the  main  focusing  field  of 
several  hundred  gauss  produced  a  leakage  field  of  a  few  gauss  at  the  cathode,  outside 
the  magnetic  field. 

In  order  to  pin  this  down  more  specifically,  an  auxiliary  winding  was  wrapped  around 
the  steel  tube  surrounding  the  gun.  With  this  it  was  possible  to  produce  fields  of  a  few 
gauss  at  the  cathode,  and  with  the  proper  polarity  the  growing  noise  wave  could  be  elimi¬ 
nated.  Subsequent  tests,  when  the  gun  was  removed  from  the  system,  indicated  that  the 
noise  was  reduced  when  the  two  fields  at  the  cathode  were  aiding,  and  the  noise  was 


enhanced  when  they  were  opposing. 

A  more  careful  experiment  is  being  prepared  in  order  to  get  a  better  picture  of  this 


phenomenon. 


C.  Fried,  L.  D.  Smullin 


D.  THE  INTERNALLY  COATED  CATHODE 

Previous  investigation  of  internally  coated  cathodes  made  from  thin  (less  than  0.  005 
inch)  nickel  with  holes  approximately  0.015  inch  in  diameter  has  shown  that  the  electrons 
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Simple  electron  sources  with  internally  cooled  cathodes. 
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Tube  allowing  determination  of  regions  of  emission. 
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originate  from  the  coating  and  are  drawn 
through  the  hole  by  the  applied  anode 
voltage.  The  more  field  that  can  be  made 
to  penetrate  the  hole,  the  greater  the 
emission.  If  such  cathodes,  then,  are  to 
be  used  as  electron  sources,  we  must  have 
"extraction"  anodes  similar  to  the  simple 
guns  shown  in  Fig.  VII- 3(a)  and  (b). 

The  source  shown  in  Fig.  VII- 3(a)  has 
produced  current  densities  of  approxi- 
mately  1  amp/cm  at  moderate  tempera¬ 
tures  and  less  than  100  volts  with  about 
10  per  cent  interception  current.  The 
source  in  Fig.  VII- 3(b)  shows  no  inter¬ 
ception  current,  and  the  construction  and  operation  are  not  so  critical.  The  presence  of 
the  extraction  anode  necessarily  produces  a  lens  combination  that  is  divergent,  and  a 
strong  field  on  the  outside  of  this  anode  is  necessary  to  remove  this  divergence.  A 
magnetic  field  might  then  hold  together  a  hollow  beam,  or  a  strongly  convergent  lens 
might  produce  a  crossover  point  to  be  imaged  by  conventional  means  in  a  cathode- ray 
tube. 

Investigation  of  the  cathode  shown  in  Fig.  VII-4  with  a  small  movable  anode  hole 
backed  by  a  collector  gives  us  an  insight  into  the  origin  of  the  electrons  by  letting  us  see 
the  effect  of  the  probe  voltage  on  the  amount  and  pattern  of  the  emission.  The  significant 
results  are  illustrated  in  Fig.  VII- 5.  Here,  the  sampled  current  through  the  anode  hole 
to  the  collector  is  plotted  as  a  function  of  the  distance  the  anode  hole  has  moved  across 
a  diameter  of  the  cathode  hole.  The  results  of  this  experiment  seem  to  be  the  same  as 
those  of  a  previous  similar  experiment;  the  electrons  originate  behind  the  hole,  but  the 
lens  action  of  the  cathode  hole  is  uncertain. 

During  this  investigation,  one  cathode  showed  a  surprising  increase  in  emission 
(about  ten  times)  after  the  coating  had  been  bombarded  by  electrons  of  approximately 
50  volts  from  an  emitter  inside  the  cavity.  This  state  of  activity  occurred  instantly  and 
has  persisted.  Possible  reasons  for  this  phenomenon  were  the  production  of  free  barium 
by  the  bombardment  or  excessive  temperature  of  the  coating  (without  excessive  tempera¬ 
ture  of  nickel). 


h 


HOLE  DIAMETER  - 


Fig.  VII- 5 

Current  density  across  beam  a  few  hole 
diameters  in  front  of  the  cathode. 
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A.  CESIUM  CLOCK 

During  the  last  quarter  the  cesium  beam  apparatus  which  has  been  under  construction 
to  provide  a  primary  frequency  standard  for  the  Atomic  Beam  laboratory  has  been  in 
successful  operation.  Much  useful  information  has  been  obtained  showing  that  such  a 

system  can  be  expected  to  provide  a  long  term  stability  approaching  1  part  in  10*°. 

133 

The  system  uses  the  field  insensitive  resonance  line  of  Cs  ,  occurring  at  approxi¬ 
mately  9192. 63197  Mc/sec,  as  the  frequency  standard.  By  exploiting  the  molecular 
beam  techniques  developed  for  measurements  on  scarce  isotopes,  which  involve  the  use 
of  a  narrow  beam  and  sensitive  detector,  an  apparatus  has  been  designed  to  give  a  good 
signal-to-noise  ratio  using  only  one  microgram  of  cesium  per  day.  The  expected  line- 
width  was  about  2  00  cps,  making  possible  a  stability  of  1  part  in  lO1^  by  splitting  the 
linewidth  to  +1  per  cent. 

For  those  not  familiar  with  the  principles  involved,  the  following  brief  description 
is  provided. 

There  are  a  number  of  magnetic  resonance  lines  associated  with  a  cesium  (or  other) 
atom,  resulting  from  transitions  between  two  quantum  states.  There  exists  a  pair  of 
states  for  which  the  energy  is  independent  (to  the  first  order)  of  the  external  magnetic 
field,  so  that  the  associated  frequency  is  similarly  field  insensitive.  (There  is,  how¬ 
ever,  a  small  term  proportional  to  the  square  of  the  field.)  The  resonance  is  excited 
by  applying  a  magnetic  field  at  the  appropriate  frequency,  the  width  of  the  resonance 
curve  being  inversely  proportional  to  the  time  the  atom  spends  in  the  field.  With  a  beam 
of  atoms  a  long  path  without  collision  can  be  obtained  with  a  high  density  of  atoms, 
because  of  the  ordered  nature  of  their  paths.  The  transition  in  question  is  associated 
with  a  reversal  in  the  magnetic  moment  of  the  atom.  If  the  beam  of  atoms  is  passed 
through  a  transverse  inhomogeneous  magnetic  field,  a  deflection  is  produced  in  a  direc¬ 
tion  determined  by  the  sign  of  the  magnetic  moment.  If  a  pair  of  such  fields  is  placed 
on  either  side  of  the  rf  field  (A  and  B  magnets),  a  distinction  can  be  made  between 
atoms  which  have  made  the  transition  (or  "flopped")  and  the  others,  by  observing  whether 
the  two  deflections  have  added  or  cancelled.  A  detector  placed  to  measure  the  number 
of  atoms  on  one  or  the  other  of  the  two  paths  will  record  a  resonance  curve  as  the 
frequency  is  varied.  If  the  rf  field  is  applied  uniformly  over  a  given  length,  a  single 
peak  is  obtained.  If  two  separated  fields  are  used  as  in  Ramsey's  method,  a  typical 
interference  pattern  is  obtained  whose  spacing  depends  on  the  velocity  of  the  atom.  If  a 
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Fig.  VIII- 1 

(a)  Schematic  diagram  of  apparatus.  The  apparatus  actually  operates  in  a 
vertical  position  with  the  top  plate  up.  (b)  Cross  sectional  view  of  deflecting 
magnets,  (c)  Schematic  view  of  detector  system.  This  assembly  goes  in  the 
space  at  the  top  of  the  can  in  (a). 


range  of  velocities  is  used,  only  the  central  peak  is  reinforced;  the  side  peaks  tend  to 
average  out.  The  use  of  a  surface  ionization  detector  (hot  wire)  serves  to  change  neutral 
atoms  into  charged  ions  which  may  be  detected  electrically.  With  an  electron  multiplier 
as  amplifier,  the  arrival  of  individual  atoms  can  be  recorded. 

The  general  arrangement  of  the  beam  apparatus  is  indicated  in  the  simplified  form 
of  Fig.  VIII- 1.  The  apparatus  is  housed  in  a  stainless  steel  can  6  feet  long  and  10 
inches  in  diameter  into  which  it  is  lowered  vertically.  All  connections  are  brought  out 
through  the  top  plate. 

The  oven  produces  two  narrow  beams  of  cesium  atoms,  each  of  initial  cross  section 
0.  0Z0  inch  by  3/8  inch  with  a  total  divergence  in  angle  of  less  than  1°.  Details  of  the 
oven  construction  and  performance  are  given  in  the  following  section,  but  with  a  total 
emission  of  10~6  gm  of  cesium  per  day  (5  •  1011  atoms/ sec)  it  was  expected  to  provide 
107  atoms/sec  reaching  the  detector,  of  which  one-eighth  can  provide  the  desired  transi¬ 
tion. 

The  A  and  B  magnets,  which  are  identical,  produce  a  transverse  inhomogeneous 
magnetic  field.  With  an  Alnico  core,  and  circular  section  pole  pieces  (Fig.  VIII- 1(a)), 

3^  field  strength  of  6000  gauss  and  a  gradient  of  more  than  10,  000  gauss/cm  is  obtained 
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in  a  gap  of  l/8  inch.  The  present  magnet  length  is  10  inches,  which  produces  a  deflec¬ 
tion  of  0.02  5  radian  on  atoms  of  the  most  probable  velocity.  This  is  excessive;  the 
magnet  length  will  be  reduced  in  the  future  to  3  l/Z  inches  and  the  gap  increased  to  l/4 
inch.  Magnetizing  coils  are  incorporated  so  that  the  magnets  can  be  set  to  the  desired 
strength.  The  rf  field  system  is  designed  to  use  the  Ramsey  method  of  separated  fields. 
For  each  beam,  two  rectangular  cavities  each  a  half  wavelength  long  and  spaced  66  cm 
apart,  produce  a  transverse  rf  field  (parallel  to  the  A  and  B  fields),  using  the  transverse 
component  of  the  TE^  field  pattern.  The  cavities  are  fed  symmetrically  with  a  wave¬ 
guide  feed.  To  get  a  symmetrical  resonance  curve  it  is  essential  that  the  phase  of  the 
fields  in  the  two  cavities  should  be  the  same.  This  is  obtained  by  providing  a  second 
waveguide  running  between  the  cavities  which  is  loosely  coupled  to  the  cavities  by  the 
slots  through  which  the  cesium  beam  passes.  If  a  probe  is  inserted  at  its  middle  point, 
a  null  signal  is  obtained  if  the  phase  and  amplitude  in  the  cavities  are  identical.  The 
phase  adjustment  is  made  by  a  tuning  screw  in  each  cavity.  This  measurement  can  be 
checked  during  normal  operation,  but  in  general  an  initial  adjustment  is  adequate. 

The  steady  field  in  the  rf  region  is  made  up  of  the  earth's  field  and  the  leakage  field 
from  the  A  and  B  magnets,  with  provision  for  some  control  by  small  coils  near  the 
cavities.  It  has  been  found  that  the  maximum  field  is  below  1  gauss.  In  the  final  form 
some  magnetic  screening  will  be  provided  to  reduce  the  contribution  of  the  earth's  field, 
and  the  desired  field  will  be  provided  by  a  suitable  coil  system. 

The  detector  is  a  0.  040-inch  tungsten  ribbon.  Cesium  atoms  falling  on  this  are 
ionized  and  are  then  pulled  off  and  deflected  electrostatically  to  one  or  the  other  of  two 
electron  multipliers,  each  having  a  gain  of  about  3  X  10^. 

The  apparatus  is  evacuated  by  a  2  l/2-inch  oil  diffusion  pump;  with  liquid  air  traps 

_7 

a  pressure  of  5  *10  mm  is  regularly  maintained. 

A  considerable  amount  of  apparatus,  which  is  described  briefly  below,  has  been 
constructed  for  use  with  the  beam  tube. 

The  most  important  item  is  the  signal  source  at  9192  Mc/sec.  This  is  produced  in 

a  silicon  crystal  tripler  driven  by  an  S-band  cavity  oscillator  (l).  The  oscillator,  with 

a  cavity  with  a  Q  of  about  50,  000  and  a  triode  amplifier,  has  proved  to  have  a  short  term 

q 

stability  better  than  1  part  in  10  .  For  convenience  in  the  experimental  operation  of  the 
beam  tube,  the  long  term  frequency  stability  has  been  improved  by  fitting  a  thermostat 
to  the  cavity  of  the  oscillator.  The  thermostat  uses  a  temperature  sensitive  resistance 
wound  over  the  whole  cylindrical  surface  of  the  cavity.  This  is  connected  in  a  sensitive 
bridge  circuit  that  provides  a  continuous  control  of  the  power  supplied  to  a  heater  winding 
outside  the  cavity,  giving  rapid  and  precise  control  of  temperature.  Changes  of  the 
thermostat  setting  provide  a  convenient  means  of  frequency  adjustment;  they  also 
provide  a  convenient  method  of  producing  a  frequency  changing  linearly  with  time. 

Other  facilities  available  include  an  automatic  frequency  sweep  circuit  that  sweeps 
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the  oscillator  over  a  predetermined  range  by  switching  the  oscillator  cavity  heating 
power  between  two  levels.  In  addition,  provision  is  made  for  the  time  of  the  occurrence 
of  the  resonance  curve  (or  any  other  event,  such  as  a  frequency  marker  from  a  crystal 
standard)  in  the  sweep  to  control  automatically  the  frequency  limit  desired  so  that  it 
continues  to  cover  the  required  range. 

The  frequency  of  the  S-band  oscillator  can  be  compared  with  the  laboratory  crystal 
standard,  using  the  beat  frequency  between  the  S-band  signal  (3064  Mc/sec)  and  the  17th 
harmonic  of  180  Mc/sec  (3060  Mc/sec)  derived  from  a  5  Mc/sec  crystal.  An  audio  fre¬ 
quency  discriminator  (with  its  peaks  located  at  1000  and  1500  cps)  has  proved  very 
convenient  for  recording  small  frequency  variations.  If  a  small  amount  of  frequency 
modulation  at  about  30  cps  is  applied  to  the  S-band  oscillator,  the  beam  current  has  a 
30  cps  component  which  is  proportional  to  the  derivative  of  the  resonance  curve.  When 
this  is  applied  to  a  phase  sensitive  detector,  an  output  having  the  characteristic  of  a 
frequency  discriminator  is  obtained  with  a  zero  value  at  the  resonant  frequency. 

A  Sanborn  two-channel  recorder  has  been  used  for  recording  the  data  and  has  proved 
of  great  value  in  a  number  of  different  types  of  measurement. 

As  soon  as  the  apparatus  was  assembled  with  magnets  of  the  quoted  specifications, 
resonance  curves  of  the  Ramsey  pattern  were  obtained.  Detailed  measurements  are 
much  too  numerous  to  quote  but  the  important  conclusions  are  as  follows. 

Resonance  curves  of  the  Ramsey  pattern  with  good  signal-to-noise  ratios  were 
obtained  (Fig.  VIII-2).  The  oscillator  proved  to  have  more  than  adequate  short-term 
stability  for  the  purpose,  and  can  produce  more  than  enough  power  at  9192  Me/ sec.  The 
oscillator  was  successfully  locked  to  the  beam  frequency,  by  using  the  beam  output  as 
an  error  signal  for  frequency  control.  With  a  very  simple  control  system,  a  peak  fre¬ 
quency  deviation  of  +1  part  109  was  reached  and  the  mean  over  longer  periods  would  be 
much  smaller. 

Probably  of  greatest  importance  is  the  fact  that  the  whole  apparatus  worked  very 
much  according  to  plan.  The  oven  has  run  for  two  extended  periods,  one  of  6  weeks. 


Resonance  curve;  axial  layout. 


Fig.  VIII- 3 

Discriminator  curve;  offset  layout. 
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and  in  each  instance  the  run  was  ended  because  the  apparatus  had  to  be  opened  to  allow 
other  unrelated  changes  to  be  made. 

The  main  feature  requiring  modification  was  that  the  linewidth  obtained,  of  300  cps, 
corresponds  to  using  atoms  of  twice  the  most  probable  velocity.  If  the  magnet  strength 
was  reduced  to  allow  the  use  of  lower  velocity  atoms,  the  background  of  "unflopped" 
atoms,  already  high,  was  further  increased  because  of  certain  features  of  the  geometry 
of  the  apparatus.  For  a  number  of  reasons,  a  more  effective  arrangement  is  a  layout 
in  which  the  oven  is  offset  to  one  side  of  the  axis  and  the  detector  is  offset  to  the  other. 
This  gives  a  loss  of  half  the  available  atoms  but  can  give  a  much  lower  background.  An 
experiment  with  an  excessive  amount  of  displacement  to  suit  the  large  magnet  deflections 
available  gave  the  expected  results,  shown  in  Fig.  VIII-3.  This  curve  is  an  example  of 
the  discriminator  type  of  characteristic  obtained  by  applying  a  small  amount  of  frequency 
modulation  to  the  oscillator.  A  linewidth  of  2  00  cps  with  negligible  background  and  very 
good  signal- to-noise  ratio  was  obtained.  The  total  available  signal  can  be  increased  by 
using  the  optimum  deflection  angle.  The  large  number  of  peaks  was  expected  from  the 
narrow  velocity  range  used  in  this  experiment.  Information  available  from  the  present 
tests  provides  the  essential  information  for  the  engineering  design  of  molecular  beam 
tubes  of  various  sorts.  The  tube  construction  is  simple  and  has  proved  very  satis¬ 
factory.  For  the  immediate  future,  the  beam  apparatus  will  be  set  up  to  provide  a  fre¬ 
quency  standard  by  direct  reference  to  the  cesium  resonance  line  with  an  accuracy  of  a 

9 

single  observation  of  about  1  in  10  .  As  soon  as  any  greater  precision  is  required,  the 
necessary  circuits  for  locking  the  oscillator  continuously  will  be  provided. 

J.  R.  Zacharias,  J.  G.  Yates 
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B.  WELL-COLLIMATED  ATOMIC  BEAM  OVENS 

To  make  possible  the  development  of  atomic  beam  clocks  which  would  run  for  long 
periods  of  time  without  maintenance,  it  was  necessary  to  develop  atomic  beam  sources 
(ovens)  that  would  conserve  as  much  material  as  possible  without  reducing,  significantly, 
the  usable  number  of  atoms.  This  report  describes  the  construction  and  test  of  such 
ovens. 

The  number  of  atoms  per  second  that  will  leave  an  oven  canal  and  reach  a  distant 
detector  which  subtends  an  infinitesimal  solid  angle  located  on  the  axis  of  the  slit  is 
dependent  only  on  the  cross  sectional  area  of  the  slit  and  on  the  pressure  inside  the  oven. 
(In  this  discussion,  the  slit  is  assumed  to  be  some  sort  of  a  closed  cylindrical  surface; 
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this  does  not  restrict  the  argument.)  On  the  other  hand,  the  total  number  of  atoms  per 
second  that  leave  the  slit  is  determined  by  the  solid  angle  which  one  end  of  the  slit  sub¬ 
tends  as  seen  from  the  other  end.  This  is  so  because  most  of  the  atoms  entering  the 
long  canal  will  hit  the  walls  on  the  oven  end  of  the  slit,  and  the  smaller  this  solid  angle, 
the  greater  the  probability  that  atoms  which  hit  the  walls  of  the  slit  will  go  back  into  the 
oven  instead  of  going  out  the  far  end  of  the  slit. 

Thus,  to  save  material  one  is  led  to  the  use  of  a  cylindrical  slit  whose  length  is  large 
relative  to  its  diameter.  However,  an  upper  limit  is  placed  on  the  length  by  the  condi¬ 
tion  that  the  length  be  short  compared  with  the  mean  free  path  for  atoms  in  the  oven 
exit  canal.  If  this  condition  is  not  fulfilled,  then  atoms  will  collide  with  each  other  and 
the  point  at  which  such  collisions  are  most  probable  will  serve  as  a  new  source  —  so 
that  the  entire  length  of  the  slit  will  not  be  contributing  to  the  material  savings.  Since 
the  forward  intensity  is  proportional  to  a  positive  power  of  the  pressure  in  the  oven  and 
since  one  is  usually  interested  in  obtaining  as  great  a  forward  intensity  as  possible,  one 
obviously  wants  to  operate  the  oven  at  as  high  a  pressure  as  possible.  But  the  mean 
free  path  is  inversely  proportional  to  the  pressure.  Therefore,  one  is  forced  to  make 
the  slit  as  short  as  possible.  This  in  turn  requires  the  slit  diameter  to  be  very  small. 
To  obtain  usable  forward  intensities  one  must  use  a  large  number  of  such  slits  with 
their  axes  parallel  and  with  as  thin  boundary  walls  as  possible. 

The  oven  slits  now  in  use  were  constructed  by  rolling  0.001  inch  thick  nickel  foil 
between  two  interlocking  bronze  rolls.  The  width  of  sheet  having  the  best  rolling 
properties  was  found  to  be  0.5  inch.  This  implies  a  canal  length  of  0.  5  inch.  The  length 
of  the  sheet  is  completely  arbitrary.  The  resulting  sheet  (nicknamed  "crinkly  sheet") 
looks  very  much  like  miniature  corrugated  iron.  The  corrugations  are  about  0.002  inch 
deep  and  about  0.  007  inch  from  peak  to  peak.  The  sheet  can  be  cut  to  the  desired  length 
with  a  pair  of  kitchen  scissors.  By  alternately  stacking  these  "crinkly  sheets"  and 
sheets  of  flat  0.001-inch  nickel  foil,  one  can  obtain  canal  assemblies  of  arbitrary 
dimensions. 

For  the  small  clock,  the  dimensions  of  the  slit  assembly  are  0.  02  5  inch  by  1. 0  inch, 
and  the  length  of  each  slit  is  0.  5  inch. 

In  the  tests  analyzing  the  directivity  of  this  oven,  a  charge  of  cesium  carbonate  and 
potassium  was  placed  in  the  oven.  The  beam  of  cesium  atoms  obtained  on  heating  the 
oven  was  detected  using  a  surface  ionization  detector  (hot  wire).  The  oven  temperature 
was  measured  with  an  iron-constantan  thermocouple  bolted  on  the  outside  of  the  oven. 
The  surface  ionization  detector  was  a  heated  tungsten  wire,  0.005  inch  in  diameter, 
located  12  inches  from  the  oven  slit.  The  useful  length  of  the  hot  wire  was  about  two 
inches.  The  oven  was  mounted  so  that  it  could  be  turned  about  a  vertical  axis  parallel 
to  the  long  dimension  (one  inch)  of  the  slit  assembly  and  perpendicular  to  the  axes  of 
the  slits.  The  angle  through  which  the  oven  was  turned  was  read  from  a  micrometer 
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Fig.  VIII- 4 

Plot  of  beam  intensity  N  vs.  oven  angle  (0).  N  is  plotted  in  units  of  Nmax, 
the  intensity  in  the  forward  direction.  The  oven  temperature  was  341  °K. 

n 

The  value  of  N  was  1.1  X  10  atoms/sec  hitting  the  hot  wire, 
max  ' 
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Fig.  VIII- 5 

Plot  of  beam  intensity  N  vs.  oven  angle  (0).  N  is  plotted  in  units  of  Nmax, 
the  intensity  in  the  forward  direction.  The  oven  temperature  was  417°K. 

g 

The  value  of  Nmax  was  1.9  X  10  atoms/sec  hitting  the  hot  wire. 
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screw  which  actuated  the  turning  lever.  The  number  of  cesium  atoms  which  hit  the 
hot  wire  in  a  second  was  measured  as  a  function  of  oven  temperature  and  angle.  Two  of 
the  curves  which  were  obtained  are  Fig.  VIII-4  and  Fig.  VIII- 5.  In  both  of  these  curves, 
the  ordinate  of  the  curve  has  been  normalized  to  unity  at  an  angle  of  zero  radian.  Since 
a  linear  assembly  of  slits  was  used  and  the  detector  was  also  linear,  one  should 
interpret  these  curves  as  the  measurement  of  angular  spread  in  one  direction  only.  Since 
each  slit  is  approximately  circular  in  cross  section,  the  curve  of  angular  spread  for  one 
slit  would  be  approximately  the  figure  of  revolution  obtained  by  rotating  Figs.  VIII-4  and 
VIII- 5  about  the  0  =  0  axis.  From  these  figures  we  see  that,  in  so  far  as  one  dimension 
is  concerned,  the  amount  of  material  within  the  half  intensity  angle  is  about  10  per  cent 
of  the  total  amount  of  material  at  a  temperature  of  341  °K.  However,  from  Fig.  VIII-5, 
the  amount  of  material  in  the  central  peak  is  only  about  5  per  cent  of  the  total  for  a 
temperature  of  417°K.  The  material  savings  in  the  other  direction  will  depend  on  the 
angle  which  the  long  dimension  of  the  detector  subtends  at  the  source.  With  Fig.  VIII-4 
as  a  guide,  we  see  that  if  this  angle  is  0.010  radian,  about  10  per  cent  of  the  beam  in 
that  direction  will  be  useful  at  a  temperature  of  341  °K.  In  the  experiments  described 
above,  the  detector  subtended  an  angle  of  approximately  0.2  radian  so  that  very  little 
material  was  wasted  in  this  direction. 

The  major  conclusions  of  the  experiments  that  have  been  completed  are: 

1.  A  definite  material  saving  results  from  the  use  of  collimated  slits.  The  curves 
obtained  show  this  when  compared  with  the  cosine  distribution  for  N/N  vs.  0  which 
would  be  obtained  for  an  infinitely  short  slit  or  when  compared  with  the  curves  given 
by  Davis  (l).  At  low  temperatures  (see  below)  this  material  saving  seems  to  be  roughly 
that  expected  on  the  basis  of  considerations  of  solid  angle  alone. 

2.  The  source  is  equivalent  to  an  assembly  of  about  300  canals  each  giving  about 
10°  atoms/sec  at  341°K  and  10  atoms/sec  at  417°K  — all  going  into  the  solid  angle 
within  the  cone  0  =  0.010  radian.  The  width  at  half  intensity  for  a  single  canal  is  about 
0.010  radian  for  both  temperatures.  The  proportion  of  atoms  going  outside  the  cone 
0  =  0.010  radian  is  a  function  of  temperature. 

3.  The  material  saving  is  a  function  of  temperature.  The  higher  the  temperature, 
the  greater  the  ratio  of  total  material  leaving  the  oven  to  "forward  intensity."  (That  is, 
the  "wings"  of  the  curves  in  Figs.  VIII-4  and  VIII-5  become  larger  as  the  temperature 
is  increased.)  This  temperature  dependence  is  probably  due  to  scattering  in  the  slit  (1). 

J.  R.  Zacharias,  R.  D.  Haun,  Jr. 
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C.  HYPERFINE  STRUCTURE  ANOMALIES  IN  CESIUM  RADIOISOTOPES 

The  ratios  of  the  nuclear  g  values  of  the  four  cesium  isotopes  —  Cs  (stable),  and 

1 34  135137  6 

Cs  ,  Cs  ,  and  Cs  (radioactive  with  half  lives  of  2.3,  3  X  10  ,  and  37  years, 

respectively)  —  have  been  measured  by  the  atomic  beam  method,  by  observing  the  doublet 

transitions  AF  =  0,  m  =  -1  —  -2  in  the  odd-even  isotopes,  and  m  =  -3/2  —  m  =  -5/2  in 

the  odd-odd  isotope  in  fields  of  8300  to  9500  gauss.  The  results  are: 


g135 

g133 


1. 05817  +  0.  00006 


gl  37 
gl  35 


1. 04004  +  0. 00005 


g137 
g  1 33 


1.  10047  +  0. 00008 


gl  34 
g133 


1. 01447  +  0. 0003 


When  these  results  are  compared  with  those  obtained  from  the  ratios  of  the  Av's  (l) 
by  means  of  the  Fermi  formula  (2), 

»  8  ,2/r>x2I+l 

Av  =3  upj  (0)  — | - 


derived  on  the  assumption  that  the  nuclear  magnetic  moment  is  a  point  dipole,  it  is  found 
(3)  that  there  is  a  discrepancy  A.  The  experimental  values  are 

A( 135- 133)  =  +0. 034  +  0. 006  per  cent 

A(  1 37-  1 35)  =  -0.  022  +  0.  005  per  cent 

a(  137- 133)  =40.  006  ±  0.  008  per  cent 

A(l34-133)  =  +0.  17  +  0.  03  per  cent 

When  these  values  are  compared  with  the  results  of  a  detailed  treatment  of  the  hfs  inter¬ 
action  according  to  reference  3,  in  which  a  departure  from  a  nuclear  point  dipole  theory 
is  considered,  one  finds,  on  the  basis  of  a  strict  independent  particle  model  with 
quenched  intrinsic  magnetic  moments,  agreement  with  experiment  for  the  three  isotopes 

Cs1  ,  Cs1  ,  and  Cs1  .  Cs  has  an  hfs  anomaly  which  has  a  sign  opposite  to  the 

13  5 

one  predicted  with  respect  to  Cs  ,  as  well  as  one  smaller  by  a  factor  of  10  with  respect 

133  137 

to  Cs  .  This  result  may  be  related  to  the  fact  that  Cs  is  a  magic  nucleus  with  a 

closed  shell  of  82  neutrons.  A  detailed  description  of  the  experiment  and  theory  will  be 

presented  in  a  separate  report. 

H.  H.  Stroke,  V.  Jaccarino,  D.  S.  Edmonds,  Jr.,  R.  Weiss 
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D.  NUCLEAR  MAGNETIC  OCTUPOLE  MOMENTS  OF  THE 
STABLE  GALLIUM  ISOTOPES 

2 

The  three  zero-field  hyperfine  structure  intervals  of  the  metastable  P,/?  state  of 

6971  ^  ^ 

Ga  7  and  Ga  have  been  remeasured  (l)  to  higher  precision  to  establish  the  existence 

of  a  nuclear  magnetic  octupole  moment. 

The  atomic  beam  magnetic  resonance  method  was  employed;  the  transitions  were 
induced  by  the  Ramsey  technique  (2)  of  separated  oscillating  fields.  One  oscillating  field 
was  phase -modulated  at  93  cps  with  respect  to  the  other  so  that  the  oscillating  fields 
were  alternately  in  phase  and  180°  out  of  phase.  The  alternately  in-phase  and  out-of¬ 
phase  Ramsey  patterns  produced  were  electronically  subtracted  by  a  synchronous 
detector  and  displayed  on  an  oscilloscope  (see  Fig.  VIII- 6). 

The  intervals  were  measured  in  several  external  fields  ranging  from  0.2  gauss  to 
1  gauss,  extrapolated  to  zero  field,  and  corrected  (3)  for  the  perturbing  effects  of  the 


Fig.  VIII- 6 

Typical  oscilloscope  photograph  showing  the  quadratic  Zeeman  splitting  of 
the  (F  =  2 ,  =  0)  — -  ( 1 ,  0)  and  (2 ,  1 )  — -  ( 1 ,  1 )  transitions  in  an  external  mag- 

netic  field  of  0.4  gauss.  The  Ramsey  technique  (2)  of  separated  oscillating 
fields  was  used.  The  two  blank  portions  of  the  trace  calibrate  the  picture 
in  frequency.  The  width  at  half  maximum  of  each  central  peak  is  500  cps. 
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Table  I 


Measured  and  Corrected  HFS  Intervals  of  the 
Metastable  ^3/2  State  of  Ga^  and  Ga71. 

Ga^  (Mc/sec)  Ga71 


(Me/  sec) 


Measured 


II 

0 

1 

II 

H— « 

128.27730  + 

0. 00020 

203. 04340  + 

0.  00020 

II 

t 

r— H 

11 

319. 06706  + 

0. 00020 

445. 46960  +  0. 00020 

F  =  2  — -F  =  3 

634. 90183  + 

0. 00020 

766. 69580  + 

0. 00020 

Corrected 

F  =  0  —  F  =  1 

128.27650  + 

0. 00030 

203. 04133  + 

0. 00040 

F  =  1  —  F  =  2 

319. 06373  + 

0. 00050 

445.46566  + 

0. 00060 

F  =  2  —  F  =  3 

634. 90597  + 

0. 00060 

766. 70181  + 

0. 00080 

neighboring  &rounc*  state  as  well  as  the  admixture  of  (4s)  (4p)  (5s)  electronic  con¬ 

figuration  (see  Table  I). 

The  magnetic  dipole,  electric  quadrupole,  and  magnetic  octupole  interaction 
constants  (4),  a,  b,  and  c,  are  calculated  from  the  corrected  intervals  to  be: 


Ga 


69 


Ga 


71 


a:  190.  794280  +  0.  000150  Mc/sec  242 . 433950  +  0.  000200  Mc/sec 

b:  62.  522470  ±  0.  000300  Mc/sec  39. 399040  +  0. 000400  Mc/sec 

c:  84+6  cps  115+7  cps 

where  the  uncertainty  is  the  root- sum- square  of  the  uncertainties  in  each  of  the  terms 
of  the  equations  that  yield  the  interaction  constants  in  terms  of  the  intervals. 

The  data  cannot  be  explained  by  the  mechanism  of  off-diagonal  terms  in  the  inter¬ 
action  matrix  giving  rise  to  octupole-like  (I*  J)^  terms,  since  (ab)^/(ab)7 1  and  (b^9)2/ 
712  6971 

(b  )  are  greater  than  unity,  and  c  /c  is  less  than  unity.  To  within  the  quoted  uncer¬ 
tainty  a69/a?1  =  c69/c?1. 

Schwartz  has  evaluated  the  nuclear  octupole  moment,  ft,  and  finds 

69  -24  2 

£2  =  (0.  107  +  0.  02)  X  10  nuclear  magneton  cm 


£271  =  (0.  146  +  0.  02)  X  10  2  nuclear  magneton  cm2 

Here,  as  in  the  corrected  intervals,  the  quoted  uncertainty  includes  the  experimental 
error  and  the  uncertainty  in  theoretical  evaluation. 

R.  T.  Daly,  Jr.,  J.  H.  Holloway 
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E.  HYPERFINE  STRUCTURE  OF  HALOGENS 

V.  Jaccarino  and  the  author  have  remeasured  the  hyperfine  structure  of  the  stable 
bromine  isotopes  to  higher  precision  and  have  found  an  octupole  interaction  that  is  well 
outside  the  experimental  uncertainties.  Further  experimental  work  must  be  done  before 
the  final  results  are  published. 

129 

The  hyperfine  structure  of  the  long-lived  radioactive  isotope  I  is  to  be  investi- 
gated.  Since  a  large  octupole  interaction  was  found  in  the  hfs  of  I  previously  studied 
at  this  laboratory,  there  is  considerable  interest  in  the  proposed  experiment.  Apparatus 
for  handling  the  small  sample  available  and  for  improving  the  signal-to-noise  ratio  is 
being  designed. 

J.  G.  King 
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IX.  MULTIPATH  TRANSMISSION 


Prof.  J.  B.  Wiesner 
E.  J.  Baghdady 


NARROW -BAND  LIMITING 


A  study  of  the  behavior  of  the  plots  relating  the  minimum  requirement  in  discrimi¬ 
nator  bandwidth  (to  reject  two-path  FM  interference)  to  the  ideal  limiter  filter  bandwidth 
used  (Quarterly  Progress  Report,  April  15,  1954,  Fig.  VIII-2)  sheds  some  interesting 
light  on  the  nature  of  a  set  of  curves  that  may  be  associated  with  those  plots.  We  refer 
here  to  the  smooth  curves  that  may  be  drawn  to  connect  the  values  of  minimum  required 
(BW)disc/(BW)^  at  the  odd  integral  values  of  (^W)^m/(BW)^  used.  These  curves  are 
identified  as  the  smooth  curves  shown  in  Fig.  IX- 1  superimposed  on  the  calculated 
broken-curve  plots.  We  will  call  these  smooth  curves  the  envelopes  of  the  calculated 
plots  since  Fig.  IX-1  shows  that  as  (BW)^m/(BW).^  grows  large,  the  broken-curve  plots 
tend  to  be  bounded  by  the  smooth  curves. 

If  we  make  the  definitions 

(BW)if 


<BW>disc 

y  =  -  at  odd  values  of  x 

(BW)if 

<BW>disc 

=  "envelope"  value  of  - 

(BW)if 


then  the  semilogarithmic  plots  of  Fig.  IX-2  show  that 


In 


1  +  a 
1  -  a 


is  linearly  related  to  x.  We  recall  that  a  is  the  ratio  of  the  undesired  to  the  desired 
signal  amplitude  at  the  input  to  the  ideal  limiter.  Simple  calculations  show  that  the  plots 
of  Fig.  IX-2  lead  to  the  expression 


y(a,  x) 


1  +  a 
1  -  a 


1  -  ^(a)  e 


where  values  of  k(a)  and  t,(a)  are  plotted  in  Figs.  IX-3(a)  and  (b).  Figure  IX-3(a)  shows 
k(a)  plotted  against  In  a  and  shows  that  the  plotted  points  fit  rather  closely  into  a  straight 
line  given  by 

k(a)  =  -0.  395  In  a 
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(a) 


Fig.  IX- 3 

Calculated  values  of  the  constants  in  Eq.  1,  resulting  from  the  plots  of 
Fig.  IX-2 .  The  deviations  of  the  plotted  points  from  the  straight  line  values 
may  be  attributed  to  small  cumulative  errors  in  computation. 
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The  plot  of  t,(a)  shown  is  based  on  k(a)  as  given  by  Eq.  2,  with  the  reasonable  assumption 
that  the  small  deviations  of  the  calculated  values  of  k(a)  from  the  values  given  by  Eq.  2 
may  be  attributed  to  small  cumulative  errors  in  the  computations.  Figure  IX- 3(b)  shows 
that  the  plotted  values  of  f(a)  fit  into  a  straight  line  given  by 

C(a)  =  0.  30  a  +  0.  44  (3) 

Equation  1  may  be  normalized  into  the  form 


4>(a)  =  y(a,  x) 


1  +  a 

1  -  a 


=  1  -  C(a)  e 


k(a)x 


(4) 


which  is  plotted  in  Fig.  IX-4.  Clearly,  the  second  term  in  the  right  member  of  Eq.  4 
gives  the  fractional  reduction  in  the  minimum  required  discriminator  bandwidth  that 
results  from  passing  the  resultant  two-path  signal  through  an  ideal  limiter  whose  band¬ 
width  is  an  odd  multiple,  x,  of  the  bandwidth  of  the  intermediate-frequency  section. 

Now,  at  the  odd  values  of  x,  it  is  recalled  that  (l) 


(BW)if 


[Aw] 

=  2  - +  1 

(BW)if 


(5) 


where  [acj]  is  the  magnitude  of  the  maximum  deviation  in  the  instantaneous  frequency 
of  the  resultant  signal  at  the  output  of  the  ideal  limiter  filter  from  the  frequency 
of  the  desired  carrier.  It  is  clear  from  Fig.  IX- 1  that  Eq.  5  also  gives  the  value 
of  (BW)^-gc/(BW).|.  over  the  horizontal  segment  following  the  odd  value  of  x  for  which 
[aco]  is  calculated.  A  combination  of  Eq.  1  and  Eq.  5  yields 


[au] 

(BW)if 


=  j(y-!) 


a 

1  -  a 


1  /I  +  a 

2  U  -  a 


e 


k(a)x 


(6) 


Further  normalization  yields 


1  +  a  r  /  \  -k(a)x 
— ;(a)e 


(7) 


The  second  term  on  the  right  of  Eq.  7  yields  the  fractional  amount  by  which  the  magni¬ 
tude  of  the  maximum  deviation  in  instantaneous  frequency  has  been  reduced  by  passing 
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(0)/f\ 
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Fig.  IX-4 

Normalized  plots  of  the  envelopes  of  the  broken  curves  in  Fig.  IX- 1  .  »|<(a) 
equals  the  normalized  minimum  required  discriminator  bandwidths  at  the  odd 
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the  resultant  two-path  signal  through  an  ideal  limiter  whose  bandwidth  is  an  odd 
multiple,  x,  of  the  bandwidth  of  the  intermediate  frequency. 

Equation  1  is  a  useful  supplement  to  the  calculated  plots  of  Fig.  IX- 1.  It  can  be  used 
to  extrapolate  the  correct  values  of  the  minimum  required  discriminator  bandwidth  when 
the  values  of  limiter  bandwidth  used  are  odd  integral  multiples,  x,  of  the  i-f  bandwidth, 
that  lie  beyond  the  range  covered  by  the  calculated  plots.  At  other  values  of  x,  Eq.  1 
will  yield  an  approximate  answer.  Equation  1  may  also  be  used  for  the  same  purpose 
when  the  capture  ratio,  a,  is  required  to  assume  a  value  that  lies  in  the  range  0.7^ 
a  0.  95,  but  is  not  covered  by  the  plots  of  Fig.  IX- 1.  The  indicated  restriction  in  the 
range  is  given  because  the  behavior  of  t,(a)  outside  this  range  is  not  completely  obvious 
from  the  plot  of  Fig.  IX-3(b). 

E.  J.  Baghdady 
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Prof.  Y.  W.  Lee  S.  G.  Margolis  S.  D.  Pezaris 

A.  G.  Bose  C.  E.  McGinnis  W.  B.  Smith 

J.  Y.  Hayase  R.  E.  Wernikoff 

A.  A  THEOREM  CONCERNING  NOISE  FIGURES 

An  expression  is  obtained  for  the  greatest  lower  bound  of  the  noise  figure  of  a  system 

consisting  of  n  identical  tubes  whose  plate  signals  are  directly  added  and  whose  inputs 

are  obtained  through  a  linear  passive  coupling  network  from  a  source  having  an  internal 

impedance  R  .  It  is  found  that  the  optimum  attainable  noise  figure  of  such  a  system  is 
s 

equal  to  the  optimum  noise  figure  of  a  circuit  consisting  of  one  such  tube  and  an  appro¬ 
priate  coupling  network. 

Introduction 

The  electron  tube,  while  providing  amplification,  is  a  source  of  noise.  Thus,  for 
the  detection  of  weak  signals  it  is  advantageous  to  design  the  circuitry  associated  with 
the  amplifier  to  maximize  the  signal-to-noise  ratio  at  the  output  of  the  amplifier  when 
it  is  fed  from  a  given  source.  That  is,  it  is  advantageous  to  minimize  the  noise  figure 
of  the  circuit.  In  the  past,  attempts  have  been  made  to  obtain  smaller  noise  figures  than 
obtainable  with  one  tube  by  using  circuits  involving  more  tubes;  the  plate  signal  of  each 
tube  was  added  directly  to  the  plate  signals  of  the  other  tubes  to  form  the  output  of  the 
circuit.  Such  attempts  have  been  motivated  by  the  hope  that  smaller  noise  figures  could 
be  obtained  by  taking  advantage  of  the  incoherency  of  noise  voltages  from  different  noise 
sources.  (The  signal  at  the  output  of  each  tube  adds  coherently  while  the  noise  arising 
from  different  tubes  adds  incoherently.)  It  is  the  purpose  of  this  work  to  prove  that  the 
noise  figure  of  a  circuit  involving  n  identical  tubes,  connected  as  indicated  above,  is  at 
best  equal  to  the  optimum  noise  figure  attainable  with  the  use  of  only  one  tube  and  an 
appropriate  coupling  network. 

Proof 

The  noise  figure  calculations  are  based  on  a  unit  bandwidth  centered  about  a  fre¬ 
quency  f .  It  is  assumed  that  each  tube  may  be  represented,  as  far  as  noise  figure  cal¬ 
culations  are  concerned,  by  a  shunt  resistance  R 

J  g 

and  a  series  equivalent  resistance  R  followed  by 

an  ideal  noiseless  tube.  According  to  this  assump¬ 
tion,  the  equivalent  circuit  is  shown  in  Fig.  X-l. 
Associated  with  each  resistance  is  a  series  thermal 
noise  voltage  source.  For  a  resistance  R  the 
open  circuit  mean  square  voltage  per  unit  band¬ 
width  is  4KTpR,  where:  K  =  Boltzmann's  constant; 


Fig.  X-l 

Noise  equivalent  circuit  of  con¬ 
ventional  amplifier  tube. 
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Fig.  X-2 

Equivalent  circuit  (for  noise  figure  calculation)  of  a 
source  coupled  through  a  general  coupling  network  to 
n  tubes  whose  plate  signals  are  directly  added. 


T  =  temperature  in  degrees  Kelvin;  p  =  a  positive  factor  used  to  take  into  account  the 
equivalent  noise  temperature  of  the  resistance.  We  shall  call  it  the  equivalent  temper¬ 
ature  factor. 

An  equivalent  circuit  (as  far  as  noise  figure  calculations  are  concerned)  of  a  system 

consisting  of  n  identical  tubes  whose  plate  signals  are  added  directly  and  whose  grid 

signals  are  obtained  through  a  general  linear  passive  coupling  network  from  a  source 

with  impedance  R  is  shown  in  Fig.  X-2.  The  ideal  noiseless  tubes  of  Fig.  X- 1  have 

s 

been  replaced  by  ideal  transformers  whose  turns  ratios  have  conveniently  been  chosen 
to  be  unity.  There  is  no  loss  of  generality  in  this  choice  since  a  change  in  the  output 
level  affects  both  signal  and  noise  and  thus  leaves  the  signal -to -noise  ratio  unaltered. 
The  circuit  shown  in  Fig.  X-2,  with  the  terminals  p-p  open  circuited,  has  the  same 
noise  figure  as  the  system  of  n  tubes  under  consideration.  The  resistances  R  through 
Rj  represent  all  the  resistors  in  the  coupling  network.  The  sources  in  series  with  each 
resistor  are  the  equivalent  noise  voltage  sources.  The  equivalent  temperature  factor 
of  each  R^  is  p.  The  equivalent  temperature  factor  of  each  of  the  resistors  in  the 
coupling  network  is  arbitrary  and  is  denoted  in  Fig.  X-2  by  the  symbol  p  with  a  sub¬ 
script  corresponding  to  the  resistor  with  which  it  is  associated. 

For  a  linear  system,  the  ratio  of  the  open  circuit  voltage  response  at  terminals  k-k 
to  a  voltage  excitation  at  terminals  j-j  can  be  expressed  in  terms  of  the  short  circuit 
admittance  parameters  as: 

(1) 

The  noise  figure  of  a  system  is  defined  as 


gm  = 

^kk 
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F  =  1  + 


"noise  power"  at  the  output  caused  by  noise  sources  in  the  system 


"noise  power"  at  the  output  caused  by  noise  generated  in  the  source  impedance 

(2) 

The  term  "noise  power"  denotes  the  integral  of  the  power  spectrum  of  the  open  circuit 
voltage  over  a  unit  bandwidth  centered  about  a  frequency  f .  Making  use  of  Eqs.  1  and  2, 
the  noise  figure  for  the  system  shown  in  Fig.  X-2  can  be  written  as 

p  R 

~(lyoil2  +  |y02l2  +  •  •  • +  |yonl2) +  nhz  (paRa  ly<J2  +  •  •  • +  ^  lyoj|2) 


n  R  + 
eq 


F  =  1  +- 


|y. 


oo 


oo 


R 


ly 


os 


s  .  ,2 

1^00 


(3) 

Certainly  the  right-hand  side  of  this  equation  will  not  increase  if  we  assume,  for  the 
purpose  of  calculation,  that  each  of  the  resistors  R  through  R.  in  the  coupling  network 
has  an  equivalent  temperature  factor  equal  to  zero.  Equation  3  then  becomes 

2  .  „  / ■  ,2  .  ,  ,2  .  .  ,  ,2' 


F  ^  1  + 


nReql^oor  +  P  Rg(>0l|2  +  ly02l  +  •**  +  lyoJ  ) 


(4) 


Rs  lyDS  I 


Notice  that  the  equality  sign  holds  in  Eq.  4  if  R  through  R.  in  Eq.  3  are  zero. 

a  J 

Let  gQO  be  the  real  part  of  the  output  admittance  yQO-  As  a  consequence  of  the 


relation  |y  I  ^  g  we  can  write 
i j oo 1  °no 


F  ^  1  + 


oo 


n  Reqgoo  +  P  Rg  ( ly01  l2  +  ly02  I2  +  •••  +  f  yon  1 2 ) 


(5) 


Rslyosl 


Notice  that  the  equality  sign  in  Eq.  5  holds  when  the  condition  stated  under  Eq.  4  holds 
and,  at  the  same  time,  the  output  admittance  is  real. 

In  order  to  reduce  the  number  of  variables  in  this  expression  for  F  (Eq.  5)  let  us 
make  use  of  power  relations  at  the  output  terminals  0-0.  Let  P  be  the  power  absorbed 
by  the  system  when  one  volt  is  applied  at  the  output  terminals  and  all  other  voltage 
sources  are  shorted. 

Then 


p  =  Rslysol2  +  Rg  (lyio 1 2  +  ly2<>l2  +  •••  +  ly: 


no 


+  P. 


(6) 
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where  P.  is  the  power  absorbed  in  the  coupling  circuit  resistors  R  through  R.  and  is 
1  a  3 

therefore  equal  to  or  greater  than  zero.  Furthermore,  for  the  particular  excitation 

mentioned  above,  P  is  just  equal  to  gQO*  Using  this  fact  and  the  relation  y.,  =  y,.  for 

a  bilateral  network  we  obtain  the  equation 

§oo =  Rs|yos|2  +  Rg  (tyoi  i2  +  |y02l 2  +  ••• +  |y0nl2)+pi  <7> 

Using  Eq.  7  in  Eq.  5  we  have 

nR  g2  +pg  ■  p  P, 
eq  &oo  &oo  r  1 

F  ^  1  -  p  + -  (8) 


In  order  to  simplify  this  expression  further  let  us  again  consider  the  system  shown 
in  Fig.  X-2.  We  shall  excite  the  system  at  the  terminals  0-0  with  a  unit  voltage  source 
For  convenience  Fig.  X-2  is  redrawn  in  Fig.  X-3  with  the  resistances  R  referred  to 
the  secondary  of  the  ideal  one-to-one  ratio  transformers.  We  desire  to  find  a  relation 
between  gQo  and  |yos  [•  In  obtaining  this  relation  it  is  convenient  to  make  use  of  the 
following  inequality  (the  proof  of  this  inequality  is  given  in  the  following  subsection): 


(9) 


For  a  one-volt  excitation  at  terminals  0-0  it  is  seen,  from  power  considerations  and 
the  use  of  Eq.  7,  that 


Notice  that  the  equality  sign  holds  if  the  coupling  network  is  symmetrical  with  respect 

2 

to  its  n  outputs.  Substitution  of  the  expression  for  |y  |  ,  obtained  from  Eq.  10,  in 
Eq.  8  yields 


n  R  g2  +  p  ■  p  P, 
eq  &oo  °oo  r  1 

F  ^  1  -  p  +  -  =  A  (11) 

g  -  — 5-  -  P, 

&oo  n  R  1 

g 

The  quantity  A  is  thus  expressed  as  an  explicit  function  of  the  two  variables  gQo  and  P^ 
We  desire  to  minimize  A  with  respect  to  these  two  variables.  It  can  be  seen  from 
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Fig.  X-3 


Fig.  X-4 


Redrawn  with  all  the  noise  voltage  sources  set 

to  zero  and  with  the  resistances  R  referred 

g 

to  the  secondary  of  the  ideal  transformers. 


Example  of  a  circuit  using  n  tubes 
and  having  the  minimum  obtainable 
noise  figure  of  one  tube. 


Eq.  7  that  gQQ  and  are,  in  general,  dependent  quantities.  If  A  is  minimized  with 
respect  to  gQQ  and  P^,  considered  as  independent  variables,  the  minimum  value  thus 
obtained  will  certainly  be  equal  to,  or  less  than,  the  minimum  obtained  when  the  con¬ 
straints  are  imposed.  Further,  it  will  be  seen  that  in  this  particular  case  the  minimums 
are  equal  since  a  coupling  network  corresponding  to  the  minimum  A  will  be  physically 
realizable . 

Proceeding  with  the  minimization,  Eq.  11  yields 


8A  nRg 


P  +  n  R  g2 
eq  &oo 


dP 


^ oo  ~  n  Rg  P1 


>  0 


(12) 


A  justification  of  the  inequality  symbol  in  this  equation  is  necessary.  Referring  to 

Eq.  10  we  see  that  the  denominator  in  Eq.  12  is  always  positive  and  greater  than  zero. 

(It  can  be  zero  only  for  the  trivial  case  in  which  R  |  y  |  is  zero,  a  case  which  clearly 

is  of  no  interest  since  it  implies  either  a  zero  impedance  source  or  no  transmission 

from  input  to  output,  that  is,  y  =0.)  Since  p  is  a  positive  quantity,  it  is  clear  that 

s  o 

the  numerator  of  Eq.  12  is  greater  than  zero  and  that  the  inequality  symbol  is  valid  as 
shown. 

Since  (OA/OPj)  >  0  and  since  P^  ^  0,  the  minimum  value  of  A  occurs  for  P^  =  0. 

We  next  find  the  optimum  gQQ  which  gives  the  minimum  A  for  Pj  =  0. 

n  Req  Soo  +  C 

A(Pj  =  0)  =  1  -  P  +  |  (13) 

1  '  n  Rg  §00 
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2  R 


n  R 


eq 


dA 


eq 


(Pj  =  0) 


R  g  o  2 

g  &oo  n  R  g 
&  g  &oo 


3g 


=  0 


oo 


1  - 


1 


n  R  g 
g  soo. 


(14) 


For  the  same  reasons  discussed  earlier  relative  to  Eq.  12,  the  denominator  of  Eq.  14 
is  always  greater  than  zero.  Thus  by  setting  the  numerator  equal  to  zero  we  obtain  the 
optimum  value  of  gQO- 


g 


1 


oo  ,  n  R 
opt  g 


p  R 


1/2 


1  +  1  + 


g 


R 


eq 


Substituting  this  optimum  value  of  gQQ  into  Eq. 


A 

min 


1  +  2 


1  + 


P  R 


1/2 


g 


R 


eq 


(15) 

13  we  obtain  the  minimum  value  of  A. 


(16) 


The  fact  that  the  value  of  A  obtained  by  this  procedure  is  a  minimum  can  be  verified  by 
considering  Eqs.  10  and  13. 

In  view  of  Eqs.  11  and  15  we  can  write 


F  5s  1  +  2 


P  R, 


1/2 


1  +  1  + 


R 


eq 


(17) 


The  fact  that  the  equality  sign  in  this  inequality  is  valid  is  most  readily  demonstrated  by 

presenting  an  example  in  which  it  holds.  The  construction  of  an  example  is  facilitated 

by  noting  the  sufficient  conditions  for  the  validity  of  the  equality  sign  in  the  relations  of 

Eqs.  4,  5,  and  10.  These  conditions  are  stated  immediately  following  their  respective 

equations.  Such  an  example  is  shown  in  Fig.  X-4.  All  of  the  n  tubes  are  in  parallel, 

and  the  general  coupling  network  is  just  an  ideal  transformer  with  an  appropriate  turns 

ratio.  A  straightforward  calculation  shows  that  this  circuit  has  the  noise  figure  given 

by  Eq.  17  with  the  equality  sign  when  each  tube  has  the  input  resistance  R  (with  an 

g 

equivalent  temperature  factor  p)  and  the  equivalent  shot  noise  resistance  R^. 

The  interesting  result  is  that  the  right-hand  member  of  Eq.  17  is  independent  of  n 

and  is  exactly  the  optimum  noise  figure  attainable  using  a  single  tube.  For  R  /R  «  1, 

e4  & 


-53- 


(x.  STATISTICAL  COMMUNICATION  THEORY) 


Eq.  17  reduces  to  the  more  familiar  form  for  the  noise  figure  of  a  single  stage  (1): 


F  5s  1  +  2 


R 


1/2 


eq 


R 


g 


(18) 


Derivation  of  Inequality  9 


Derivation  of  the  inequality: 

k=  1 


(9) 


Since  the  magnitude  of  the  sum  of  complex  numbers  is  less  than  or  equal  to  the  sum 
of  the  magnitudes  we  can  write 


I 


But 


1 

n 


n 


k=  1 


1 

n 


•<  (n- 1 ) 


n 


k=l 


k  I  (leil  -  lekl)2s  I 

i,  k  k=  1 

i*k 


(19) 


(20) 


Thus,  combining  Eq.  19  and  Eq.  20  we  have  the  desired  result: 


A.  G.  Bose,  S.  D.  Pezaris 


References 

1.  G.  E.  Valley,  Jr.  and  H.  Wallman,  M.I.T.  Radiation  Laboratory  Series,  Vol.  18 
(McGraw-Hill  Book  Company,  Inc.,  New  York,  1948)  p.  640. 


-54- 


(X.  STATISTICAL  COMMUNICATION  THEORY) 


B.  THE  WIENER  THEORY  OF  NONLINEAR  SYSTEMS 


Introduction 

Briefly  stated,  the  objectives  of  Dr.  Wiener's  method  are  to  characterize  nonlinear 
systems  independent  of  their  input  and  to  present  a  method  for  synthesizing  the  systems 
from  their  characterizing  parameters.  An  operator  relating  the  input  and  output  of  the 
nonlinear  system  is  defined  in  such  a  way  that  the  parameters  which  characterize  the 
system  can  be  evaluated  experimentally. 

The  method  is  confined  to  those  nonlinear  systems  in  which  the  remote  past  becomes 
less  and  less  relevant  to  the  behavior  of  the  system  as  we  push  it  back  in  time.  It  is 
further  restricted  to  those  systems  in  which  fast  changes  in  the  output,  relative  to  the 
input,  do  not  occur.  The  reason  for  these  restrictions,  as  will  become  apparent,  is  to 
avoid  an  excessive  number  of  terms  in  the  Laguerre  and  Hermite  expansions  respec¬ 
tively. 

As  Dr.  Wiener  has  shown,  the  proper  probe  for  the  general  investigation  of  non¬ 
linear  systems  is  the  output  of  a  shot-effect  generator  (gaussian  noise  with  a  flat  power 
spectrum).  The  shot  effect  is  the  most  general  probe  and  enables  nonlinear  systems 
to  be  classified  for  all  inputs  (1). 


Definitions 

To  simplify  the  description  of  the  method,  it  is  convenient  at  this  point  to  define 
certain  quantities  and  relations  that  are  useful  in  the  development  of  the  method. 

A.  The  nth  Laguerre  polynomial, 

L  (x)  =  —  eX  — —  (xne~X) 
nv  '  n!  ,  n  v  ’ 

dx 


B.  The  normalized  Laguerre  functions  hn(x)  are  defined 


h  (x)  =  e~X/2  L  (x) 
n'  '  nv  ' 


(1) 


The  following  orthogonality  relation  exists  for  these  functions: 


00 


h  (x)  h  (x)  dx  = 
m  '  n'  ' 


'0 


1  if  m  =  n 


0  if  m  *  n 


(2) 


C.  The  nth  Hermite  polynomial 

2 


Hn(x)  =  (-Dnex  (£)  e 


n  2 

-x 


D.  The  normalized  Hermite  polynomials  4>  (x)  are  defined 
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<j>  (x)  = 


H  (x) 
n  ' 


n 


1/2 


(2n  n!  Ml/Z) 

E.  The  normalized  Hermite  functions  are  defined 


(3) 


4 i  (x)  =  e"X  <b  (x) 
n  '  n  ' 


(4) 


These  functions  form  a  normal  orthogonal  set  over  the  interval  -«  to  «.  Conse¬ 
quently,  we  have  the  relation 


,00 


'-0C 


4>  (x)  <!)  (x)  e~X  dx  = 

mv  '  nv  ' 


1  m  =  n 


0  m  4-  n 


(5) 


Description  of  the  General  Method  of  Characterization  and  Synthesis 

Now  let  us  consider  a  nonlinear  system  whose  input  and  output  are  the  stationary 
time  series  x(t)  and  y(t),  respectively.  In  general,  y(t)  is  some  function  of  the  past  of 
x(t).  Consequently,  it  is  of  interest  to  obtain  a  set  of  parameters  which  characterize 
the  past  of  x(t).  The  past  of  x(t)  can  be  expanded  in  a  series  of  Laguerre  functions  as 
indicated  in  Eq.  6 


x<-‘>=  X  unhn(t)  (6) 

n 

The  Laguerre  functions  form  a  complete  set  over  the  interval  zero  to  infinity  and  the 

coefficients  u^  form  a  set  of  variables  which  characterize  the  past  of  x(t). 

These  Laguerre  coefficients  of  the  past  of  a  time  function  x(t)  are  easily  generated 

in  practice.  They  are  simply  the  output  of  a  Laguerre  network  (a  network  whose  impulse 

response  is  hn(t))  whose  input  is  x(t).  This  can  be  seen  as  follows.  Referring  to 

Fig.  X-5,  the  output  at  the  present  time  (t=0)  is  given  by  the  convolution  of  x(t)  with 

h  (t).  That  is, 
n'  ' 

00 

x(-z)  hn(z)  dz  (7) 

Let  us  compare  this  expression  for  y  (0)  with  the  expression  for  un  which  is  obtained 


y(o)  = 


:(z)  hn(0-z)  dz  = 


INPUT 

LAGUERRE  NETWORK 

OUTPUT 

IMPULSE  RESPONSE  IS 
h„(t) 

Fig.  X-5 

Block  diagram  of  Laguerre  network. 
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from  Eq.  6  by  making  use  of  the  orthogonality  property  of  Eq.  2.  The  latter  expres¬ 
sion  is 


00 

x(-t)  hn(t)  dt  (8) 

Comparison  of  Eqs.  7  and  8  shows  that  y  (0)  =  u^.  Thus  we  have  demonstrated  that  the 
coefficients  un  can  be  generated  by  the  method  shown  in  Fig.  X-5.  The  actual  form  of 
the  Laguerre  network  and  the  method  of  synthesis  of  this  network  will  not  be  discussed 
here.  For  a  complete  discussion  of  these  topics  the  reader  is  referred  to  reference  2. 

Since  the  set  of  coefficients  un  characterize  the  past  of  x(t),  any  quantity  dependent 
only  on  the  past  of  x(t)  may  be  represented  as  a  function  of  these  coefficients.  Thus  for 
the  nonlinear  system  with  input  x(t)  and  output  y(t)  we  can  write 


y(t)  =  F 


lV 


u 


2’ 


(9) 


For  reasons  that  will  soon  become  apparent  it  is  convenient  to  express  this  function  F 
as  a  Hermite  function  expansion  of  the  variables  u^  through  ug.  We  then  have 


y(t)  "  X  X  '  '  ‘  X  ai,  j,  •  •  • ,  h  ‘  •  -  \(Us)  e 


2  ,  2  2 

u .  +u0  +  .  .  .  +u 
12  s 


1  3 


(10) 


In  order  to  simplify  this  equation  let  $(<*)  represent  the  polynomial  ^(Uj)  .  .  •  <J>h(us),  and 

A  represent  the  corresponding  value  of  a.  .  h •  Then  Eq.  10  becomes 

a  r  ^  J<  •  ■  •  >  11 


y(t)  =  X  A°  4(a)  e 


2  ,  .2 

1  s 


(11) 


If,  in  Eq.  10,  the  indices  i,  j,  .  .  . ,  h  range  from  0  to  N  -  1,  then  a  ranges  from  0  to 


NS  -  1. 


To  characterize  the  nonlinear  system,  it  is  necessary  to  evaluate  the  coefficients 

A  inEq.  11.  To  this  end,  consider  the  time  average 
a 


2  ,  2 
u.  +  .  .  .  +u 
1  s 


(12) 


y(t)  «(P)  =  ^  ao  *(°)  *(P)  e 

a 

Let  us  take  advantage  of  the  ergodic  theorem  in  evaluating  the  right-hand  member  of 
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this  equation.  We  then  have 


y(t)  *(P)  = 


I  00 


2  2 

u .+...+  u 
1  s 


S(a)  $(p)  e 


P(u  ,  . . . , u  )  du  ...  duc 

J.  SI  5 


-00 


(13) 


As  mentioned  earlier,  the  probe  for  the  investigation  of  the  nonlinear  system  is  the  out¬ 
put  of  a  shot-effect  generator.  Thus,  for  the  evaluation  of  the  coefficients  Aa  in  Eq.  13, 
the  input  x(t)  to  the  nonlinear  system  is  a  gaussian  signal  with  a  flat  power  spectrum. 

Now  recall  that  the  parameters  u^  through  ug  which  characterize  the  past  of  x(t)  can  be 
obtained  as  the  outputs  of  the  linear  passive  Laguerre  network  whose  input  is  x(t).  Thus, 
if  x(t)  has  a  gaussian  distribution,  then  u^  through  ug  all  have  gaussian  distributions. 
In  fact,  it  can  be  shown  from  the  properties  of  the  Laguerre  network,  or  from  the  defini¬ 
tion  of  the  Laguerre  functions,  that  u.  through  u  all  have  independent  normalized 
gaussian  distributions  when  x(t)  is  gaussian  with  a  flat  power  spectrum  and  the  proper 
variance.  Thus  we  have 


P(ur 


’Us} 


(2tt)"s/2  e 


2  2 

U.+...+U 

1  s 


(14) 


so  that  Eq.  13  becomes 

y(t)  *(P)  =  (2*)'s/z  Aa 

a 


,oo 


'-00 


,  oo 


'-00 


/  2^  2v 

-(u.  +  .  .  .  +u  ) 
§(a)  $((3)  e  S  du 


.  du 


(15) 


The  reason  for  the  choice  of  Hermite  functions  for  the  expansion  of  Eq.  9  now 
becomes  apparent.  The  joint  gaussian  probability  density  (Eq.  14)  supplies  the  neces¬ 
sary  exponential  factor  in  the  Hermite  function  expansion  (Eq.  12)  to  enable  us  to  take 
advantage  of  the  orthogonality  relation  of  Eq.  5  in  evaluating  the  coefficients  A  .  Taking 
advantage  of  this  orthogonality  we  have: 

Ap  =  (2ir)s/2  y(t)  *(p)  (16) 

The  Ag's,  together  with  a  knowledge  of  the  number  of  Laguerre  coefficients  used,  serve 
to  characterize  the  nonlinear  system  under  test.  The  setup  for  the  experimental  evalua¬ 
tion  of  the  coefficients  Aq  is  shown  in  Fig.  X-6. 

Now  let  us  consider  the  problem  inverse  to  that  discussed  above.  Suppose  we  are 
given  a  set  of  coefficients  Aq  and  asked  to  find  a  nonlinear  system  corresponding  to 
these  coefficients.  Equation  1 1  is  the  guide  for  the  synthesis  problem.  This  equation 
tells  us  that,  for  each  c,  we  must  form  $(c)  and  multiply  it  by  Aq  and  the  exponential 
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Fig.  X-6 

Block  diagram  of  the  circuitry  for  the  classification  of  nonlinear  systems. 


Fig.  X-7 

Block  diagram  for  the  synthesis  of  nonlinear  systems. 

exp -(u*  +  .  .  .  +  u^j/Z.  Then  each  such  product  must  be  added.  In  practice,  the  number 
of  multipliers  is  reduced  if  we  first  form  the  sum  of  the  products  Aq  $(g)  and  then  mul¬ 
tiply  by  the  exponential  function. 

The  exponential  function,  exp  -  (u ^  +  .  .  .  +  u^/2,  can  be  obtained  as  the  product  of 

s  exponential  function  generators  whose  inputs  are  respectively  u.  through  u  .  Such 

2  is 

generators  give  an  output  of  exp  -(x  /2)  when  the  input  is  x.  They  are  realizable,  among 
other  ways,  in  the  form  of  a  small  cathode -ray  tube  with  a  special  target  to  generate 
the  exp  -(x^/2)  function. 

The  block  diagram  of  the  circuitry  used  for  the  synthesis  procedure  is  shown  in 
Fig.  X-7. 

Example 

In  order  to  fix  ideas,  let  us  consider  a  simple  example  which  is  particularly  suited 
to  characterization  and  synthesis  by  Dr.  Wiener's  method.  Let  the  given  nonlinear 
system  contain  no  storage  elements.  Further  let  its  input -output  characteristic  be  given 
by  the  equation 

y(t)  =  e'x2(t)/2  (17) 

where  x(t)  is  the  input  and  y(t)  is  the  output. 

In  both  the  characterization  and  synthesis  procedures  described,  the  Laguerre 
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network  is  used  to  introduce  dependence  of  the  output  upon  the  past  of  the  input.  The 
nonlinearity  is  brought  about  by  the  Hermite  polynomial  generator.  For  the  simple 
example  under  consideration,  there  is  no  dependence  upon  the  past  and  thus  we  can  by¬ 
pass  the  Laguerre  network.  In  the  experimental  procedure,  the  fact  that  this  nonlinear 
system  has  no  storage  can  easily  be  determined  from  the  results  of  a  priori  tests  made 
on  the  system. 

To  proceed  with  the  characterization  of  the  simple  example,  we  notice  that  as  a 
result  of  by-passing  the  Laguerre  network  (Fig.  X-6)  the  variables  u^  through  ug  are 
replaced  by  the  single  variable  x(t).  Equation  10  then  becomes: 


y(t)  =  ^  ai  ^(x)  e"x 


Equation  16  becomes 

ai  =  (2tt) */2  y(t)  ^(x) 


(18) 


(19) 


Let  us  make  use  of  the  ergodic  theorem  to  evaluate  this  time  average  as  an  ensemble 
average.  Using  Eq.  17  we  can  write 


ai  =  (2tt) 


1/2 


cj^(x)  e-X  /2  p(x)  dx 


(20) 


But  since,  in  the  test  setup,  x(t)  is  the  output  of  an  ideal  shot-effect  generator,  the 
probability  density  of  x  is 

p(x)  =  T/2  e"X  ^  (21) 

(2*)1/2 

Thus 


a.  = 
l 


2 

cJk(x)  e~X  dx 


(22) 


Referring  to  Eq.  3  and  the  definition  of  the  Hermite  polynomial,  it  is  seen  that 

<|>o(x)  =  TT-1/4  (23) 

With  this  result,  Eq.  22  becomes 


a.  =  tt 
l 


1/4 


,  00 


4>t(x)  4>o(x)  e"x  dx 


—  00 


(24) 
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As  a  consequ 
1 

a  =  it 
o 

a.  =  0 
1 

The  relations  of  Eq.  25  together  with  the  knowledge  that  the  Laguerre  network  has 
been  by-passed  serve  to  characterize  the  given  nonlinear  system. 

Now  let  us  consider  the  inverse  problem  with  respect  to  the  same  example.  Suppose 
we  are  given  the  coefficients  of  Eq.  25  along  with  the  associated  information  in  the 
preceding  paragraph.  Our  job  is  to  synthesize  the  corresponding  nonlinear  system. 

The  guide  for  the  synthesis  is  Eq.  18  which  corresponds  to  Eq.  1 1  in  the  more  com¬ 
plicated  case.  Following  the  general  procedure  for  the  synthesis  described  earlier  it 
is  readily  seen  that  the  resulting  system  is  that  shown  in  Fig.  X-8.  Figure  X-8(a) 
shows  how  the  synthesis  takes  form  according  to  the  formal  dictates  of  Eq.  18.  The 
simplified  equivalent  system  obtained  from  Fig.  X-8(a)  by  inspection  is  shown  in 
Fig.  X-8(b).  We  see  that  for  the  simple  example  considered  the  synthesized  system 
consists  solely  of  the  "function  generator,"  a  component  which  in  the  more  complicated 
case  will  form  only  a  part  of  the  synthesized  system. 

At  this  point  it  is  of  interest  to  note  that  we  have  at  our  disposal  the  two  scale  factors 
of  the  Laguerre  and  Hermite  functions.  For  simplicity  in  the  presentation  of  the  basic 
concepts  of  Dr.  Wiener's  method,  these  scale  factors  have  been  chosen  to  be  unity.  In 
practice,  they  can  be  used  to  reduce  the  number  of  parameters,  A  ,  necessary  to  char¬ 
acterize  the  system.  They  are  to  be  judiciously  chosen  on  the  basis  of  a  priori  knowl¬ 
edge  of,  or  tests  made  on,  the  nonlinear  system. 


ence  of  the  orthogonality  relation  of  Eq.  5  we  have  the  result 

/4 
i  *  0 


(25) 


(a) 


x(t) 


FUNCTION 

y  (t)  =  e  8 

GENERATOR 

x» 

e  8 

(b) 


Fig.  X-8 

2 

(a)  Synthesis  of  the  nonlinear  characteristic  y  =  exp  -(x  /2) 
according  to  the  general  procedure  of  Fig.  X-7.  (b)  Equivalent 
reduced  form  of  part  (a). 
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In  concluding  this  discussion  of  Dr.  Wiener's  method  of  nonlinear  systems  the  ques¬ 
tion  naturally  arises  as  to  the  number  of  parameters,  A  ,  that  may  be  necessary  to 
characterize  a  practical  system.  In  the  simple  example  worked  earlier,  only  one  param¬ 
eter  was  necessary.  However,  suppose  we  have  a  system  which  requires  s  Laguerre 
coefficients  to  sufficiently  characterize  that  portion  of  the  past  of  x(t)  which  influences 
the  behavior  of  the  system.  Further  suppose  we  decide  that  a  sufficient  approximation 
to  the  performance  of  the  system  is  obtained  by  letting  the  indices  i,  j,  .  .  .,h  in  Eq.  10 

g 

range  from  0  to  N  -  1.  Then  there  are  N  parameters  (Aq)  necessary  to  characterize 
the  system.  Without  a  doubt  this  number  can  become  quite  large  in  many  practical 
cases.  However,  with  the  freedom  that  exists  in  nonlinear  systems,  one  can  hardly 
expect  to  apply  such  a  general  approach  without  considerable  effort. 

A.  G.  Bose 
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PROPERTIES  OF  SECOND -ORDER  AUTOCORRELATION  FUNCTIONS 


In  the  study  of  second-order  autocorrelation  functions  the  following  two  definitions 
are  used: 


4>1i1(Ti»  Tz}  =  fi(t)  fi(t  +  Ti>  fi(t  +  Ti  +  Tz) 


<(>lll(Tl.  T2)  =  f  2(t)  f  x(t  +  Tj)  f  x(t  +  T2) 


(1) 


(2) 


where  f  ^(t)  is  a  stationary  time  series. 

By  shifting  the  time  origin  in  Eq.  1  we  obtain  the  following  relations  between  the  two 


definitions: 


+llltTl'  t2>  "WvV 

hiFV’b  =<t>in(TrT2> 


(a) 


(b) 


<f,Ill(Tl*  "t2>  =  “t2)  (c) 


hll(-Tl'-T2>  =  +lll(-T2’  Tl) 


(d) 


(3) 


From  the  definition  of  ^(t  ,  t^)  it  is  clear  that 
4,m<Ti’T2)  =  ,»,Tii(t2-  ti> 

Making  use  of  this  relation,  we  obtain  from  3(a)  and  3(d)  the  relations 

‘hll^l’  T 2)  =  t2^ 

+  111(_t1’~t2>  =  4>111(VT2) 


(4) 


(3a1) 


(3d1) 


The  relations  between  t^)  an<3  <t>^i(Ti’  rz ^  as  £iyen  by  Eqs.  3b,  c,  a',  and  d'  are 

illustrated  in  Fig.  X-9.  From  these  relations  it  can  be  shown  that  the  <j>^^(T^,  T?)  plane 
is  the  4b  11  (Ti’  Tz)  Plane  rotated  about  the  axis  by  180°.  That  is,  the  4>111(T1>  plane 
is  the  mirror  image  of  the  4)111(T1>  t^)  plane  about  the  t7  axis.  Consequently,  given 
the  properties  of  4)111(T1,  T2)«  we  can  easily  establish  the  corresponding  properties  of 
4>1 1 1  ( t  1 ,  r  )  and  vice  versa. 

By  a  process  similar  to  that  used  to  derive  Eq.  3,  we  obtain  from  Eq.  1  the  relation 

t2^  =  ^  1 1 1  ^ “ T2 *  “Ti^  3  ^11  i^"Tr  Ti  +  t2^  3  ^  1 1 1  ^ T2 ’  _Ti  “  Tz)  ^ 

This  shows  that  a  knowledge  of  one  quadrant  of  4)111(T1»  T2)  is  sufficient  to  determine  the 
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4  2 

\^TZ 

(ri  >T2'  =  " 

\ 

0 

r 

V(-T 

I’  2 7 


^,n(r|,r2)- PLANE 


remaining  three  quadrants,  although 
such  a  procedure  as  indicated  by  this 
equation  would  not,  in  general,  be  the 
most  practical  method  for  constructing 
the  remaining  quadrants. 

Now  let  us  consider  symmetry 
properties  of  second-order  correla¬ 
tion  functions.  From  Eq.  4  we  see 
that  ^  ,  t^)  is  always  symmetri¬ 

cal  about  the  line  =  r^.  It  follows 
from  this  fact,  and  the  mirror  image 
relationship  between  the 
and  <j>111(T1,T2)  planes,  that  4>111(T1>T2) 
is  always  symmetrical  about  the  line 

T1  =  -t2‘ 

It  is  of  interest  to  show  that  if  f(t) 
is  a  stationary  time  series  such  that 
the  statistical  behavior  of  f(t)  is  iden¬ 
tical  with  the  statistical  behavior  of 
the  time  series  f(-t),  then  T 2) 

is  symmetrical  about  the  line 
in  addition  to  being  symmetrical  about 
the  line  =  ~Tz'  Symmetry  about  the 

line  t ^  means  that  T ^ =  ^lll^T2’  Tl^'  Sufficient  conditions  for  this  equality 

can  be  established  by  expressing  T2^  as  an  ensemble  average. 


Fig.  X-9 

Relation  between  4>^  ^  -t^  ,  t^)  and  ,  t^)  . 

$ 

The  4»211  ('r1»  'r2)  Plane  is  the  <l>iii(Ti»  T  z) 
plane  rotated  180°  about  the  axis. 


^lll^l’  rz) 


xyz  P(x,  y  ,  z  +  j  dx  dy  dz 


(6) 


In  this  expression  P  (x,  y  ,z  ,  ]  is  the  third  probability  density  of  f(t).  The  sub- 

\  1  T1  T2/ 

scripts  of  the  variables  denote  the  time  of  observation  relative  to  the  time  t  =  0  when 
the  variable  with  no  subscript  (x  in  this  particular  expression)  is  observed.  Using  this 
notation  we  can  write 


^111^T2’  t1>  - 


x’y’z'  P^x'.y^.  z^+T^  dx'  dy'  dz' 


(7) 


Now,  in  Eq.  7  let  us  change  the  dummy  variables  of  integration.  Let  x1  =z,  y'  =  y, 
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and  z'  =  x.  Then  Eq.  7  becomes: 


xyz  P  ( x 


vv 


z  )  dx  dy  dz 


(8) 


To  put  this  equation  into  a  form  convenient  for  comparison  to  Eq.  6,  let  us  translate 
the  t  axis  to  establish  the  reference  axis  (t=0)  at  the  time  when  the  variable  x  is 
observed.  The  result  is: 


hn( 


xyz  P(x,y_T 


"Tl"T2, 


dx  dy  dz 


(9) 


Comparing  Eqs.  6  and  9  we  see  that  a  sufficient  condition  for  the  existence  of  sym¬ 
metry  about  the  line  is 


P 


z 


VT 


=  P 


(10) 


This  relation  clearly  holds  for  all  stationary  time  series  f(t)  which  have  the  same 
statistical  behavior  as  the  corresponding  time  series  f(-t).  If  f(t)  is  a  periodic  function, 
Eq.  10  will  hold,  and  thus  there  will  be  symmetry  about  the  line  t  =  t7  for  all  functions 
having  a  time  origin  about  which  the  function  is  symmetrical. 

It  has  been  shown  that  t^)  a^ways  symmetrical  about  the  line  and 

is  sometimes  symmetrical  about  the  line  =  t^.  The  question  naturally  arises  as  to 
the  existence  of  symmetry  about  the  and  axes.  It  is  next  shown  that,  with  the 
exception  of  functions  which  give  a  constant  second-order  autocorrelation  function, 
symmetry  about  the  and  axes  is  not  present.  This  result  will  be  proved  by  the 
method  of  contradiction. 

By  performing  a  translation  of  the  time  origin  in  Eqs.  1  and  2,  it  can  be  shown  that 
if  symmetry  exists  about  either  the  t  or  axis,  then  it  exists  about  both  axes.  Hence, 
for  the  proof  it  is  sufficient  to  show  that  the  assumption  of  symmetry  about  the  axis 
leads  to  a  contradiction. 

The  assumption  is 

4,iii(TrT2)  =  4,in('TrT2)  (11) 


Using  this  relation  and  Eq.  3a1  we  obtain  the  relation 


^111^  Tr  T2^  ^111^  Tl’  T2^ 


(12) 
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Expressing  -t^  ,  t^)  in  terms  of  definition  2, 


<l>lll('TrT2)  =  +  lll('TllT2  ■  Tl> 


(13) 


Hence,  from  Eqs.  12  and  13 


♦in 


< -T1  ’  t2>  =  + 


* 

111 


(14) 


If  Tj  is  held  fixed  in  Eq.  14,  (j^  (-t  .t^)  Peri°dic  with  respect  to  t ^  with  a  period 
equal  to  the  value  of  t^.  This  result  is  a  contradiction  when  f^(t)  is  a  random  stationary 
time  series,  for  the  second-order  autocorrelation  function  of  a  random  function  is  not 
periodic  with  respect  to  either  or  t^.  The  result  is  also  a  contradiction  when  f^(t) 
is  a  periodic  function,  for  the  second-order  autocorrelation  of  a  periodic  function  is 
periodic  with  respect  to  for  a  given  but  the  period  is  independent  of  Hence 

assumption  11  leads  to  a  contradiction.  This  proves  that  T2)  i-s  n°t  symmetri¬ 

cal  with  respect  to  the  or  axes  except  for  the  degenerate  case  in  which 
is  a  constant. 

Y.  W.  Lee,  A.  G.  Bose,  J.  Y.  Hayase 


D.  FIELD  MAPPING  BY  CROSSCORRELATION 


The  problem  of  determining  the  location  of  n  sources  of  independent  random 
stationary  time  series  by  the  method  of  crosscorrelation  has  been  considered  under 
restricted  conditions. 

The  sources  lying  on  a  plane  are  designated  by  S.,  where  j  =  1,  .  .  .  ,  n,  as  shown  in 

J 

Fig.  X-10.  On  the  same  plane,  three  points,  L^,  i  =  1,  2,  3,  forming  an  equilateral  tri¬ 
angle,  are  chosen  as  listening  posts.  The  separation  between  these  points  is  d.  The 
signals  received  at  L^  are 

n 

gi(t)  =  ^  fj(t  "  kdij)  (!) 

3=1 

where  f.(t)  ^  0  is  the  signal  generated  by  S.,  d..  is  the  distance  between  the  points  S.  and 
J  J  J 

L.,  and  k  is  the  reciprocal  of  the  velocity  of  light.  The  second-order  crosscorrelation 

function  of  g^(t)  is 


The  definition  of  the  crosscorrelation  function  is  similar  to  the  definition  of  the 
autocorrelation  function,  <t>m(Tj,  r2)'  discussed  in  subsection  C. 
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Fig.  X-10 

Distribution  of  sources  and  the 
location  of  listening  posts. 


Fig.  X-ll 
the  areas  ii 

lie  and  the  areas  in  which  points  (t^  ,  x^^lie. 


Relation  between  the  areas  in  which  sources  S. 

.(i) 


<l>123(Tr  T2>  =  gl(*)  §2^  +  Tl>  §3(t  +  T2> 


n  n  n 


=  Z  X  X  fj(t  '  kdlj)  fi(t  +  T1  •  kd2i>  fp<‘  +  T2  -  kd3p> 
j=li=l  P=1 

As  a  consequence  of  the  independence  of  the  sources,  Eq.  2  is  equivalent  to 

n 

<t’l23(Tr  T2>  =  X  fj(t  “  kdlj)  fj<t  +  Tl  "  k d2j>  fj(t  +  T2  -  k d3j) 
j=l 

n  3  3 


(2) 


+  )  )  /  a.  f.(t+T  -kd...)  f.(t  +  T  -kd  .) 

Jip  V  1-1  ly  Y  P-1  PT 


3=1  i=l  P=1 
i*P 

+  constant  terms  (n  >  1) 

* 


(2a) 


where  =  0.  The  value  of  <1)123(ti»  T2^ 

T<ll)  =  k<d2i  '  V 


T2*  =  k(d3i  -  dli> 


(3) 


IS 
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=  f3  (t  -  kd. .)  +  A  f2(t  -  kd.  .) 
1  '  lr  iv  lr 


n  3  3 

+  III 

1=1  1=1  D=1 


JiP  J 


f.  t+  T 


(i) 

i-1 


f. 

J 


t  +  T 


(i) 

P-1 


+  other  constant  terms  (n  >  1) 


In  the  simple  case  where  n  =  1, 

^123  {k(d21  "  dll^’  k^d31  ~  dll^}  =  M*  ”  kdll^ 


Since  it  can  be  shown  that 


(4) 

(5) 


yt)  ^  fj(t)  f  j(  t  +  tx)  L(t  +  t2) 
the  time  displacements  and  t^, 


(6) 


which  are  related  to  the  distances  d^,  d^,  and 


31 


by 


"l  ^  -  k^d2 1  dll^’  T2  ^  ”  k^d3 1  dll^ 


(7) 


can  be  determined  by  plotting  T 2)  vs'  Tl’  T2  the  aid  of  a  correlator. 

When  n  >  1,  the  summation  terms  in  Eq.  4  may  introduce  an  error  in  the  measured 


.(i) 


.(i) 


'  and  r2‘  K  can  be  shown  that  this  error  depends  on  the  relative  position  of  the  points 

S.,  and  the  nature  of  the  random  function  f .( t) .  If  the  summation  terms  introduce 

J  (i)  J  (i) 

negligible  error  on  the  measured  '  and  ’ ,  then  those  determine  the  location  of  S^. 

The  difference  in  the  distances  d?.  -  d...  d,.  -  d..,  and  d9.  -  d~.  can  be  expressed 

\  1  .\  ul  II  ol  il  c.1  jl 


in  terms  of  r\x/  and  t\v  by  the  equations 


1 


.(i) 


2i 


d, . 

3i 


2i 


li 


dli 


d_ . 
3i 


.(i) 


T(i)  -  T(i) 
1  2 


(8) 


J 
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Three  families  of  hyperbolas,  whose  respective  foci  are  L^  and  L^,  L^  and  L^,  and  L^ 
and  L^,  are  available  for  determining  the  location  of  S..  Any  two  families  of  hyperbolas 
are  sufficient  to  determine  the  location  of  S^  which  is  given  by  the  intersection  of  a 
branch  from  each  family.  The  choice  of  the  two  families  is  dictated  by  the  geometrical 
configuration  of  S^  and  the  fixed  points,  L^ ,  L^,  and  L^.  The  differences  in  the  dis¬ 
tances  from  the  two  foci  to  S^  for  a  given  family  are  given  by  Eq.  8.  This  hyperbolic 
grid  system  is  the  same  as  that  of  Loran. 

With  the  restriction  that 


and 


lTli}|  <  kd 
It2^  I  ^  kd^ 


(9) 


the  relation  between  the  areas  in  which  the  sources,  S^,  and  the  points,  lie 

can  be  shown  as  in  Fig.  X-ll.  The  correspondence  between  the  areas  in  the  xy  plane 
and  the  plane  is  shown  by  primed  integers  and  unprimed  integers. 

J.  Y.  Hayase 


E.  CORRELATION  DETECTION  OF  A  PERIODIC  SIGNAL  MULTIPLIED  BY  NOISE 


If  a  periodic  signal  S'(t)  is  multiplied  by  statistically  stationary  random  noise,  the 
product  may  be  written 

S'(t)  N'(t)  =  [s(t)  +  S]  [N(t)  +  N]  (1) 

where  S  and  N  are  the  mean  values  of  S'(t)  and  N'(t),  respectively.  If  N  is  not  zero, 
the  periodic  signal  can  be  detected  by  autocorrelation,  but  if  it  is  zero,  detection  by 
autocorrelation  alone  is  not  possible. 

By  filtering  in  the  frequency  domain,  Eq.  1  may  be  reduced  to  S(t)  N(t).  For  this 
case,  the  autocorrelation  function  of  the  product  is  the  product  of  the  autocorrelation 
functions 

4>sn(t)  =  ^S^  <Mt) 


Although  4*g(r)  is  periodic,  <j>^(r)  tends  to  zero  for  t  large,  so  that,  in  general,  the  peri¬ 
odic  signal  may  not  be  detected  by  autocorrelation.  If  the  factors  of  Eq.  2  are  separated, 
or,  if  4)]vj(T)  Is  made  to  approach  a  constant  other  than  zero,  detection  by  autocorrelation 
is  possible. 

The  factors  may  be  separated  by  taking  the  logarithm,  thus 

log  ]SN |  =  log  is |  +  log  |N |  (3) 
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<f>(T)  =  ^log  S(t)  +  ^log  N(t)  +  2  log  log  (4) 

The  periodic  signal  is  revealed  by  the  first  term  of  Eq.  4;  the  other  terms  approach 
constants  for  r  large. 

The  noise  factor  of  Eq.  2  may  be  made  nonzero  for  t  large  by  a  nonlinear  operation 
on  the  product  S(t)  N(t).  If  the  product  S(t)  N(t)  is  applied  as  an  input  to  a  nonlinear 
device  the  transfer  characteristic  of  which  may  be  written 

00 

«  =  (5) 

0 

then  the  output  autocorrelation  function  may  be  written 


00  00 

<|>(t)  =  ^  An  [s(t)  N(t)]n  ^  A.  [s(t+T)  N(t+r) 
n=0  j=0 


j 


=  [S(t)  S(t+r)]n  [N(t)  N(t+r)]n 

n=0 


00  00  _ 

+  An  A  S(t)n  S(t+T)j  N(t)n  N(t+T)j  (6) 

n=0  j  =  0 

n^j 


Some  noise  terms  having  nonzero  means  are  multiplied  by  signal  terms.  Since  it  is  not 
possible  for  a  random  function  raised  to  an  even  power  to  have  a  zero  mean,  any  non¬ 
linear  device  having  even-power  terms  in  the  power  series  expansion  of  its  transfer 
characteristic  permits  the  detection  of  the  periodic  signal.  If  the  probability  density 
function  of  N(t)  has  even  symmetry  about  the  origin,  then  the  odd-power  terms  of  Eq.  5 
do  not  aid  in  the  detection  of  the  periodic  signal;  but  if  not,  the  odd-power  terms  may 
aid. 

A  commercially  available  logarithmic  device  and  a  full-wave  linear  rectifier  provided 
the  more  interesting  transfer  characteristics  for  the  experimental  work.  The  character¬ 
istics  are 

log  k^x  a  ^  x  <  b 

k^x 


Limiter 

(odd) 


y  =  S 


-a  <  x  <  a 


(?) 


3 


71 


Fig.  X- 14  Fig.  X-  1 5 

(}>(t)  for  full-wave  log  SN  and  full-wave  SN.  <j>(T)  for  full-wave  SN. 
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Fig.  X- 1 6 

Oscilloscope  photographs:  time  scale,  1  msec;  sine  wave,  5.  1  kc/sec; 

noise,  4  kc/sec  centered  at  21.9  kc/sec. 
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Full-wave 

(even) 


r 


K  lQg 


k2x 


a  ^  |  x  |  <  b 
0  <  |  x  |  ^  a 


(8) 


Full-wave  y  =  |kAx|  (9) 

(even) 

If  the  input  of  an  ideal  logarithmic  device  is  E  sin  co^t,  then  the  output  is 
y  =  B  log^  |E  sin  co^t|  ,  and  the  autocorrelation  function  of  the  output  is 


4>(t)  =  (B  loge  E)2  + 


(10) 


where  B  is  determined  by  the  characteristics  of  the  log  device  and  the  gain  of  the  output 
amplifiers.  The  autocorrelation  function  for  the  output  of  a  full-wave  linear  rectifier, 
y  =  |  E  sin  go  1 1  is 


♦w  -4  [(! - 


^1T 


cos  co  t  +  sin  go^t 


(11) 


Calculated  and  experimental  results  are  shown  in  Figs.  X-12  and  X-13  for  these  two 
cases. 

Figure  X-14  shows  the  experimentally  obtained  autocorrelation  function  for  full-wave 
log  SN  and  full-wave  SN  for  identical  correlator  control  settings.  S  is  a  sine  wave  and 
N  is  band-limited  gaussian  noise.  Figure  X-15  shows  the  autocorrelation  function  for 
full-wave  SN  for  increased  correlator  sensitivity  (and  increased  error).  The  similarity 
of  the  curves  in  Fig.  X-12  and  X-14  suggests  that  the  factors  did  separate  as  indicated 
by  Eq.  3.  The  autocorrelation  function  for  limiter-log  SN  does  not  reveal  the  periodic 
component. 

Oscilloscope  photographs  such  as  those  in  Fig.  X-16  proved  to  be  a  valuable  aid. 

C.  E.  McGinnis 
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F.  AN  IMPROVEMENT  IN  THE  DIGITAL  ELECTRONIC  CORRELATOR 

The  correlator  employs  a  compensating  circuit  to  minimize  drift  error.  When 
processing  low-frequency  information  (dc  to  500  cps),  this  compensating  circuit  acts 
upon  a  reference  voltage  of  zero  volts  which  is  applied  to  the  inputs  of  the  correlator  in 
intervals  alternate  to  those  during  which  the  signals  to  be  correlated  are  applied.  The 
intervals  are  such  that  samples  taken  by  the  machine  are  alternately  signal  and  refer¬ 
ence  samples.  Difficulty  in  establishing  a  constant  reference  voltage  of  zero  volts  during 
the  reference  sample  intervals  has  led  to  the  development  of  a  very  satisfactory  circuit 
for  this  purpose. 

The  circuit  previously  in  use  consisted  of  two  tubes  in  series  working  between 
+  150  volts  and  -150  volts.  These  tubes  were  alternately  gated  on  and  off,  and  their 
junction  was  connected  to  the  input  tube  of  the  correlator.  A  potentiometer  in  series 
with  the  two  tubes  adjusted  the  voltage  at  the  junction,  during  the  on  time,  to  zero. 
During  the  off  time,  the  grid  of  the  correlator  input  tube  was  at  high  impedance  and  free 
to  receive  the  signal  fed  to  it  through  a  series  resistor.  This  method  suffered  from  the 
obvious  difficulty  of  being  very  dependent  upon  the  tube  characteristics  and  upon  the 
stability  of  the  plus  and  minus  B  supplies. 

The  improved  circuit  is  shown  in  Fig.  X-17.  During  the  intervals  for  which  both 
sections  of  the  tube  are  on,  there  is  a  very  small  resistance  (about  150  ohms)  between 
points  a  and  b  and  the  signal,  f(t),  passes  freely  to  the  grid  of  the  correlator  input  tube. 
During  the  intervals  for  which  both  sections  of  the  6AS7  are  off,  the  grid  of  the  corre¬ 
lator  tube  is  isolated  from  the  source  f(t),  and  it  returns  to  zero.  This  circuit  has  the 
advantage  of  being  almost  completely  independent  of  changes  in  the  6AS7  characteristics. 


f(t) 


Fig.  X-17 

Improved  circuit  for  the  digital  correlator.  The  time  function  f(t) 
is  the  function  to  be  correlated  (10  volts  maximum  either  side  of 
zero).  It  is  to  be  obtained  from  a  low  impedance  dc  coupled  source. 
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Since  there  is  no  B  voltage  applied  to  the  6AS7,  and  since  the  filament  voltage  is  reduced 
to  minimize  hum,  the  tube  life  is  long.  This  circuit,  in  conjunction  with  the  compensa¬ 
tor,  enables  the  machine  to  process  low-frequency  information  with  an  accuracy  of  the 
same  order  of  magnitude  as  that  obtained  when  the  machine  is  ac  coupled  and  processing 
higher  frequency  information. 

A.  G.  Bose 


G.  MEASUREMENT  OF  INDUSTRIAL  PROCESS  BEHAVIOR 

Experimental  determinations  of  the  transfer  function  H(gj)  or  impulse  response  h(t) 
of  industrial  processes  are  complicated  by  the  presence  of  noise  and  by  the  difficulty 
of  applying  suitable  sinusoidal  or  impulsive  driving  functions  to  massive  industrial 
systems.  The  statistical  method  of  determining  system  characteristics  overcomes  these 
difficulties  (1).  This  method  uses  as  the  driving  function  the  random  fluctuations  already 
existing  in  the  input  variable.  The  result  is  not  usually  affected  by  random  fluctuations 
originating  within  the  process. 

A  heat  transfer  process  frequently  encountered  in  industry  has  been  chosen  for 
experimental  study.  In  this  process,  a  liquid  flowing  through  a  copper  tube  is  heated 
by  water  flowing  through  a  concentric  copper  shell.  Thermocouples  are  used  to  meas¬ 
ure  the  temperatures  of  the  liquids.  When  cold  water  from  the  building  mains  flows 
through  the  tube  and  a  mixture  of  hot  and  cold  water  from  the  same  mains  flows  through 
the  shell,  the  shell  and  tube  temperatures  fluctuate  as  shown  in  Fig.  X-18. 


Fig.  X-18 

Temperature  fluctuations  in  shell  and  tube  of  heat  exchanger. 
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Preliminary  results  obtained  with  the  low-frequency  multisignal  correlator  (2), 
while  encouraging,  indicate  that  practical  application  of  the  method  must  await  further 
refinement  of  the  measuring  system. 

S.  G.  Margolis 
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Prof.  M.  Halle 
G.  W.  Hughes 
J.  -P.  A.  Radley 


A.  RECOGNITION  OF  SILENCES 

The  class  of  speech  sounds  known  as  stops  or  plosives,  e.  g.  ,  /p/,  /k/,  /t/,  /d/, 
/b/,  and  so  on,  are  produced  by  a  complete  closure  of  the  vocal  tract  followed  by  a 
sudden  opening.  During  this  movement  the  vocal  cords  may  be  vibrating.  This  vibra¬ 
tion  distinguishes  voiced  stops  like  /b/  and  /d/  from  voiceless  ones  like  /p/  and  /t/. 

The  acoustic  correlate  of  this  articulatory  process  is  an  absence  of  energy  in  a 
broad  frequency  band  followed  by  a  rapid  onset  of  energy  in  all  frequencies.  For  the 
voiceless  stops,  the  absence  of  energy  is  more  or  less  total;  for  the  voiced  stops, 
strong  components  are  usually  present  in  a  region  below  300  cps  but  above  this  region 
there  is  hardly  any  energy.  This  relative  silence  is  thus  one  of  the  cues  that  permit  us 
to  identify  stops. 

The  operation  to  be  performed  is  that  of  determining  when  the  input  to  a  device  falls 
below  a  certain  low  level  and  remains  there  for  some  minimum  time.  To  recognize  a 
critical  level,  a  Schmitt  amplitude  discriminator  is  used  after  the  speech  signal, 
suitably  filtered  and  amplified,  has  been  rectified  (full-wave). 

In  the  first  method  tried,  the  output  of  the  Schmitt  circuit  was  given  a  certain 
amplitude  which  controlled  the  onset  and  timing  of  a  highly  unbalanced  multivibrator. 
Impulses  from  this  stage  were  fed  to  one  side  of  a  flip-flop.  When  the  silence  was 
ended,  the  flip-flop  was  returned  to  its  first  state  and  the  multivibrator  turned  off.  If 
a  silence  did  not  last  long  enough,  the  multivibrator  was  turned  off  before  emitting  its 
first  impulse.  In  fact,  at  every  half -pitch  period,  the  signal  could  be  expected  to  reach 
zero  for  a  short  time.  It  was,  however,  not  possible  to  reset  the  multivibrator  quickly 
enough,  which  led  to  cumulative  timing  errors.  A  new  approach  was,  therefore, 
developed. 

In  the  second  method,  the  output  of  the  Schmitt  stage,  a  rectangular  wave,  is  used 
to  generate  a  negative -going  triangle  whose  slope  can  be  varied  and  which  resets  very 
quickly  (with  respect  to  audio  frequencies)  to  its  origin  when  the  Schmitt  tells  it  to.  This 
triangle  is  coupled  to  a  "simplified"  Multiar  (1).  A  pentode,  normally  heavily  con¬ 
ducting,  is  cut  off  by  the  end  of  the  triangle  excursion.  In  the  Multiar  a  transformer  is 
used  in  a  feedback  path  for  rapid  cutoff  action,  but  this  is  unnecessary  in  this  applica¬ 
tion,  since  the  tube  is  cut  on  rapidly  enough  by  the  essentially  vertical  return  of  the 
triangle.  It  is  this  rapid  return  that  marks  the  end  of  a  silence.  The  location  of  the 
beginning  of  the  silence  does  not  have  to  be  known  with  great  accuracy  for  present 
purposes. 

The  circuit  was  tested  with  short  sentences  read  by  four  English  speakers  (one 
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Fig.  XI- 1 

Lower  trace  signals  stops  by  "end  of  silence"  mark.  Note  that  /n/  is  wrongly  identified. 

female  and  three  male).  Because  of  the  relatively  gentle  slope  of  the  filter  available, 
the  signal  was  highpassed  at  1500  cps  or  2000  cps  to  eliminate  the  voicing  component. 
The  critical  timing  was  set  at  about  16  msec.  This  is  the  minimum  time  for  which  the 
amplitude  must  remain  below  the  critical  level  before  a  silence  can  be  noted.  In  other 
words,  no  silence  shorter  than  the  critical  timing  is  registered.  The  hysteresis  effects 
of  the  Schmitt  circuit  were  arranged  to  be  least  for  sine  waves  between  30  cps  and 
300  cps.  A  filmed  record  was  made  with  a  dual  beam  cathode-ray  oscilloscope  showing 
a  rectified  trace  of  the  speech  on  one  trace  and  the  output  signal  of  the  "silence  recog¬ 
nizer"  on  the  other.  See  Fig.  XI- 1. 

As  expected,  the  two  criteria  (level  and  time)  used  were  necessary  but  not  sufficient 
for  the  identification  of  stops.  Different  settings  will  mark  silences  before,  or  during, 
different  classes  of  sounds.  Almost  all  stops  in  the  samples  were  marked  by  the  end 
of  a  silence.  Exceptions  to  this  were  a  few  very  weak  stops  that  did  not  reach  above  the 
critical  amplitude  (set  just  above  the  noise  region)  and  the  flapped  / 1/  as  in  "water," 
in  which  the  silence  is  extremely  short. 

On  the  other  hand,  because  of  the  very  high  cutoff  frequency  of  the  filter,  sounds 
such  as  the  nasals  and  the  semi-vowels  in  which  most  of  the  energy  is  located  in  the 
low  frequencies  were  judged  to  be  silences.  Better  filtering  might  eliminate  this. 

J.  -P.  A.  Radley 
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B.  THE  NASAL  CONSONANTS 

The  nasal  consonants  have  in  common  with  the  vowels  a  spectrum  that  is  dominated 
by  a  small  number  of  resonance  peaks,  the  so-called  formants.  In  both  the  nasals  and 
the  vowels  the  lowest  resonances  (the  first  formant)  contribute  most  of  the  energy. 
The  vowels  differ  from  the  nasals  in  having  higher  first  formants.  Furthermore,  if  the 
spectrum  contains  formants  above  the  first  formant,  these  are  much  less  attenuated  in 
the  vowels  than  in  the  nasals.  The  first  formant  of  the  nasals  is  generally  below  300  cps 
The  lowest  first  formant  in  vowels  is  found  in  [i]  and  [u],  where  it  is  located  at  about 
300  cps.  The  first  formant  in  the  other  vowels  is  well  above  300  cps.  These  facts 
suggested  that  a  comparison  of  the  energy  below  300  cps  with  that  above  might  yield 
information  useful  in  distinguishing  between  vowels  and  nasals. 

Our  speech  samples  were  isolated  syllables  consisting  of  a  vowel  preceded  or 
followed  by  a  nasal.  These  syllables  had  been  recorded  on  tape  by  one  female  and  two 
male  speakers.  Portions  50  msec  long  were  gated  out  of  the  nasal  and  the  vowel.  The 
gated  portion  was  passed  through  a  Spencer  Kennedy  302  variable  electronic  filter  in 
which  both  sections  were  set  either  to  300  cps  highpass  or  300  cps  lowpass.  This  gave 
us  an  attenuation  of  36  db  per  octave.  We  found  that  a  lesser  attenuation  (18  db  per 
octave)  yielded  no  results  of  interest.  The  output  of  the  filter  was  passed  through  a 
variable  precision  attenuator,  integrated  and  measured. 

The  following  results  were  obtained:  (The  numbers  in  the  table  are  the  extreme 
values  found  when  the  lowpass  output  was  subtracted  from  the  highpass  output.  ) 

vowels  (113  measured)  -8  db  to  +15  db 

[i]  -8  db  to  +4 . 5  db 

other  vowels  -2  db  to  +15  db 

[m]  (37)  -15.  5  db  to  -8  db 

[n]  (37)  -17.  5  db  to  -3  db 

If  we  omit  the  vowel  [i]  from  consideration  it  is  possible  to  separate  the  vowels  and 
the  nasals  quite  well  by  this  criterion:  the  vowels  give  positive  and  moderately  negative 
values,  while  the  nasals  give  strongly  negative  values.  The  vowel  [i]  could  be  made  to 
give  more  positive  values  if  frequencies  above  1000  cps  were  boosted.  Such  a  weighting 
would  affect  the  nasals  but  little,  because  they  have  very  little  energy  in  that  region; 
however,  it  would  greatly  improve  our  chances  of  separating  [i]  from  the  nasals.  In  the 
next  quarter  we  intend  to  investigate  this  further. 

An  attempt  was  also  made  to  establish  the  difference  between  the  two  nasal  conso¬ 
nants.  For  this  purpose  we  measured  spectra  of  these  sounds  (1).  No  consistent  dif¬ 
ference  could  be  found  between  spectra  of  [n]  and  of  [m].  See  Fig.  XI-2. 

We  pursued  this  question  further  and  attempted  to  synthesize  some  vowels  and  the 
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Spectra  in  the  first  row  were  obtained  from  the  syllables  /mi/  and  /ni/. 
Spectra  in  the  second  and  third  rows  are  of  the  nasal  and  the  vowel  in  the 
syllables  /mu/  and  /ni/,  respectively.  Note  the  predominance  of  the  low 
frequencies  in  all  spectra. 
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nasals  by  means  of  oscillators.  A  good  [m]  could  be  synthesized  by  a  single  oscillator 
tuned  at  150  cps.  When  a  second  oscillator  tuned  to  300  cps  was  added,  no  appreciable 
difference  was  noted  until  the  output  of  the  second  oscillator  was  substantially  equal  to 
that  of  the  first.  Under  these  conditions  [u]  was  perceived.  In  order  to  obtain  [n]  it 
was  necessary  to  introduce  a  very  weak  component  (-40  db  below  the  fundamental)  at 
about  2300  cps.  When  this  component  was  increased  in  intensity  to  -17  db  a  sound  much 
like  the  German  [ii]  was  heard.  The  best  [n]  was  synthesized  by  introducing  two  fairly 
strong  components  (-20  db)  at  600  cps  and  750  cps  in  addition  to  the  above-mentioned 
weak  component  at  2300  cps.  It  was,  however,  observed  that  the  2300  cps  component 
was  much  more  essential  in  the  perception  of  [n]  than  the  600  and  750  cps  components. 

M.  Halle,  G.  W.  Hughes 
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Dr.  C.  M.  Witcher 
L.  Washington,  Jr. 

A.  ACOUSTIC  OBSTACLE  DETECTION 

In  the  Quarterly  Progress  Report  of  April  15,  1954,  we  described  a  high-frequency 
sound  pulse  projector  incorporating  automatic  scanning.  The  results  obtained  with  this 
device  were  of  sufficient  interest  to  warrant  the  construction  of  a  second  model  incor¬ 
porating  several  improvements.  Figure  XII- 1  is  a  photograph  of  the  new  projector;  its 
circuit  is  shown  schematically  in  Fig.  XII-2.  This  model,  considerably  reduced  in  size, 
is  about  8  l/2  inches  long  and  3  5/8  inches  in  diameter.  It  weighs  slightly  less  than 
2  pounds.  As  in  the  earlier  model,  scanning  is  accomplished  by  means  of  a  deflecting 
ellipse,  but  the  motion  of  the  ellipse  in  the  new  model  takes  the  form  of  an  angular 
oscillation  rather  than  a  rotation.  This  motion  is  achieved  with  a  Scotch  crosshead 
which  is  driven  by  the  motor  near  the  bottom  of  the  device.  The  crosshead  carries  a 
toothed  rack  which  engages  a  small  spur  gear  on  the  shaft  of  the  ellipse.  This  scanning 
arrangement  is  considerably  superior  to  that  previously  used,  since  the  latter  allowed 


Fig.  XII-2 

Sound  pulse  projector  schematic. 


Fig.  XII- 1 

Sound  pulse  projector. 
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the  sound  beam  to  emerge  from  the  projector  for  only  about  half  of  each  scanning  cycle. 
With  the  oscillating  system,  the  beam  is  available  for  obstacle  detection  at  all  times. 

The  angular  oscillation  of  the  scanning  ellipse  is  almost  strictly  sinusoidal,  having  an 
angular  amplitude  of  60°  and  a  period  of  about  1. 0  sec.  The  reduced  dimensions  of  the 
new  model  result  in  one  disadvantage,  however;  the  sound  beam  is  more  divergent.  The 
motor  in  the  new  projector  is  of  the  same  general  type  as  that  used  in  the  older  model, 
but  it  was  made  specially  for  us  by  the  H.  D.  Brailsford  Company,  of  Rye,  New  York. 

The  operation  of  the  electrical  circuit  is  as  follows.  The  inductive  "kicks"  appearing 
across  the  0.  1 -henry  coil  at  each  break  of  the  motor  contacts  produce  shock  excitations 
of  the  tuned  circuit  associated  with  the  crystal  transducer  through  the  type  NE-2  neon 
lamp.  The  neon  lamp  functions  as  a  switch  that  disconnects  the  tuned  circuit  from  the 
exciting  coil  during  the  second  half-cycle  of  oscillation  after  the  arrival  of  an  exciting 
pulse,  so  that  the  tuned  circuit  is  then  left  free  to  ring  until  the  next  exciting  pulse 
arrives.  This  scheme  has  given  much  more  satisfactory  sound  pulses  than  were 
previously  obtained.  The  resonant  frequency  of  the  tuned  circuit  and  transducer  is 
9.  0  kc/sec,  and  its  Q  is  of  the  order  of  10.  The  amplitude  of  the  sound  pulses  produced 
by  the  new  projector  has  not  been  measured,  but  it  is  apparently  greater  than  that 
obtained  with  the  old  model  (104  db  above  a  sound  pressure  of  0.  0002  dyne/cm  at  dis¬ 
tance  of  18  inches  from  the  projector). 

We  have  described  this  projector  in  considerable  detail  because  it  appears  to  repre¬ 
sent  a  fairly  close  approach  to  what  can  be  ideally  expected  from  this  type  of  obstacle 
detector.  Numerous  tests  conducted  with  it  during  the  past  few  months  apparently  lead 
to  the  conclusion  that  it  cannot  be  expected  to  serve  as  a  genuinely  useful  travel  aid  to 
the  blind  except  under  highly  specialized  conditions.  All  obstacles  with  rough  or  irreg¬ 
ular  surfaces  yield  clearly  audible  sound  reflections  at  distances  up  to  15  or  20  ft,  but 
obstacles  such  as  small,  smooth,  round  pipes  or  columns,  as  well  as  smooth  walls  on 
which  the  sound  beam  is  incident  obliquely,  do  not  yield  reliably  detectable  sound  reflec¬ 
tions.  A  further  increase  of  perhaps  10-15  db  in  sound  intensity  might  partially  elimi¬ 
nate  these  difficulties,  but  then  the  intensity  would  be  so  great  as  to  render  the  user  of 
the  instrument  unduly  conspicuous.  We  plan  to  attempt  this,  however,  if  and  when  a 
more  rugged,  but  sufficiently  small,  transducer  can  be  found  or  developed. 

The  suggestion  has  been  made  that  a  similarly  constructed  ultrasonic  projector 
together  with  a  pair  of  ultrasonic  receivers  (to  provide  binaural  localization)  might 
provide  an  improved  obstacle  detection  system.  The  highly  satisfactory  scanning  mech¬ 
anism  developed  for  the  present  model,  together  with  its  reasonably  low  total  power 
consumption  (about  150  mw),  definitely  warrants  further  work  along  one  of  the  two  lines 
just  mentioned. 
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B.  NEW  STEP-DOWN  DETECTOR 

Some  months  ago  it  was  concluded  that  a  complete  redesign  of  the  step-down  detector 
was  needed  if  reliable  results  were  ever  to  be  obtained.  The  first  step,  a  redesign  of 
the  optical  system,  was  carried  through,  but  a  delay  of  several  months  followed  because 
of  difficulties  encountered  in  obtaining  the  necessary  components.  The  new  receiver  and 
light  source  have  now  been  completed  and  will  be  described  briefly. 

The  principal  objectives  underlying  the  new  design  have  been:  to  achieve  a  consider¬ 
able  improvement  in  signal -to-noise  ratio,  and  to  reduce  substantially  the  blocking 
effects  caused  by  sunlight-to-shadow  transitions,  as  discussed  in  the  Quarterly  Progress 
Report  of  January  15,  1954.  In  designing  the  optical  system,  three  steps  have  been  taken 
to  further  these  objectives:  (a)  Much  faster  optics  have  been  used  in  the  receiver,  per¬ 
mitting  a  considerable  reduction  in  photocell  area;  (b)  The  vibration  frequency  of  the 
light  source  has  been  increased  from  30  cps  to  160  cps;  (c)  By  use  of  a  specially 
designed  lamp  for  the  source,  together  with  a  Schmidt  optical  system,  the  beam  from 
the  source  has  been  shaped  to  coincide  closely  with  the  receiver  image  space  in  the 
region  where  the  latter  "looks  at"  the  terrain. 

The  optical  element  in  the  receiver  is  a  parabolic  mirror  on  a  lightweight  plastic 

base,  which  was  generously  donated  to  us  by  the  Polaroid  Corporation  of  Cambridge, 

Massachusetts.  This  mirror  has  a  focal  length  of  1.  06  inches  and  a  circle  of  confusion 

less  than  0.25  mm.  The  receiver  aperture  is  2.45  inches,  thus  giving  a  speed  of 

F-0.43.  The  PbS  photocell,  produced  for  us  by  Photoswitch  Corporation,  is  of  the 

evaporated  type,  enclosed  in  an  evacuated  bulb,  and  the  dimensions  of  its  active  area 

are  0.  005  by  0.  050  inch.  This  is  about  1  / 1 5  the  area  of  our  previous  cells,  and  hence, 

for  the  same  amount  of  incident  light  energy,  it  should  provide  an  improvement  in  signal  - 

1  /2 

to-noise  ratio  (assuming  same  signal  frequency)  of  15  '  ,  or  11.8  db. 

Since,  for  a  given  bandwidth,  the  noise  power  from  a  PbS  cell  varies  inversely  with 
frequency,  the  increase  of  5.  3 -fold  in  source  vibration  frequency  should  provide  an  addi¬ 
tional  improvement  of  7.25  db.  Thus  we  may  reasonably  expect  a  total  improvement  in 
signal-to-noise  ratio  of  about  19  db  from  the  two  factors  so  far  considered.  Measure¬ 
ments  made  by  the  Photoswitch  Corporation  on  the  new  cells  quantitatively  confirm  this 
improvement. 

The  increased  source  frequency,  combined  with  the  use  of  a  more  sharply  tuned 
amplifier,  must  certainly  provide  some  relief  from  the  blocking  effects  produced  by 
light-to-shadow  transitions.  For  a  walking  speed  of  5  ft/sec  and  for  a  receiver  image 
space  extending  along  the  surface  of  the  terrain  a  distance  of  1 .  5  inches  in  the  direction 
of  motion  of  the  detector,  the  rise  times  of  the  pulses  produced  by  light-to-shadow  tran¬ 
sitions  will  be  of  the  order  of  1.  5/60  =  25  msec,  so  that  most  of  the  energy  in  these 
pulses  will  be  at  frequencies  below  320  cps,  that  is,  twice  the  vibration  frequency  of  the 
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source.  The  use  of  a  Polaroid  type  XRX-40  infrared  filter  in  front  of  the  receiver  (short 
wavelength  cutoff  8600  angstroms)  further  serves  to  minimize  the  influence  of  sunlight. 

The  greatly  improved  light  concentration  from  the  source  should  result  in  a  further 
improvement  in  signal-to-noise  ratio,  although  this  improvement  has  been  sacrificed 
to  some  extent  by  the  use  of  a  lower  power  lamp  in  the  source.  The  lamp  now  used  was 
developed  especially  for  us  by  the  lamp  division  of  the  General  Electric  Company,  and 
we  are  greatly  indebted  to  Clifton  Found,  of  that  organization,  for  contributing  several 
of  these  lamps.  The  tungsten  filament  in  the  new  lamp  is  in  the  form  of  a  helix  only 
0.  005  inch  in  diameter  and  approximately  0.  05  inch  long.  It  is  designed  to  operate  at 
a  nominal  temperature  of  2400°K  with  an  applied  voltage  of  1 . 5  volts  and  a  current  of 
0.  25  amp.  These  lamps  have  also  proved  capable  of  withstanding  the  vibrational  accel¬ 
erations  of  the  order  of  250g  to  which  they  are  subjected  when  the  source  is  in  operation. 

The  receiver  is  housed  in  a  cylindrical  bakelite  case  2.  75  inches  in  diameter  and 
2.  25  inches  long;  its  total  weight  is  3.  5  oz.  The  photocell  mounting  embodies  provisions 
for  accurately  centering  the  cell  with  respect  to  the  axis  of  the  mirror,  and  a  special 
optical  alignment  tool  has  been  developed  for  this  purpose. 

A  photograph  of  the  source  assembly  is  shown  in  Fig.  XII-3.  As  may  be  seen,  the 
frame  consists  of  a  bakelite  box  with  aluminum  corner  supports  to  provide  added  rigid¬ 
ity.  Note  that  the  aluminum  ring  which  holds  the  Schmidt  corrector  plate  fits  snugly  into 
the  front  end  of  the  box  to  afford  still  more  rigidity  in  the  neighborhood  of  the  supports 

for  the  vibrating  element.  This  element  is  tuned 
to  mechanical  resonance  at  a  frequency  of  160  cps 
and  is  driven  from  a  small,  balanced  armature 
type  of  headphone  unit  mounted  in  the  collet 
assembly  on  top  of  the  box.  The  system  operates 
in  the  manner  of  an  ordinary  buzzer,  utilizing  a 
make-and-break  contact  arrangement  and  a  small 
transformer  powered  from  the  same  1.5  volt 
source  which  furnishes  power  for  the  lamp.  The 
position  of  the  lamp,  at  the  end  of  the  vibrating 
element  toward  the  mirror,  is  accurately  deter¬ 
mined  to  insure  that  the  filament  moves  in  an  arc 
lying  on  the  focal  surface  of  the  Schmidt  system. 
We  wish  to  express  our  sincere  thanks  to 
Harrison  J.  Merrill,  of  the  Evans  Signal  Labora¬ 
tory,  Fort  Monmouth,  New  Jersey,  for  providing 
us  with  the  optical  components  for  the  Schmidt 
system. 

Considerable  development  work  has  been 


Fig.  XII-3 

Step-down  detector  light  source. 
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required  to  obtain  a  vibratory  system  with  sufficient  rigidity  and  a  sufficiently  large 
mechanical  Q  to  give  a  vibration  amplitude  of  4-5°  at  160  cps  with  a  power  input  of  not 
more  than  100  mw.  These  conditions  have  now  been  achieved  satisfactorily. 

As  this  is  written,  the  complete  assembly  for  the  source  and  receiver  has  not  been 
finished;  we  have  therefore  not  had  an  opportunity  to  make  measurements  of  the  system 
performance. 
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Dr.  J.  S.  Barlow*  Dr.  D.  H.  Raab  L.  S.  Frishkopf 

Dr.  J.  R.  Hughes 

A.  A  PROBABILISTIC  APPROACH  TO  CERTAIN  NEUROELECTRIC  PHENOMENA 
1.  Basic  Data 

An  acoustic  stimulus  presented  at  the  ear  of  an  anesthetized  cat  evokes  an  electrical 
response  from  the  auditory  nerve  of  that  animal.  Our  primary  stimuli  to  date  have  been 
clicks  produced  by  feeding  rectangular  pulses  of  duration  0.  1  msec  into  a  PDR-10  ear¬ 
phone.  These  short  stimuli  evoke  spike-like  responses,  which  are  detected  by  a  fine 
wire  electrode  placed  within  the  bullar  cavity.  Suitable  electrode  placement  makes 
these  click  responses  substantially  free  of  microphonic  contamination  (1,  2).  A  typical 
response  is  shown  in  Fig.  XIII- 1. 

The  earliest  neural  component,  labeled  N^,  can  be  assumed  to  represent  the  syn¬ 
chronous  and  summated  firing  of  first  order  auditory  neurons.  The  peak-to-peak  ampli¬ 
tude  of  Nl  (App),  taken  as  a  measure  of  total  activity,  when  plotted  against  the  stimulus 
intensity  (S),  yields  a  monotonic  function  (Fig.  XIII-2)  referred  to  hereafter  as  the  inten¬ 
sity  function  (A  =  f(S))  (2).  Each  datum  point  represents  the  average  of  ten  or  more 
PP 

measurements . 

If  the  click  is  presented  in  a  background  of  continuous  white  noise  (limited  by  the 


Fig.  XIII- 1 

Typical  peripheral  response  to  click  (-30  db  re  1.29  volts)  with 
very  little  microphonic  content.  M,  microphonic;  N^,  first 
neural  component;  A  ,  peak-to-peak  amplitude. 


* 


From  the  Neurophysiological  Laboratory  of  the  Neurology  Service  of  the  Massachu¬ 
setts  General  Hospital. 


-87- 


Fig.  XIII- 2 

Intensity  function  obtained  from  Cat  C-307.  This  is  fairly 
typical  of  the  functions  obtained  from  animals  whose  thresh¬ 
olds  seem  to  be  normal. 


Fig.  XIII -3 

Masked  and  unmasked  intensity  functions. 
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Fig.  XIH-4 

Standard  deviation  of  response  amplitudes  to  a  given  stimulus  intensity 
plotted  as  a  function  of  stimulus  intensity.  Each  point  is  based  on  about 
100  responses.  Plotted  with  it  (dashed  line)  is  the  intensity  function 
obtained  from  the  same  data. 


(a) 


Fig.  XIII- 5 

(a)  A  possible  threshold  probability  distribution  function,  D(S). 
The  shaded  area  equals  the  probability  of  finding  the  threshold 
in  the  interval  to  S£.  The  total  area  under  the  curve  is  1. 

(b)  p(S),  the  probability  of  a  unit  responding  to  stimulus  S. 
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characteristics  of  the  system  that  delivers  the  acoustic  stimulus),  the  peak-to-peak 
amplitude  of  N^  is  reduced.  Plotting  this  amplitude  vs.  click  intensity  for  a  fixed  noise 
level  yields  (3)  a  modified,  or  masked,  intensity  function,  =  f(S,  N).  Several  such 

masked  functions  for  several  noise  levels,  together  with  the  corresponding  unmasked 
functions,  are  presented  in  Fig.  XIII-3.  Each  point  represents  the  average  of  the  peak- 
to-peak  amplitude  of  ten  or  more  responses  to  identical  stimuli. 

We  have  studied  the  variability  of  these  response  amplitudes  and  have  shown  (4)  that 
many  more  than  ten  responses  to  a  single  stimulus  intensity  were  needed  to  obtain  a 
statistically  meaningful  measure  of  variability.  From  a  sequence  of  100  responses  to 
a  given  stimulus  setting,  the  standard  deviation  of  the  peak-to-peak  amplitudes  was 
determined  for  a  number  of  stimulus  values.  The  resultant  function,  c  =  f(S),  is  plotted 
in  Fig.  XIII-4. 

2.  The  Model 

The  neural  unit  is  the  elemental  concept  of  our  model.  Its  relation  to  the  basic 
physiological  entities,  the  neuron,  the  sense -cell,  and  the  junction  between  them,  has 
been  discussed  in  previous  reports  (2,  3). 

The  neural  unit  has  the  following  properties: 

a.  It  obeys  the  all-or-none  law;  that  is,  it  is  characterized  by  only  two  states, 
maximal  response  or  no  response. 

b.  Its  threshold,  or  "triggering  level"  of  stimulation,  fluctuates  in  time.  No  attempt 
is  made  to  specify  the  detailed  time  course  of  this  fluctuation.  Rather,  the  description 
is  limited  to  a  specification  of  the  fraction  of  the  time  the  unit  spends  in  any  threshold 
interval,  i.  e.,  by  a  threshold  probability  distribution.  Fig.  XIII-5.  It  is,  in  addition, 
assumed  that  the  characteristic  time  of  fluctuations  is  very  small  compared  to  the  inter¬ 
val  (about  two  seconds)  between  stimulus  presentations. 

c.  The  firing  of  a  unit  alters  its  threshold  distribution,  which  then  becomes  a  func¬ 
tion  of  the  time  At  since  firing.  For  times  At  that  are  sufficiently  long  (greater  than 

0.  1  sec)  the  unit  returns  to  the  resting  distribution.  For  At's  that  are  less  than  2  msec 
the  unit  will  respond  to  no  stimulus,  whatever  its  intensity.  This  time  interval  is  called 
the  absolute  refractory  period  of  the  elemental  unit. 

Our  model  requires  a  second  concept:  a  population  of  neural  units;  that  is,  an 
aggregate  of  neural  units  whose  properties  a,  b,  and  c  involve  identical  parameters. 

The  present  model  assumes  that  at  least  two  populations  of  neural  units  are  contri¬ 
buting  to  the  peak-to-peak  amplitude  of  the  responses.  This  classification  is  based 
mainly  on  differences  in  mean  threshold  values,  with  each  unit  labeled  as  either  sensi¬ 
tive  or  insensitive.  It  may  eventually  be  found  necessary  to  subdivide  one  or  both  of 
these  populations  into  subpopulations. 

No  interaction  between  neural  units  has  been  introduced  explicitly  into  the  model. 
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Moreover,  the  threshold  fluctuations  occurring  in  any  two  units  are  assumed  to  be 
statistically  independent. 


3.  The  Model  and  the  Data;  Predictions 


A  single  neural  unit,  presented  with  a  stimulus  S,  will  respond  if  its  threshold  is 
less  than  or  equal  to  S  and  not  otherwise.  If  the  threshold  probability  distribution 
(Fig.  XIII-5(a))  is  represented  by  D(S),  then  the  probability  of  finding  the  unit  with 
threshold  less  than  S,  and  therefore  the  probability  that  the  unit  will  respond  to  a  stim¬ 
ulus  S,  is  given  by  (Fig.  XIII- 5(b)) 


S 


(1) 


D(s)  ds 


oo 

A  single  population  of  neural  units,  all  of  which  have  the  same  threshold  probability 
distribution  D(S),  gives  rise  to  an  intensity  function 


S 


D(s)  ds 


(2) 


■00 


where  r  is  the  magnitude  of  the  response  of  a  single  unit  and  N  is  the  number  of  units 
in  the  population.  If  a  single  population  were  involved,  the  function  D(S)  could  be 
inferred  at  once  by  differentiating  the  observed  intensity  function  and  making  use  of  the 
normalization 


(3) 


D(s)  ds  =  1 


If  we  are  dealing  with  two  independent  populations,  we  obtain 


where  the  subscripts  1  and  2  refer  to  the  two  populations.  The  situation  is  particularly 
simple  if  D^S)  and  D^S)  are  disjoint  (Fig.  XIII-6).  Then  again  differentiation  plus  the 
normalization  conditions 


permit  immediate  inference  of  D^(S)  and  of  D^S).  If  the  distributions  overlap. 
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Fig.  XIII- 6 

The  type  of  intensity  function  obtained  from  two 
disjoint,  or  nearly  disjoint,  populations. 


separation  in  the  region  of  overlap  on  the  basis  of  the  observed  intensity  function  is  no 
longer  possible . 

The  unmasked  intensity  functions  of  Figs.  XIII-2  and  XIII-3  are  of  the  type  shown  in 
Fig.  XIII-6  and  are,  therefore,  interpretable  within  a  framework  of  fluctuating  threshold 
and  two  populations.  Indeed,  it  was  data  of  this  sort  that  originally  suggested  the  two 
populations  hypothesis. 

We  now  apply  the  model  to  the  problem  of  masking.  A  population  characterized  by 
fluctuating  thresholds  and  a  period  of  nonresponsiveness  after  firing  could  be  rendered 
nearly  incapable  of  responding  to  a  stimulus,  such  as  a  click,  by  a  continuous  stimulus 
such  as  noise,  even  at  a  very  low  level.  This  conclusion  is  reached  in  the  following 
manner  (3):  neural  units  continually  drop  in  threshold  below  the  average  level  of  the 
masking  noise,  and  fire.  The  "location 11  of  the  noise  in  the  distribution  determines  the 
rate  at  which  such  firing  occurs.  The  weaker  the  noise,  the  less  often  a  unit  fires.  If 
threshold  fluctuation  is  sufficiently  rapid,  nearly  all  units  will  have  fired  in  response  to 
even  a  weak  noise  in  the  2  msec  preceding  the  click.  Hence  practically  all  units  will  be 
absolutely  refractory  at  the  time  the  click  is  presented,  and  none  will  respond  to  it. 

If  a  second  population  is  present  as  well,  whose  threshold  range  lies  far  above  the 
noise  level,  it  will  not  be  affected.  As  the  noise  level  is  further  raised,  the  number  of 
units  belonging  to  this  population  that  respond  to  the  click  will  also  decrease.  The 
data  of  Fig.  XIII-3  support  these  predictions.  They  show  that  a  weak  noise  eliminates 
the  sensitive  population  completely,  leaving  the  growth  of  the  insensitive  portion  of  the 
intensity  function  substantially  intact.  A  still  weaker  noise  permits  the  appearance  of 
a  partial  response  from  the  sensitive  population.  As  suggested  above,  the  noise  level 
at  which  the  masked  response  is,  for  example,  one  half  as  large  as  the  unmasked 
response,  depends  on  the  rate  of  threshold  fluctuation.  On  the  basis  of  a  simple  model 
of  fluctuations,  involving  random  transitions  between  threshold  states,  this  calculation 
has  been  carried  out  and  the  rate  determined  (3).  The  characteristic  time  of  occupancy 
of  a  threshold  state  is  found  from  these  data  to  lie  between  0.  25  msec  and  0.  5  msec. 

The  significant  features  of  the  masking  curves  are  thus  seen  to  follow  from  the 
assumptions  of  the  model. 
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Finally,  let  us  inquire  as  to  the  variability  in  response  amplitudes  expected  on  the 
basis  of  this  model.  At  a  given  level  of  stimulation,  S,  the  probability  of  a  unit  firing 
is  p(S).  The  probability  that  exactly  n  units  out  of  a  population  of  N  will  fire  is 

P(„)=(£V(1-p)N-n  (4) 

Plotted  as  a  function  of  n,  P(n)  is  peaked  at  n  =  Np;  for  sufficiently  large  N,  P(n)  may 
be  fitted  by  a  gaussian  distribution.  Such  a  fitting  procedure  yields  the  result  that 

o-  =  [Np  (1-p)]1/2  (5) 

If  we  introduce  now  an  explicit  dependence  on  S, 

o-(S)  =  [Np(S)  (1  -p(s))]l/2 

For  A  =  rn,  P  (A)  =  P(n),  where  P  (A)  is  the  probability  of  a  response  amplitude  A. 
Since  the  distributions  P(n)  and  P  (A)  are  related  by  a  scale  factor  r,  it  follows  that 
the  standard  deviation  of  response  amplitudes  to  a  stimulus  S  is  given  by 

o*(S)  =  ro-(S)  =  r  [Np (S )  (1  -  p(S))]l/2 
The  average,  and  most  probable,  response  is  given  by 


Aay(S)  =  rNp(S) 


(6) 


This  expression  is,  of  course,  the  intensity  function  (Eq.  2).  Figures  XIII-7(a)  and 
XIII-7(b)  show  an  intensity  function  for  one  population  and  the  corresponding  function 

sfc  ^ 

cr  (S).  Note  that  for  both  very  strong  and  very  weak  stimuli,  cr  (S)  =  0,  since  for  a  maxi¬ 
mal  stimulus  (p(S)  =  1)  all  units  respond  all  of  the  time,  while  for  a  minimal  stimulus 
(p(S)  =  0)  no  unit  will  ever  respond.  <r(S)  is  maximal  for  S  such  that  p(S)  =  l/2.  The 
relative  variability,  given  by 


a*(S) 


A  (S) 
av  ' 


N' 


1  -p(S)“ 


1/2  L  P(S) 


1/2 


(?) 


is  determined  primarily  by  the  population  size. 

If,  in  addition,  a  stimulus -independent  source  of  variability  is  present  (for  instance, 
amplifier  noise)  the  total  standard  deviation  of  the  response  is 


o"T(S) 


(a*(S))2  +  cr2 


nl/2 


A 


(8) 


Here  <r^  represents  all  additional  sources  of  variation  and  is  assumed  to  be  independent 
of  S.  Figure  XIII-7(c)  shows  the  function  cr^  for  a  given  value  of  <r^. 

If  two  more  or  less  disjoint  populations  are  involved,  one  would  expect  (^(S)  to  have 
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Fig.  XIII-7  Fig.  XIII- 8 

(a)  A  one -population  intensity  function  for  (a)  A  two-population  intensity  function, 

a  gaussian  threshold  probability  distribu-  (b)  A  corresponding  function  crrp(S)  of  the 

tion.  (b)  The  corresponding  function  cr*(S)  kind  predicted  by  the  model, 
predicted  from  the  model,  (c)  The  function 
crrp(S)  obtained  from  the  model  for  a  given 
value  of  a  .  . 


two  peaks.  The  peaks  occur  at  p^S)  =  l/2  and  at  p2(S)  =  l/2;  here  the  subscripts  1  and 
2  refer  to  the  two  populations.  Figure  XIII-8  shows  both  cr^  and  the  corresponding  inten¬ 
sity  function  for  the  two -population  case. 

If  we  compare  from  Fig.  XIII-8  with  the  variability  data  of  Fig.  XIII-4  there  is 
good  qualitative  agreement  with  the  prediction  of  the  model  in  the  range  of  the  sensitive 
population.  The  experimental  curve  cr(S)  exhibits  a  peak  at  a  point  where  on  the  basis 
of  the  intensity  function  we  conclude  that  approximately  50  per  cent  of  the  units  of  the 
sensitive  population  have  contributed.  cr(S)  for  the  insensitive  population  does  not  show 
the  expected  behavior.  This  may  be  due  to  our  inability  to  deliver  a  truly  maximal 
stimulus.  Such  a  stimulus  might,  of  course,  overload  the  mechanical  system  of  the 
cat's  ear. 

L.  S.  Frishkopf,  W.  A.  Rosenblith 
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Fig.  XIII- 9 

Block  diagram  of  analog  correlator. 


B.  ANALOG  CORRELATOR  FOR  ELECTROENCEPHALOGRAPHY 

With  the  completion  of  the  recycler  and  magnetic  delay  drum  (with  associated  elec¬ 
tronics)  the  analog  correlator  (Fig.  XIII-9)  has  successfully  produced  a  correlogram 
from  a  tape-recorded  EEG.  The  delay  drum  has  already  been  described  (1,  2).  The 
recycler  plays  through  the  desired  sample  of  tape,  rewinds  it  to  the  beginning  of  the 
sample  and  then  automatically  repeats  the  play  and  rewind  operations  until  the  correlo¬ 
gram  is  completed.  During  each  rewinding  operation,  the  recycler  advances  the  movable 
arm  of  the  drum  through  one  time  delay  unit  (At).  A  piece  of  silvered  tape  stuck  to  the 
back  of  the  magnetic  tape  at  both  ends  of  the  sample  to  be  correlated  reflects  light  into 
a  photocell  providing  signals  for  the  cycling  operations. 

A  technical  report  describing  the  detailed  functioning  of  the  complete  correlator  will 
be  published. 

J.  S.  Barlow,  M.  A.  B.  Brazier,  R.  M.  Brown 
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Prof.  E.  A.  Guillemin  N.  DeClaris  Y.  C.  Ho 

Prof.  R.  E.  Scott  T.  E.  Stern 


TIME  DOMAIN  SYNTHESIS  BY  DELAY  LINE  TECHNIQUES 


Research  in  the  past  quarter  has  resulted  in  the  development  of  an  experimental  type 
of  discrete  filter  and  its  application  to  time  domain  equalization  problems  in  both  open- 
and  closed-loop  systems. 

It  has  been  shown  (1,2)  that  a  discrete  filter  of  the  general  type 


N 


dVsT) 


ak€ 


,-skT 


k=0 


(1) 


can  be  used  to  approximate  the  transfer  function 


F(s) 


n 

X 

k=0 


b,  s 
k 


m 


k 

s 


(2) 


The  chief  advantage  of  discrete  filters  (3)  is  that  mathematically  they  can  be  handled 
easily  in  time  domain,  which  places  in  evidence  the  transient  response  of  a  system. 
Furthermore,  discrete  filters  are  more  versatile  than  their  conventional  counterparts, 
in  the  sense  that  their  transfer  characteristics  can  be  changed  at  will  after  their  con¬ 
struction  by  merely  changing  the  coefficients  a^  s  in  Eq.  1.  This  property  is  very  useful 
in  practical  time  domain  equalization  problems,  since  experimental  adjustment  over 
quite  a  wide  range  of  transfer  characteristics  is  provided. 

Preliminary  theoretical  study  shows  that  discrete  filters  can  be  used  to  equalize 
closed  loop  systems  with  good  results.  The  approach  used  is  in  the  time  domain  and 
is  based  on  the  impulse  response  of  the  system. 

An  experimental  discrete  filter  is  being  constructed  with  a  magnetic  recorder  having 
delayed  pick-up  heads  along  a  continuous  tape  loop.  The  unit  will  be  used  in  the  control 
loop  of  a  servomultiplier .  This  will  serve  not  only  as  an  experimental  check  but  also 
as  a  test  of  the  practicability  of  the  proposed  technique  for  time  domain  equalization. 

Y.  C.  Ho 
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output  signals;  thus  it  differs  from  the  type  described  in  sampled-data  systems, 
although  mathematically  they  are  almost  equivalent. 
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XV.  NETWORK  SYNTHESIS 


Prof.  E.  A.  Guillemin  P.  M.  Lewis  II  R.  A.  Pucel 

Prof.  F.  M.  Reza  M.  Strieby 


VOLTAGE  TRANSFER  SYNTHESIS  -  RLC  LATTICE 


The  work  on  voltage  transfer  synthesis  has  been  extended  to  cover  the  RLC  lattice. 
(The  RC  lattice  was  discussed  in  the  Quarterly  Progress  Report,  April  15,  1954.)  The 
necessary  and  sufficient  conditions  for  a  transfer  function  to  be  realizable  as  a  lattice 
have  been  found,  and  methods  of  synthesis  have  been  developed. 

1 .  Theorem 

The  necessary  and  sufficient  conditions  that  a  voltage  transfer  function 
K(sn  +  a  .  Sn_1  +  .  .  .  +a,  S  +  an) 

V  n-1  1  0)  KN 

i  D 

Sm  +  b  .  Sm  +  .  .  .  +  b,  S  +  bn 
m-1  1  0 


be  realizable  as  an  RLC  lattice  are:  1.  D  is  a  Hurwitz  polynomial;  2.  D  +  KN  and 
D  -  KN  are  Hurwitz  polynomials.  The  order  of  these  two  polynomials  can  differ  by  no 
more  than  two.  (This  condition  specifies  the  maximum  value  of  the  parameter  K.) 
Proof: 

(a)  The  proof  of  condition  1  is  well  known  and  will  not  be  repeated  here. 

(b)  Consider  the  lattice  shown  in  Fig.  XV- 1. 

The  voltage  transfer  function  is 


A  = 


KN 

D 


zh  -  z  1  -  <W 

b  a  _ _ 

Zb  +  Za  1  +  (Za/Zb) 


Solving  this  equation 


Za  _  1  -  A  _  D  -  KN 

Z,  1  +  A  D  +  KN 
b 


Since  the  ratio  of  two  positive  real  impedances  must  have  Hurwitz  polynomials  in 
both  its  numerator  and  denominator,  the  theorem  is  proved. 


Fig.  XV- 1 
An  RLC  lattice. 
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2.  Synthesis 

Given  a  voltage  transfer  function  which  satisfies  the  above  conditions,  one  can 
immediately  write 


Za  1  -  A  D  -  KN 

Z,  1  +  A  D  +  KN 

b 


The  problem  is  then  to  split  Z^/Z^  into  Z&  and  Z^.  Three  different  ways  of  doing  this 
are  suggested. 

1.  Enough  surplus  factors  can  always  be  added  so  that  Z  /  Z^  can  be  immediately 
split  into  Z&  and  Z^.  However,  after  this  is  done,  the  methods  of  Brune  or  Bott  and 
Duffin  must  be  used  to  synthesize  the  resulting  networks.  This  direct  approach  may 
therefore  be  impractical. 

2.  If  either  the  odd  or  the  even  parts  of  D  +  KN  and  D  -  KN  are  proportional,  a 
particularly  easy  synthesis  is  possible.  The  details  will  be  presented  for  the  odd  parts 
proportional.  The  extension  to  the  other  case  is  obvious. 

If 


D  -  KN  ml  +  nl 
D  +  KN  m^  +  kn^ 


where  m  and  n  are  the  even  and  odd  parts  of  the  polynomials,  then  we  can  write 


1  +  (m1/n1 ) 
k  +  (m2/n1 ) 


Z 

a 


m. 


1  + 


n . 


and 


=  k 


m 

H - 

n 


2 

1 


Then  Z  and  Z,  are  realized  as  lossless  networks  in  series  with  resistors.  Part 
a  b 

of  this  resistance  can  be  removed  from  the  lattice  to  give  a  series  resistor.  If  a  paral 
lei  resistance  termination  is  desired,  the  function  can  be  split  into 

1  1 

Z  = -  and  Z,  =  - 

k  +  (n^/nj)  1  +  (mj/n^ 


This  type  of  synthesis  will  be  possible  when  N  is  an  even  or  odd  polynomial,  when 
D  is  of  the  second  order,  and  in  certain  other  particular  cases. 

3.  If  the  maximum  possible  gain  is  not  required,  a  synthesis  in  terms  of  partial 
fractions  is  possible.  Split  Z&/Z^  into 
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Z  =  1  -  A  and  Z,  =  1  +  A 
a  b 

(or  Za=TTA  and  Zb  =  T^a) 

Expand  A  in  partial  fractions.  The  resulting  terms  may  or  may  not  be  realizable 
as  simple  RLC  circuits.  However,  those  terms  that  are  not,  can  be  made  realizable 
by  the  addition  of  a  suitable  resistance.  Consider  the  1  in  the  expressions  1  -  A  and 
1  +  A  as  a  one-ohm  resistance  reservoir  that  can  be  used  for  this  purpose.  Pick  K 
small  enough  so  that  the  one  ohm  is  sufficient  to  make  all  the  terms  realizable.  The 
resulting  RLC  partial  fraction  canonic  forms  are  well  known. 

If  K  is  made  smaller  than  is  required,  some  extra  resistance  will  be  left  over  to 
remove  as  a  series  or  parallel  resistor.  Incidentally,  the  input  impedance  of  this  lattice 
is 

Z.  =  \  (Z  +  Zh)  =  +  (1+A)  =  i  ohm 

m  2  v  a  b7  2 

This  is  particularly  convenient  for  some  applications. 

P.  M.  Lewis  II 
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Prof.  S.  H.  Caldwell  Prof.  W.  K.  Linvill 

Prof.  D.  A.  Huffman  E.  J.  McCluskey,  Jr. 

There  is  a  general  algebraic  method  that  permits  the  design  of  switching  circuits 
to  perform  according  to  given  specifications  (1).  For  any  reasonable  set  of  specifica¬ 
tions,  it  is  possible  to  obtain  several  appropriate  switching  circuits. 

Often,  the  design  object  is  to  obtain  the  circuit  which,  while  satisfying  the  specifi¬ 
cations,  uses  the  fewest  number  of  switch  contacts  or  relay  springs.  The  algebraic 
design  method  is  a  great  help  toward  realizing  this  object;  but,  since  this  method  always 
results  in  a  series-parallel  type  of  network,  it  does  not  take  advantage  of  possible  con¬ 
tact  economies  that  nonseries -parallel  (bridge  or  nonplanar)  networks  might  afford. 
This  report  describes  a  preliminary  investigation  of  switching  circuits  of  the  bridge 
type. 

DEFINITIONS 

(1)  A  bridge  circuit  is  defined  as  a  circuit  that  contains  at  least  one  bilateral  branch. 

(2)  A  bilateral  branch  is  defined  as  a  branch  that  contains  at  least  one  switch  contact 
and  is  traversed  in  both  directions  in  tracing  out  all  possible  paths  from  input  to  output 
terminals  of  the  circuit  containing  the  branch.  A  path  is  considered  to  exist  even  if  the 
occurrence  of  two  complementary  transmissions  in  series  prevents  signal  flow  under 
normal  conditions.  In  other  words,  the  conditions  for  a  branch  to  be  bilateral  depend 
only  on  the  topological  properties  of  the  associated  network  and  as  such  can  be  deter¬ 
mined  from  the  graph  (2)  of  the  network. 

(3)  A  bilateral  transmission  is  the  transmission  function  associated  with  a  bilateral 
branch. 

(4)  A  bridge  transmission  is  any  transmission  corresponding  to  a  bridge  circuit. 

(5)  A  simple  bridge  circuit  is  a  bridge  circuit  in  which  neither  the  bilateral  trans¬ 
mission  nor  any  of  its  components  appear  elsewhere  in  the  circuit. 

Theorem  1.  Any  bridge  circuit  may  be  realized  in  the  form  shown  in  Fig.  XVI-1, 
where  B  is  the  bilateral  transmission  of  one  of  the  bilateral  branches. 


—  '  1  — 

T 

B 

~  T3  - 

-  A  — 

~r~ 

B 

—  D  - 

o - 

■ -  o  o  < 

1 

-  T  2  — 

— 1— 

-  L - 

2  1 

-  C  — 

— 1— 

—  E  - 

T  - -  - -  T 


Fig.  XVI-1 

General  form  of  a  bridge  circuit. 
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Fig.  XVI -2 

General  form  of  a  simple  bridge  circuit 
having  bilateral  transmission  B. 
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Proof:  Any  bridge  circuit  must  contain  at  least  one  bilateral  branch.  Let  the 
associated  bilateral  transmission  be  denoted  by  B,  and  the  branch  terminals  by  3  and 
4.  Denote  the  input  and  output  terminals  of  the  circuit  by  1  and  2,  respectively.  By 
definition  there  must  be  a  path  from  1  to  2,  traversing  B  from  3  to  4;  thus  there  must 
be  a  path  from  1  to  3  and  from  4  to  2.  Denote  the  corresponding  transmissions  by  Tj 
and  T  . .  Similarly  a  path  from  1  to  2  traversing  B  from  4  to  3  must  exist,  and  hence 
there  must  also  be  transmissions  T^  from  1  to  4,  and  T^  from  3  to  2. 

Corollary  1.  The  transmission  of  any  bridge  circuit  containing  a  bilateral  trans¬ 
mission  B  can  be  written  in  the  form: 

t  =  b(t1t4+t2t3)  +  t1t3+t2t4 

where  TjT^  ^  0  or  ^2^4  ^  and  BT2T3  ^  0  or  BTiT4  ^ 

(If  TjT^  =  T^T4  =  0,  then  T  =  B(TjT4  +  T2T^),  and  T  may  be  realized  as  economi¬ 
cally  in  series-parallel  form  (without  the  need  for  a  bilateral  transmission  B)  as  in 
bridge  circuit  form.  Therefore,  the  requirement  that  either  T^T^  or  T^T^  ^  0  will 
henceforth  be  assumed.  If  BT^T^  =  0  and  BT^T^  =  0,  B  is  unnecessary;  therefore  the 
requirement  that  BT2T^  ^  0  or  BT^T^  ^  0  will  also  be  assumed.) 

Theorem  2.  If  T  is  the  transmission  of  a  simple  bridge  circuit  with  bilateral  trans¬ 
mission  B,  T  can  be  written  in  the  form 

T  =  AD  +  CE  +  B(AE  +  CD) 

where  BAE  £  0,  or  BCD  ^  0,  AD  or  CE  ^  0,  and  A,  C,D,E  are  not  functions  of  B. 

The  corresponding  circuit  is  shown  in  Fig.  XVI- 2. 

We  now  proceed  to  determine  under  what  conditions  a  transmission  T  =  f  +  Bg  can 
be  realized  as  a  simple  bridge  circuit  with  bilateral  transmission  B. 

Expanding  A,  C,  D,  E  in  a  standard  sum: 

A  =  Art  X'X' .  .  .X'  +  A.  XIX'.  .  .X  +  ... 

012  n  112  n 

C  =  Cn  X'X'  .  .  .X'  +  C.  X'X'  .  .  .X  +  .  .  . 

012  n  112  n 

D  =  D„  X'X'  .  .  .X'  +  D.  X'X'  .  .  .X  +  .  .  . 

012  n  112  n 

E  =  E0  XI  XL.  .  .X'  +  E1X'X'...X  +  ... 

012  n  112  n 

AD  =  AnDn  X'X'  .X'  +  A.D.  X'X»  .  .X  +  .  .  . 

0012  n  1112  n 

CE  =  C„En  X'X'  .X'  +  C.E.  X’XL.  .  .X  +  ... 

0012  n  1112  n 

AD  +  CE  =  (A-D-  +  CnE  J  X'X'  .  .  .X'  +  (A  .D .  +  C  ,E  . )  X'X'  .  .  .  X  +  .  .  . 
v00  0  0/  12  n'll  ll7  12  n 
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f  =  f_  X'X'  .  .X'  +  f.  X'X'  .  .X  +  .  .  . 

012  n  112  n 

g  =  g„  X'X'  .  .X'  +  g.  X'X'  .  .X  +  .  .  . 

6  s0  1  2  n  &1  1  2  n 

It  is  assumed  that  T  has  been  simplified  by  eliminating  any  terms  which  appear 

in  both  f  and  g  through  use  of  the  theorem  XjX^.  •  .X^  +  (X^X^.  .  .  X^jB  =  X^X^.  .  .  X^. 

(The  notation  X^  is  used  where  either  X^  or  X^  may  apply.)  Thus  the  situation  that 

f.  =  1,  and  g.  =  1  will  not  occur.  In  other  words,  we  have  either  f.  =  g!  or  f.  =  0  and 
i  °i  i  toi  i 

g^  =  0.  However,  whenever  we  have  L  =  1,  g^  may  also  be  set  equal  to  1  in  order  to 
simplify  the  resulting  circuit  without  changing  T.  Thus  it  is  sufficient  to  consider  the 
three  possibilities:  L  =  0,  g^  =  0;  L  =  0,  g^  =  1;  and  f.  =  0,  g^  =  <}>,  where  <j)  signifies  that 
either  a  0  or  1  may  be  assumed  to  apply. 

Now  we  assume 

T  =  f  +  Bg  =  AD  +  CE  +  B(AE  +  CD) 
therefore 


f  =  AD  +  CE,  g  =  AE  +  CD 

/  j{c  >{<  \ 

Equating  coefficients  of  like  products  (X^X^.  •  . X  j  yields 


A.D.  +  C.E.  =  f. 
ii  ii  l 


A.E .  +  C.D.  =  g. 

ii  i  i  bi 

(1)  Assume  L  =  g^  =  0 


A.D.  +  C.E.  =  0  either  A.  =  C.  =  0 
1111  ii 


or  D.  =  E .  =  0 
l  i 


A.E.  +  C.D.  =  0 

ii  ii 


(2)  Assume  L  =  1,  g^  =  0 


A.D.  +  C.E.  =  1  either  A.  =  D.  =  1 
1111  ii 


or  C.  =  E.  =  1 

l  l 


A.E.  +  C  .D .  =  4> 

ii  ii 


(3)  Assume  L  =  0,  g_.  =  1 


A.D.  +  C.E.  =  0  either  A.  =  E.  =  1,  C.  =  D.  =  0 
1111  1111 


A.E.  +  C.D.  =  1 

ii  ii 


or  C.  =  D.  =  1,  A.  =  E.  =  0 
1111 


Consider  the  requirement  BAE  ^  0  or  BCD  0. 


AE  =  A0Eq  X'jX!,. 


.X'  +  A.E,  X'X'  .  .X  + 
n  1112  n 


This  requires  AE  ^  0  or  CD  4  0.  But 
So  that,  if  AE  £  0,  then  A^  =  E^  =  1 
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for  at  least  one  i.  Similar  reasoning  for  BCD  ^  0  leads  to  the  requirement  that 

A.  =  E.  =  1  orC.  =  D.  =  1  for  at  least  one  i. 
ii  ii 

Now  consider  the  requirement  that  AD  or  CE  ^  0.  This  in  turn  leads  to  the  require¬ 
ment  that  A.  =  D .  =  1  orC.  =  E .  =  1  for  at  least  one  i. 

ii  ii 

In  summary,  we  require  that 
(i)  A.  =  E.  =  1  or  C.  =  D.  =  i] 

'  ii  ii 


(ii) 


> 

A.  =  D.  =  1  or  C.  =  E.  =  1 

ii  ii! 


for  at  least  one  i 


Study  of  cases  1,  2,  and  3  shows  that  (i)  requires  that  g^  ^  0,  and  (ii)  requires  that  L  ^  0. 


Theorem  3.  If  T  =  f  +  gB  where  f  and  g  are  not  functions  of  B  and  f  ^  0  and  g  ^  0, 
then  T  may  be  realized  as  a  simple  bridge  circuit  which  contains  B  only  as  the  bilateral 
transmission. 


Theorem  4.  If  a  transmission  T  is  expanded  in  a  standard  sum: 

T  =  Tn  X'X'  .  .X'  +  T.  X'X’  .  .X  +... 

012  n  112  n 

and  if  whenever  T^  =  1  for  i  even  (odd)  T^  +  1  =  1  (T^  -1  =  1)  also  (for  at  least  one  such 
i),  and  if  T^  =  1  for  i  odd  (even)  with  T^  -  1  =  0  (T.  +  1  =  0)  (for  at  least  one  such  i),  then 
T  can  be  realized  as  a  simple  bridge  circuit  having  X^(X^)  as  the  bilateral  transmission. 
Example  1. 

T  =  xy  +  x'w  +  z(w+y) 

If  we  let  f  =  xy  +  x'w,  g  =  w  +  y,  it  is  clear  that  T  can  be  realized  as  a  simple  bridge 
circuit  with  z  bilateral.  For  the  order  wxy,  L  =  1  for  i  =  3,  4,  5,  7  and  g^  =  1  for  1,  6. 
Thus  for  i  =  1,  6 


r-H 

II 

w 

II 

!— H 

<< 

ci 

=  D 

A6  “  E6  ” 

C6 

=  D 

for  i  =  3,  4,  5,  7 

A^  =  D^  =  1  and  CL  =  <j>,  E^  =  4>  or  A^  =  D^  =  4>,  and  CL  =  E^  =  1 


for  i  =  0,  2 

A.  =  C.  =  0  and  D.  =  <b,  E.  =  cborD.  =  E.  =  0,  A.  =  <b,  C .  =  <b 

Figure  XVI-3  shows  how  the  A,  C,D,E  functions  may  be  simultaneously  minimized 
by  the  use  of  Karnaugh  charts  in  accordance  with  the  above  constraints.  The  first  row 
of  tables  shows  the  assignment  for  the  =  1  case;  the  second  row,  the  assignment  for 
the  g.  =  f  =  0  case;  the  third  row,  the  assignment  for  the  f.  =  1  case;  and  the  final  row. 
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Fig.  XVI-4 

Bridge  circuit  realization  of  transmis¬ 
sion  T  =  xy  +  x'w  +  z(w+y). 


Fig.  XVI -3 

Karnaugh  charts  used  to  minimize  A,  C, 
D,  and  E,  simultaneously. 


the  selection  of  the  <j>  condition  used  in  order  to  complete  the  specification. 

The  resulting  bridge  circuit  is  shown  in  Fig.  XVI-4. 

Another  method  by  which  this  circuit  could  have  been  realized  will  be  illustrated. 
By  Theorem  2,  if  z  is  the  bilateral  transmission  of  a  simple  bridge  circuit  having 
transmission  T  then  T  can  be  written:  T  =  AD  +  CE  +  z(AE  +  CD).  We  now  expand 
the  transmission  T  being  considered  about  z,  simplify,  and  equate  the  two  expressions. 
Thus  we  have 

T  =  yx  +  wx'  +  z(y+w)  =  AD  +  CE  +  B(AE  +  CD) 

where  B  =  z,  and  yx  +  x'w  =  AD  +  CE .  Set 
yx  =  AD;  y  =  A,  x  =  D 
wx'  -  CE;  w  =  C,  x'  =  E 

Now  form  AE  +  CD  =  yx'  +  wx.  This  clearly  is  not  equal  to  w  +  y,  but  when  the  complete 
expression  for  T  is  formed 

T  =  xy  +  x'w  +  z(xw  +  x'y) 

and  it  is  noted  that 

w(x'  +  xz)  =  w(x'  +  z)  and  y(x  +  x'z)  =  y(x+z) 
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it  is  evident  that  T  =  xy  +  x'w  +  z(xw  +  x'y)  does  equal  xy  +  x'w  +  z(w+y),  as  required. 
This  assignment  of  A,  B,  C,  and  D  leads  to  the  circuit  shown  in  Fig.  XVI-4.  Another 
possible  assignment  would  have  been 

A  =  X,  D  =  y,  E  =  w,  C  =  x' 

This  leads  to  the  circuit  of  Fig.  XVI-4  with  terminals  1  and  2  interchanged. 

E.  J.  McCluskey,  Jr. 
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